Re: [R] crossed random fx nlme lme4

2005-07-17 Thread Emilio A. Laca
random = pdBlocked( list( pdIdent( ~ y), pdIdent( ~ observer - 1), pdIdent( ~ set - 1) ) gave me output only when the data was a groupedData. The results are different depending on whether I specify observer or set as grouping factor. I am missing something ... Simon, Prof Bates, thanks

Re: [R] crossed random fx nlme lme4

2005-07-15 Thread Douglas Bates
On 7/13/05, Simon Blomberg [EMAIL PROTECTED] wrote: At 09:35 AM 14/07/2005, Emilio A. Laca wrote: 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

[R] crossed random fx nlme lme4

2005-07-13 Thread 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

Re: [R] crossed random fx nlme lme4

2005-07-13 Thread Simon Blomberg
At 09:35 AM 14/07/2005, Emilio A. Laca wrote: 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. Does