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.