Search
Question: Should I contrast 2 treatment groups using the control group as reference or directly against each other?
2
gravatar for serpalma.v
5 months ago by
serpalma.v30
Germany
serpalma.v30 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 5 months ago by Michael Love18k • written 5 months ago by serpalma.v30
2
gravatar for Michael Love
5 months ago by
Michael Love18k
United States
Michael Love18k 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 5 months ago • written 5 months ago by Michael Love18k
Please log in to add an answer.

Help
Access

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