EdgeR analysis with 3 conditions
Entering edit mode
9.4 years ago
Matina ▴ 230

Hi guys,

I have human samples coming from 3 conditions (control, breast cancer and endometrial cancer). I want to do pairwise DE analysis between conditions using EdgeR but I'm not sure how to set up the samples table. Should I use a table with counts from all conditions and compare the conditions of my interest specifying the pair in the exact test function or set up a different table each time with the counts from the conditions i want to compare? Ive tried both ways and they seem to give me different lists of DE genes so I was wondering which one is the correct way to do it. My guess is that the differences come from the estimation of dispersion for each gene.

Sorry if my question is a bit basic but im a beginner.

Thanks a lot for any comments or input you might have.

RNA-Seq edgeR • 3.5k views
Entering edit mode
9.4 years ago

If you're going to make pairwise comparisons anyway, then keep everything in one table with a single design matrix and then just use a contrast with the last comparison. There's an example of exactly this in section 2.9.3 of the edgeR manual. You're correct that the difference between doing things this way and making separate pairwise tests with different count tables is in the dispersion estimation, which should be more reliable with the GLM.

Entering edit mode

Thank you so much!


Login before adding your answer.

Traffic: 1925 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6