| 查看: 420 | 回复: 0 | ||
[求助]
已经得到基因差异的数据，是.txt格式的，不知用limma如何分析。
|
|
有关limma进行数据分析的问题。我是R的新手，现在正头疼怎样用limma进行基因矩阵数据的分析，已经得到基因差异的数据，就是不知道如何分析。 以下是站内的一组代码： library(limma) eset<-read.table("eset.txt",sep="\t",header=T,row.names=1)#基因表达谱 design<-read.table("design.txt",header=T,sep="\t"  ##对样本两类样本（正常和疾病）的设计fit<-lmFit(eset,design) fit<-eBayes(fit) DERes<-topTable(fit,coef=colnames(design),number=length(fit)) write.table(DERes,"DEG.txt",sep="\t",quote=F) 想问下“design.txt”是怎么设计的？能具体说下吗？非常感谢，急需答复，谢谢高手指导！ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
