[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