我有一个ExpressionSet
对象,我想要分组。例如,
> str(ESet)
Formal class 'ExpressionSet' [package "Biobase"] ..
..@ assayData :..
..@ phenoData :
.. .. .. ..$ STATUS : num [1:210] 1 1 1 1 1 1 1 1 1 1 ...
....
我想提取STATUS==0
的子集。我试过了:
exprs(ESet@phenoData$STATUS==0)
但它不起作用。
答案 0 :(得分:2)
你快到了。猜测你的数据结构,我认为以下内容应该有效:
exprs(ESet)[ESet@phenoData$STATUS==0,]
如果你看一下这个paper,就会有很多根据应变和时间对数据进行子设置的例子。我应该承认我是本文的作者。