[R-sig-ME] Help understanding residual variance
Ben Bolker
bbolker at gmail.com
Thu Mar 29 20:25:27 CEST 2012
Joshua Wiley <jwiley.psych at ...> writes:
>
> Hi Ista,
>
> To me the onis is on the statistician consultant to explain *why* you
> cannot have both random intercepts and slopes. Does the consultant
> have papers to reference or proofs?
>
> In any case, this is hardly exclusive to 'R doing something strange'.
> SAS and Stata happily join the gang. See the attached file for code
> and output from all three using a minidataset simulated in R.
>
> I suppose one could bicker over whether a random intercept and slope
> is a good idea, but possible it certainly is. You might suggest that
> it is poor fare to voice strong opinions about matters which one does
> not understand.
>
> Cheers,
>
> Josh
Very nice example. Do note that while R::lme and Stata give
the same point estimates, Stata also provides estimated confidence
intervals, which are enormous -- suggesting that, while one can
do this, the resulting model might be unidentifiable.
Doug Bates commented off-list that:
> I believe that the estimates for such a model are at least ill-defined
> if not unidentifiable.
cheers
Ben
More information about the R-sig-mixed-models
mailing list