R中的共色散
我尝试使用edge R软件包来分析我的数据,如下所示R中的共色散,r,R,我尝试使用edge R软件包来分析我的数据,如下所示 Warning message: In estimateGLMCommonDisp.default(y = y$counts, design = design, : No residual df: setting dispersion to NA 这是我的数据帧,叫做子数据 A B C
Warning message:
In estimateGLMCommonDisp.default(y = y$counts, design = design, :
No residual df: setting dispersion to NA
这是我的数据帧,叫做子数据
A B C D E
1 13707 13866 12193 12671 10178
2 0 0 0 0 1
3 7165 5002 1256 1341 2087
6 8537 16679 9042 9620 19168
10 19438 25234 15563 16419 16582
16 3 3 11 3 5
genotype=factor(c("LMS1","LRS1","MS3","MS4","RS5"))
y=DGEList(counts=data.matrix(subdata),group=genotype)
design=model.matrix(~0+group,data=subdata)
y=estimateGLMCommonDisp(y,design, verbose=TRUE)
我试图计算共同的离散度来估计上调和下调基因,但我得到如下错误信息
Warning message:
In estimateGLMCommonDisp.default(y = y$counts, design = design, :
No residual df: setting dispersion to NA
请任何人帮我解决这个问题,我真的很感激所以你对每个小组都有一个单独的观察结果?是的,我对每个小组都有一个单独的观察结果。请阅读张贴指南。特别是,请发布返回警告消息的确切命令。我还强烈建议您使用
dput
发布您的数据,使我们能够轻松重现您的问题。最后,警告不是错误消息。您的计算已成功完成,但如前所述,由于缺少残差,因此无法计算离散值。在示例中,计数和数据参数不一样,这是否有帮助<我非常感谢你的帮助