Controlling the number of interactions of a lme
cleberchaves <cleberchaves <at> gmail.com> writes:
[snip]
My model have many response variables and when i run the anova, the number of interactions (up to six) is great and the p-values of all variables not appear. I wanted to know if i could to control the number of interactions of the
[snip]
v.is<-lme(is~direction*envir*region*hour*estom*esl, random=~1|ind/dir/reg, tabela) anova(v.is,test="F")
You probably want something like is ~ (direction+envir+region+hour+estom+esl)^2 for example, which would include the main effects and all two-way interactions. See the "Details" section of ?formula for a (terse) description.