edgeR差异表达分析报错求解

用edgeR做有重复的差异表达分析时出现如下报错:Error in estimateDisp.default(y = y$counts, design = design, group = group,  : 
  object 'prior.n' not found

求解

所用脚本如下:

rawdata<-read.csv("deseq2.csv",fill=TRUE,header = T)
rownames(rawdata)<-rawdata$gene_id
rawdata<-rawdata[,-1]
head(rawdata)
condition <- factor(c(rep("BWRF",3),rep("HRF",3),rep("JF",3),rep("RF",3),rep("WRF",3)),levels=c("BWRF","HRF","JF","RF","WRF"))
condition 
#构建DGEList矩阵
colData <- data.frame(row.names=colnames(rawdata),condition)
colData
y <- DGEList(counts = rawdata, group =colData$condition)
y
#过滤low counts数据
keep <- rowSums(cpm(y)>1)>=24
y <- y[keep,keep.lib.sizes=FALSE]
#标准化(以TMM)
y <- calcNormFactors(y,method = "TMM")
#差异表达分析
design <- model.matrix(~0+ condition,data =colData)
colnames(design) <- levels(colData$condition)
head(design)
design
y<- estimateDisp(y,design)
 

求解决方法,请各位指教!

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值