a****m 发帖数: 693 | 1 【 以下文字转载自 Biology 讨论区 】
发信人: affymm (no), 信区: Biology
标 题: 怎样用R 来提取 micrarray 中的regressor variable 的值?
发信站: BBS 未名空间站 (Mon Nov 2 11:14:36 2009, 美东)
我们用Lmfit 来做回归,原理上是所有的点都有Signal = u + Line + Treat + e, 的
linear model, 我们只知道signal的值,其他的variable的值只能通过这个linear Model得
到, 问题是怎样在R中把细胞系(line)和处理效应(treat)的值取从Linear model得出来,只
知道 怎样取 residue, like fit=lmFit(mymean,design),
Res=residuals(fit,mymean),不知道这个
LmFit对mix model合适不合适? Signal=u+Line+Treat+Line*Treat+e (mixed model)?
谢谢 |
|