Re: [R] random effects with lme() -- comparison with lm()

2004-01-16 Thread Douglas Bates
Jerome Asselin [EMAIL PROTECTED] writes: On Thu, 2004-01-15 at 16:30, Douglas Bates wrote: ...snip... (BTW, I wouldn't say that this is equivalent to a fixed effects model. It is still a random effects model with two variance components. It just doesn't have well-defined estimates for

[R] random effects with lme() -- comparison with lm()

2004-01-15 Thread Jerome Asselin
Hi all, In the (very simple) example below, I have defined a random effect for the residuals in lme(). So the model is equivalent to a FIXED effect model. Could someone explain to me why lme() still gives two standard deviation estimates? I would expect lme() to return either: a) an error or a

Re: [R] random effects with lme() -- comparison with lm()

2004-01-15 Thread Douglas Bates
Jerome Asselin [EMAIL PROTECTED] writes: Hi all, In the (very simple) example below, I have defined a random effect for the residuals in lme(). So the model is equivalent to a FIXED effect model. Could someone explain to me why lme() still gives two standard deviation estimates? I would

Re: [R] random effects with lme() -- comparison with lm()

2004-01-15 Thread Jerome Asselin
On Thu, 2004-01-15 at 16:30, Douglas Bates wrote: ...snip... (BTW, I wouldn't say that this is equivalent to a fixed effects model. It is still a random effects model with two variance components. It just doesn't have well-defined estimates for those two variance components.) Agreed.