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
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
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
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