[Bioc-devel] coerce ExpressionSet to SummarizedExperiment
    Vincent Carey 
    stvjc at channing.harvard.edu
       
    Sat Sep 20 16:38:33 CEST 2014
    
    
  
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]]
    
    
More information about the Bioc-devel
mailing list