[Bioc-devel] coerce ExpressionSet to SummarizedExperiment

Hector Corrada Bravo hcorrada at gmail.com
Sat Sep 20 23:45:12 CEST 2014


Would be very useful in epivizr. Right now we have a little bit of code
that can be made much more general:
https://github.com/epiviz/epivizr/blob/master/R/register-methods.R#L73


On Sat, Sep 20, 2014 at 5:19 PM, Michael Love <michaelisaiahlove at gmail.com>
wrote:

> On Sep 20, 2014 2:15 PM, "Martin Morgan" <mtmorgan at fhcrc.org> wrote:
> >
> > On 09/20/2014 10:43 AM, Sean Davis wrote:
> >>
> >> Hi, Vince.
> >>
> >> Looks like a good start.  I'd probably pull all the assays from
> >> ExpressionSet into SummarizedExperiment as the default, avoiding data
> >> coercion methods that are unnecessarily lossy.  Also, as it stands, the
> >> assayname argument is not used anyway?
> >
> >
> > I think there will be some resistance to uniting the 'Biobase' and
> 'IRanges' realms under 'GenomicRanges'; considerable effort has gone in to
> making a rational hierarchy of package dependencies [perhaps Herve will
> point to some of his ASCII art on the subject].
> >
> > I have some recollection of (recent) discussion related to this topic in
> the DESeq2 realm, but am drawing a blank; presumably Michael or Wolfgang or
> ... will chime in.
> >
>
> Andrzej was working on a  conversion function for  DESeqDataSet <=>
> DGEList. I don't think it was for SummarizedExperiment and eSet.
>
> Mike
>
> > Martin
> >
> >
> >>
> >> Sean
> >>
> >>
> >> On Sat, Sep 20, 2014 at 10:38 AM, Vincent Carey <
> stvjc at channing.harvard.edu>
> >> wrote:
> >>
> >>> do we have a facility for this?
> >>>
> >>> if not, we have
> >>>
> >>> https://github.com/vjcitn/biocMultiAssay/blob/master/R/exs2se.R
> >>>
> >>>
> https://github.com/vjcitn/biocMultiAssay/blob/master/man/coerce-methods.Rd
> >>>
> >>> it occurred to me that we might want something like this in
> GenomicRanges
> >>> (that's where SummarizedExperiment is managed, right?) and I will add
> it
> >>> if there are no objections
> >>>
> >>> the arguments are currently
> >>>
> >>>       assayname = "exprs",    # for naming SimpleList element
> >>>       fngetter =
> >>>             function(z) rownames(exprs(z)),   # extract usable feature
> names
> >>>       annDbGetter =
> >>>            function(z) {
> >>>                clnanno = sub(".db", "", annotation(z))
> >>>                stopifnot(require(paste0(annotation(z), ".db"),
> >>> character.only=TRUE) )
> >>>                get(paste0(annotation(z), ".db"))  # obtain resource for
> >>> mapping feature names to coordinates
> >>>                },
> >>>       probekeytype = "PROBEID",   # chipDb field to use
> >>>       duphandler = function(z) {    # action to take to process
> duplicated
> >>> features
> >>>            if (any(isd <- duplicated(z[,"PROBEID"])))
> >>>                return(z[!isd,,drop=FALSE])
> >>>            z
> >>>            },
> >>>       signIsStrand = TRUE,   # verify that signs of addresses define
> strand
> >>>       ucsdChrnames = TRUE    # prefix 'chr' to chromosome token
> >>>
> >>>          [[alternative HTML version deleted]]
> >>>
> >>> _______________________________________________
> >>> Bioc-devel at r-project.org mailing list
> >>> https://stat.ethz.ch/mailman/listinfo/bioc-devel
> >>>
> >>
> >>         [[alternative HTML version deleted]]
> >>
> >> _______________________________________________
> >> Bioc-devel at r-project.org mailing list
> >> https://stat.ethz.ch/mailman/listinfo/bioc-devel
> >>
> >
> >
> > --
> > Computational Biology / Fred Hutchinson Cancer Research Center
> > 1100 Fairview Ave. N.
> > PO Box 19024 Seattle, WA 98109
> >
> > Location: Arnold Building M1 B861
> > Phone: (206) 667-2793
> >
> >
> > _______________________________________________
> > Bioc-devel at r-project.org mailing list
> > https://stat.ethz.ch/mailman/listinfo/bioc-devel
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> 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