I need to specify a model similar to this

lme.formula(fixed = sqrt(lbPerAc) ~ y + season + y:season, data = cy,
     random = ~y | observer/set, correlation = corARMA(q = 6))

except that observer and set are actually crossed instead of nested.

observer and set are factors
y and lbPerAc are numeric

If you know how to do it or have suggestions for reading I will be  
grateful.


eal

ps I have already read Pinheiro & Bates, the jan 05 newsletter, and  
several postings.

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Reply via email to