Question: Should I contrast 2 treatment groups using the control group as reference or directly against each other?
gravatar for serpalma.v
5 days ago by
serpalma.v20 wrote:

Dear community,

I am analyzing RNAseq data set with DESeq2. Samples are grouped in three groups: one control group (n=16) and two case groups (case 1 and 2, n=6 and n=8, respectively).

I want to find differentially expressed genes of both case groups against the control group, but also against each other (case1 vs case2).

When fitting a GLM, I can easily get the coefficients for case1_vs_control and case2_vs_control, and to compare case1 against case2, I could do case1_vs_control *vs* case2_vs_control, and that would produce differentially expressed genes whose log2FC over the control group differ between case1 and case2.

However, because of the degrees of freedom, I know a different result will come out if I directly compare case1 and case2, ignoring the control group. This also would require fitting another GLM for a two group comparison.

My question is which approach is more robust to compare case1 and case2.


Thanks in advance!

ADD COMMENTlink modified 2 days ago by Michael Love15k • written 5 days ago by serpalma.v20
gravatar for Michael Love
2 days ago by
Michael Love15k
United States
Michael Love15k wrote:

The recommended way in DESeq2 would be to run DESeq() with a design ~condition, where condition has levels "control", "case1", and "case2", and then to use the results() function to build the tables:

res1 <- results(dds, contrast=c("condition","case1","control"))
res2 <- results(dds, contrast=c("condition","case2","control"))
res3 <- results(dds, contrast=c("condition","case2","case1"))
ADD COMMENTlink modified 2 days ago • written 2 days ago by Michael Love15k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 410 users visited in the last hour