[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