This is a test version of Biostars. For the public version, visit https://www.biostars.org.
EdgeR analysis with 3 conditions

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
• 4,387 views • link
updated by Ram Ram • written by Matina Matina

1 answer

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.

• 0 views • link
written by Devon Ryan Devon Ryan

Thank you so much!

• 0 views • link
updated by Ram Ram • written by Matina Matina

Log in to answer this question.

More posts like this