I have problem with pairwise comparison in estimate dispersion functions in edgeR pipeline
1
0
Entering edit mode
2.3 years ago
seda ▴ 10

Hi everyone!

I have silly questions, I checked edgeR user's guide and I didn't find any answer to my question.

I have two different samples. I want to compare with each other and find differential expressions. So, I run this code : y1 <- estimateCommonDisp(y, verbose = TRUE) and it gave me this message :

Warning message: In estimateCommonDisp.DGEList(y, verbose = TRUE) : There is no replication, setting dispersion to NA.

I get the same warning when I use estimateDisp() and estimateTagwiseDisp() functions.

What does it mean 'there is no replication' because I already have two samples, how can I create the replication?
Do you have any suggestion to find the differential expression for two samples?

Thank you all!

edgeR dispersion comparison estimateDisp estimate • 1.0k views
ADD COMMENT
1
Entering edit mode

Hi, If you have no replicate for your samples, then you can also use NOISeq package. The Pairwise Differential Expression Analysis (Without Replicates) functionality offers a strategy for analyzing RNA-seq datasets that do not have replicates. It is based on the software package NOISeq, which belongs to the Bioconductor project. NOISeq is a novel nonparametric approach for the identification for differentially expressed genes from RNA-Seq count data.

ADD REPLY
3
Entering edit mode
2.3 years ago
Gordon Smyth ★ 7.2k

See Section 2.12 "What to do if you have no replicates" in the edgeR User's Guide.

Replication means having at least two samples in at least one of your groups. Having two samples, one in each group, is not replication. Replication is a fundamental statistical concept that allows statistics programs like edgeR to estimate the amount of biological variablity that is present.

ADD COMMENT
0
Entering edit mode

I think you meant Replication means having at least two samples in each one of your groups.

ADD REPLY
1
Entering edit mode

No, I meant "in at least one group". It is not necessary that group sizes be balanced or that every group have replicates.

ADD REPLY

Login before adding your answer.

Traffic: 2123 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