[R] Controlling the number of interactions of a lme

[Ben Bolker]

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.