[R] crossed random fx nlme lme4

[Emilio A. Laca]

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.