EdgeR Warning message In estimateCommonDisp.
0
0
Entering edit mode
22 months ago
Dr.Animo ▴ 10

Hi everyone, I am trying to perform DEG analysis. I've two samples each has 3 replicates (one is XR and other is YR) replicates are XR1,XR2,XR3 and YR1,YR2,YR3..

d<- DGEList(counts=as.matrix(raw.data[,2:7]), lib.size=library.sizes, group=expt)
d<- calcNormFactors(d)

At the following step, I got a warning message

d<- estimateCommonDisp(d)
Warning message:
In estimateCommonDisp.DGEList(d) :
  There is no replication, setting dispersion to NA.
d$common.dispersion
[1] NA

How to resolve this?

R rna-seq edgeR DEG • 809 views
ADD COMMENT
1
Entering edit mode

what's in ur expt variable ?? how many replicates ??

ADD REPLY
1
Entering edit mode

(for the replicates: read the question)

ADD REPLY
1
Entering edit mode

Does expt contain the "XR" / "YR" status, or does it encode the "XR1", "XR2", ... status. If the latter, (as I suspect) then it will be assumed that there are six separate groups, each with a single sample.

ADD REPLY
1
Entering edit mode

Thanks for the response. I've fixed it the problem is due to expt

ADD REPLY
0
Entering edit mode

If you have these groups like XR1, XR2, XR3 you would need to indicate them like:

Name    Group
XR1 XR
XR1 XR
XR3 XR
YR1 YR
YR2 YR
YR3 YR

and then make the design as ~Group.

ADD REPLY

Login before adding your answer.

Traffic: 2811 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6