[Bioc-devel] chromosome lengths (seqinfo) for supported BSgenome builds into GenomeInfoDb?

Kasper Daniel Hansen kasperdanielhansen at gmail.com
Mon Jun 15 15:17:38 CEST 2015


Hector: we still need genome version in whatever call we make.

I assume we will get an update when this has been implemented; I am itching
to use it.

Kasper

On Fri, Jun 5, 2015 at 5:04 PM, Michael Lawrence <lawrence.michael at gene.com>
wrote:

> That's already possible, basically:
>
> seqinfo(gr) <- seqinfo(Mus.musculus)
>
> Anyway, I'm with Tim's last suggestion. Just support both. Have an
> argument to genome<- like standardize=FALSE for low-level
> manipulation.
>
>
> On Fri, Jun 5, 2015 at 1:56 PM, Hector Corrada Bravo <hcorrada at gmail.com>
> wrote:
> > Sorry if I'm too late for this...
> >
> > I think having "hg19" in 'gr <- standardizeSeqinfo(gr, "hg19")' be
> > responsible for so much, can make this new function hard to maintain.
> > Might it be better if that argument takes an object corresponding to a
> > canonical annotation? E.g.,
> >
> >
> > library(Mus.musculus)
> > gr <- standardizeSeqinfo(gr, Mus.musculus)
> >
> >
> > On Fri, Jun 5, 2015 at 4:49 PM, Hervé Pagès <hpages at fredhutch.org>
> wrote:
> >>
> >> On 06/05/2015 01:48 PM, Gabe Becker wrote:
> >>>
> >>> I dunno, standardizeSeqInfo just seems really long for a function name
> >>> users are going to have to call.
> >>>
> >>> At the risk of annoying Herve further, what about
> >>>
> >>>     gr <- castSeqInfo(gr, "gh19")
> >>
> >>
> >> grrrrr!
> >>
> >>>
> >>> ?
> >>>
> >>> ~G
> >>>
> >>> On Fri, Jun 5, 2015 at 1:46 PM, Tim Triche, Jr. <tim.triche at gmail.com
> >>> <mailto:tim.triche at gmail.com>> wrote:
> >>>
> >>>     maybe standardizeSeqinfo or fixSeqinfo is clearer after all
> >>>
> >>>     Statistics is the grammar of science.
> >>>     Karl Pearson <http://en.wikipedia.org/wiki/The_Grammar_of_Science>
> >>>
> >>>     On Fri, Jun 5, 2015 at 1:41 PM, Gabe Becker <becker.gabe at gene.com
> >>>     <mailto:becker.gabe at gene.com>> wrote:
> >>>
> >>>
> >>>
> >>>         On Fri, Jun 5, 2015 at 1:39 PM, Tim Triche, Jr.
> >>>         <tim.triche at gmail.com <mailto:tim.triche at gmail.com>> wrote:
> >>>
> >>>             how about just
> >>>
> >>>             gr <- addSeqinfo(gr, "hg19")
> >>>
> >>>
> >>>         Add sounds like it's, well, adding rather than replacing (Which
> >>>         it sometimes would do.
> >>>
> >>>         gr <- fixSeqInfo(gr, "hg19")
> >>>
> >>>         instead?
> >>>
> >>>         ~G
> >>>
> >>>
> >>>         --
> >>>         Gabriel Becker, Ph.D
> >>>         Computational Biologist
> >>>         Genentech Research
> >>>
> >>>
> >>>
> >>>
> >>>
> >>> --
> >>> Gabriel Becker, Ph.D
> >>> Computational Biologist
> >>> Genentech Research
> >>
> >>
> >> --
> >> Hervé Pagès
> >>
> >> Program in Computational Biology
> >> Division of Public Health Sciences
> >> Fred Hutchinson Cancer Research Center
> >> 1100 Fairview Ave. N, M1-B514
> >> P.O. Box 19024
> >> Seattle, WA 98109-1024
> >>
> >> E-mail: hpages at fredhutch.org
> >> Phone:  (206) 667-5791
> >> Fax:    (206) 667-1319
> >>
> >> _______________________________________________
> >> Bioc-devel at r-project.org mailing list
> >> https://stat.ethz.ch/mailman/listinfo/bioc-devel
> >
> >
>
> _______________________________________________
> Bioc-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/bioc-devel
>

	[[alternative HTML version deleted]]



More information about the Bioc-devel mailing list