Dear list,

I am modelling repeated measures data (called B) from M days (called
A) in a datebase, using lme()

from the nlme package. I have the following call for two random effecs
each with an AR1 error structure:*
*
(1) between days errors  a_i are assumed follow AR(1) structure (that
is corr(a_i,a_j)=\rho^{|i-j|})

(2) within day error b_ij are nested in days. b_i=[b_i1,...,b_iK] is the
within day error vector for day i.
     b_i ~MVN(0,Ó), for all i. It also follows AR(1) structure. For
different i i', b_i, b_i' are independent.

I don't know how to specify the mix model with two random effects in R. I
refer to lmer() in lme4 package. It can deal with two random effects, but
without AR(1) structure.

Many thanks if somebody can show me the correct form of random
arguments in lme().

Best Regards,
Bob Huang

	[[alternative HTML version deleted]]


