Question: Comparisons of comparisons in DESeq2
1
gravatar for Iora
12 months ago by
Iora10
Iora10 wrote:

Hi all,

I have an experiment where two treatment conditions TR1 and TR2 and their respective controls CL1 and CL2.

With the pairwise comparisons, I can find the differences within the respective Treatment conditions compared to their respective controls. I am also doing the pairwise comparisons to find the differences between the controls too. 

However, I would also like to compare the differences between the respective treatments after it is already compared with their respective controls. What I mean here is that I would like the make a final comparison of.(TR1 vs. CL1) vs. (TR2 vs. CL2)

Do the following design and contrast captures this treatment:

trtcontrast <- list(c('TR1','CL1'), c('TR2','CL2') )

res <- results( dds, contrast=trtcontrast , listValues=c(1/2,-1/2), alpha=0.05 )

As I understood, I do not think this is the right way as is it comparing the average of samples (TR1, CL1) and (TR2, CL2).

What would be the right way to do that specific comparison?

Thanks in advance.

deseq2 • 156 views
ADD COMMENTlink modified 12 months ago by Michael Love24k • written 12 months ago by Iora10
Answer: Comparisons of comparisons in DESeq2
2
gravatar for Michael Love
12 months ago by
Michael Love24k
United States
Michael Love24k wrote:

If you name the variables condition with levels (ctrl, trt) and group with levels (A, B), then you can use a design of ~group + group:condition. This will produce a group A trt vs ctrl effect and a group B trt vs ctrl effect. You pull these out with results(dds, name="...").

Then to contrast the two effects you can use:

results(dds, contrast=list("...", "..."))

where you fill in the ... with the two treatment effects for the two groups.

ADD COMMENTlink written 12 months ago by Michael Love24k

Thanks a lot, Michael. I get the idea now and understood what you mean. 

ADD REPLYlink written 12 months ago by Iora10
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 16.09
Traffic: 317 users visited in the last hour