Skip to content

Instantly share code, notes, and snippets.

View jwaageSnippets's full-sized avatar

jwaageSnippets

View GitHub Profile
GTF2TranscriptDB <- function(gtf.file, out.file = NULL, verbose = TRUE)
{
require(rtracklayer)
require(GenomicRanges)
require(GenomicFeatures)
min.info <- c("gene_id", "transcript_id", "exon_number")
if (verbose) message("Importing ", gtf.file)
gtf <- import.gff(gtf.file, asRangedData = FALSE)