Hi,
I'm trying to understand why FeatureSet objects behave slightly different
than eSet objects.
Here's the one example I'm trying to work out:
if (!require(pd.hugene.1.0.st.v1)){
library(BiocInstaller)
biocLite('pd.hugene.1.0.st.v1')
}
library(oligoData)
data(affyGeneFS)
affyGeneFS
data(sample.ExpressionSet)
sample.ExpressionSet
## subset ExpressionSet
## everything ok
sample.ExpressionSet[c(1, 1),]
## subset FeatureSet
## error: featureNames differ between assayData and featureData
affyGeneFS[c(1, 1),]
But FeatureSets are derived from NChannelSet objects... so:
example('NChannelSet-class')
obj
obj[c(1, 2),] ## OK
obj[c(1, 1),] ## not OK
I was wondering why/if this is intended (i.e., it works on "single channel"
eSets, but fails on NChannelSets)?
Thank you so much for any insight,
benilton
[[alternative HTML version deleted]]
_______________________________________________
Bioc-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel
