Followup to both of youse:

1) I've been using import() a fur piece lately and have found things in it
that don't make sense to me -- e.g. Brad Bernstein's recent chromatin marks
paper -- some of the segmented BED files load, some don't.  Clearly I need
to read the spec better, or something, but a bit of fussing with 'cat' and
'cut' "solved" the problem and let me load them into the DB.

2) Steve -- your bits worked like a charm -- I ended up adding generics for
getSeq(x='GRanges') and getSnpBases(x='GRanges') since rsidsToGRanges is so
handy.  This is particularly handy since I can now scan enhancers for SNP
variant potential impact and scan TSSs and putative cryptic TSS or SS
regions for motifs easily.  Using installed.genomes() and strsplit, the
first thing the function does is to see if the appropriate genome is
installed, and if so, load it.  If that fails, it calls available.genomes()
to see if the genome is available, and if so, tells the user to install it.
 If neither option works it complains that the genome is unknown and throws
an error.  For the time being I'm putting this and several other
ball-of-mud type hacks into a small package, but eventually it might be
nice for the generics to go into base?

thanks all


On Thu, Feb 9, 2012 at 2:16 PM, Steve Lianoglou <
mailinglist.honeypot@gmail.com> wrote:

> Hey Michael,
>
> I'll see if I can take a crack at it this w/e.
>
> Also I'm in complete agreement w.r.t i/o stuff and your thoroughness
> w/ rtracklayer is very much appreciated.
>
> If we do have a centralized track i/o package that is outside of
> rtracklayer, I nominate we name it bIO ... or something ;-)
>
> ... yes, I can hear the crickets ...
>
> -steve
>
> On Thu, Feb 9, 2012 at 4:48 PM, Michael Lawrence
> <lawrence.michael@gene.com> wrote:
> > Hi Steve,
> >
> > I had thought someone would want something like this, and it would be
> nice
> > to able to parse a GTF file into something more useful than a GRanges, as
> > rtracklayer supports now.
> >
> > So I'd definitely encourage you to go ahead, while recommending that it
> be
> > based on the import.gtf function in rtracklayer.
> >
> > We really need to have a centralized track I/O package. I'm not sure why,
> > but everyone seems intent on writing the same parser over and over again.
> > Duplicated effort frustrates me. It doesn't have to be rtracklayer, but
> > somewhere in the infrastructure we need reliable I/O that handles all of
> the
> > corner cases. I've embarked on a time consuming effort of implementing a
> > comprehensive test suite for rtracklayer, which should make it at least a
> > base for such a package.
> >
> > Thanks a lot for volunteering this contribution,
> > Michael
> >
> > On Thu, Feb 9, 2012 at 7:55 AM, Steve Lianoglou
> > <mailinglist.honeypot@gmail.com> wrote:
> >>
> >> Hi,
> >>
> >> I thought it would be handy to make a GenomicFeatures::TranscriptDb
> >> from a gtf file and was somehow surprised to see that I couldn't find
> >> such a function in GenomicFeatures.
> >>
> >> I'm happy to whip up such a function, but wanted to check in to make
> >> sure I'm not missing something like (1) you can already do it; or (2)
> >> it's not as straightforward as I think; (3) maybe it's there and I'm
> >> not looking hard enough.
> >>
> >> Right now I just want to build it on the gff that the knew versions of
> >> tophat build when you pass in a gtf/gff file of known transcripts (the
> >> --G/--GTF cmd line arg), but I think it'd be handy overall.
> >>
> >> I don't think gff/gtf's have any column for exon ordering, though, in
> >> which case I'd just assume the exons are ordered linearly (bye bye
> >> trans-splicing)).
> >>
> >> Good idea? Bad idea? Already done?
> >>
> >> -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
> >>
> >> _______________________________________________
> >> Bioconductor mailing list
> >> Bioconductor@r-project.org
> >> https://stat.ethz.ch/mailman/listinfo/bioconductor
> >> Search the archives:
> >> http://news.gmane.org/gmane.science.biology.informatics.conductor
> >
> >
>
>
>
> --
> 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
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor@r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>



-- 
*A model is a lie that helps you see the truth.*
*
*
Howard Skipper<http://cancerres.aacrjournals.org/content/31/9/1173.full.pdf>

	[[alternative HTML version deleted]]

