由买买提看人间百态

topics

全部话题 - 话题: micrarray
(共0页)
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)?
谢谢
b******s
发帖数: 1089
2
做生物的如果按部就班得做,其实就那么几板斧。
你做的是不是Brachypodium?还是水稻?都已经有表型了,应该很好做了。
可以从表型分析开始,根据你对project的期望来做。决定生物表型的已经因素并不太
多。最容易着手的是plant hormone。有表型的,hormone多少会有反应。
如果是直接的,恭喜你,你有cns的潜力。间接的,也可能发个小文章。
具体做法无非是首先从生理上确定有没有hormone involvement,然后genetics杂交,
然后做到细胞学的机理。
从表型开始做,需要很多genetics和physiology的经验。
要么从gene本身做起。你已知蛋白是ER膜定位,那就看其在secretion或者protein
trafficking中有没有作用。整个secretion system有很多marker可供细胞学研究。也
可以做parmacological的实验,看看有没有哪个pathway被影响了。下面的手段根据实
验室条件,要是细胞学强,可以多做细胞学,要是生化强,可以做blot。从gene本身开
始还可以,如你说的,Y2H screen这个蛋白... 阅读全帖
s*r
发帖数: 2757
3
来自主题: Statistics版 - 用R画boxplot的问题
时间一定要是factor吗才能画四个并排的boxplot?
如果是的话,一个笨办法是把时间变成char/string,用040替代40,用080替代80
然后再转化为factor用于作图
是micrarray data吗, 我记得rma有内置的画boxplot的方法
(共0页)