View Single Post
Old 07-13-2011, 10:23 PM   #2
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

Hi

you test for a rather odd contrast:
Quote:
Code:
fit1 <- fitNbinomGLMs(cdsFull, count ~ type + condition + type:condition)
fit0 <- fitNbinomGLMs(cdsFull, count ~ type + type:condition)
The standard ANOVA or ANODEV approach is to test first for main effects without any interaction terms present in the fit, and then for the interaction, with all main effect present in both fits, i.e.,

Code:
fit1 <- fitNbinomGLMs(cdsFull, count ~ type + condition + type:condition)
fit0 <- fitNbinomGLMs(cdsFull, count ~ type + condition)
For the pasilla data, this shouldn't give many hits, though. (If there were interactions, this would be rather worrying.)
Simon Anders is offline   Reply With Quote