最近想用roast()函数,看了文档,没看明白
roast(y, index=NULL, design=NULL, contrast=ncol(design), set.statistic="mean",
gene.weights=NULL, array.weights=NULL, weights=NULL, block=NULL, correlation,
var.prior=NULL, df.prior=NULL, trend.var=FALSE, nrot=999)
这个里面的design参数解释是design matrix,我不明白什么是design matrix. 看例子:design <- cbind(Intercept=1,Group=c(0,0,1,1))
这个里面Intercept和Group是设置什么的
还有那个index直接是设置行号吗?
我的数据值20*40的矩阵,前20列为一类样本,剩下的为2类样本
哪位帮帮忙.多谢啦 |