24小时热门版块排行榜    

CyRhmU.jpeg
查看: 340  |  回复: 0

xo5547

铁虫 (小有名气)

[求助] 已经得到基因差异的数据,是.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”是怎么设计的?能具体说下吗?非常感谢,急需答复,谢谢高手指导!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xo5547 的主题更新
信息提示
请填处理意见