[BioC] Combining expressionSets from GEO

Francois Pepin fpepin at cs.mcgill.ca
Wed Jan 30 18:41:22 CET 2008

Hi everyone,

I'm getting an error message when trying to combine two parts of a GSE

> tmp2<-combine(tmp[[1]],tmp[[2]])
Error in alleq(levels(x[[nm]]), levels(y[[nm]])) && alleq(x
[sharedRows,  :
  invalid 'x' type in 'x && y'

Checking to make sure that I should be able to combine them (from the
eSet documentation):

#eSets must have identical numbers of 'featureNames'
> all(featureNames(tmp[[2]])==featureNames(tmp[[2]]))
[1] TRUE

#must have distinct 'sampleNames'
> any(sampleNames(tmp[[1]])%in%sampleNames(tmp[[2]]))

#and must have identical 'annotation'.
> annotation(tmp[[2]])==annotation(tmp[[2]])
[1] TRUE

> sessionInfo()
R version 2.6.0 (2007-10-03)


attached base packages:
[1] tools     stats     graphics  grDevices utils     datasets  methods
[8] base

other attached packages:
[1] GEOquery_2.2.0 RCurl_0.8-1    Biobase_1.16.0

loaded via a namespace (and not attached):
[1] rcompgen_0.1-15

Does anyone know why that is happening and if there would be any way
around it?


More information about the Bioconductor mailing list