[BioC] extract introns

Steve Lianoglou mailinglist.honeypot at gmail.com
Fri Nov 11 04:19:10 CET 2011


Also, just to make sure the i's are dotted, and t's are crossed.

This code I gave:
R> gr.exons <- with(xcripts, {
 GRanges(chromosome_name,
   IRanges(exon_chrom_start, exon_chrom_end),
   ifelse(strand == 1, '+', '-'))
})

Assumes you've already loaded the GenomicRanges library via a call like so:

R> library(GenomicRanges)

You mentioned you're still a bit green with bioconductor, so I just
wanted to make sure you knew where to find "GRanges" ... sorry if
that's too obvious.

-steve

-- 
Steve Lianoglou
Graduate Student: Computational Systems Biology
 | Memorial Sloan-Kettering Cancer Center
 | Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact



More information about the Bioconductor mailing list