View Single Post
Old 04-28-2016, 05:38 AM   #8
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Code:
results(ABpairwise, contrast=c(1, -1, -1, 1))
You want the additional effect due to the interaction of treatment and genotype B. In other words, you need to take B.treat - B.ctl - (A.treat - A.ctl) (represented by a vector of "1, -1, -1, 1").

Regarding the LRT, you can do that too, but you'd need to be a bit more comfortable with what that's actually testing.
dpryan is offline   Reply With Quote