Question: Use interactions to get the difference in the effect of two treatments
gravatar for rrcutler
2 days ago by
rrcutler30 wrote:

Hello all,

I have two experimental conditions and one control, I want to compare the difference of these effects in reference to the control. A simple way to do this would to get the gene sets that result from comparing each experimental level to the control separately, and then take the difference of these sets. However, this might miss where the difference in effects is direction or magnitude. 

However, I am wondering if this can also be legitimately done using an interaction term? One would duplicate the control condition so that the setup would look like so (ignoring replicates here):

design(dds) <- ~Condition * Experimental_Condition
Sample  Condition Experimental_Condition
cond1 Experimental cond1
cntr Control cond1
cond2 Experimental cond2
cntr(duplicated) Control cond2

Note that the second control is just a duplicate of the first. 

The resulting genes from the interaction term, conditionExperimental.Experimental_Conditioncond2, would be where the effect was significantly different between the condition effects in reference to the control. I realize this is similar to the LRT test, where we would be looking to see if there are any differences at all between all the levels. But, we will still not get only the genes where the effect is different. Please correct me if I'm wrong. 

I also realize that this might mess with the dispersion estimates, as there is now a pseudo replicate.



ADD COMMENTlink written 2 days ago by rrcutler30
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: 409 users visited in the last hour