[R] Mixed effects model help

Martin Maechler maechler at stat.math.ethz.ch
Fri Dec 10 09:13:57 CET 2004


>>>>> "Chris" ==   <crevell at rice.edu>
>>>>>     on Thu, 9 Dec 2004 22:59:13 -0000 writes:

    Chris> I am attempting to set up a mixed effects model, with
    Chris> one factor fixed (put in as.factor, with three factor
    Chris> levels) and one factor random(put in as.factor, with
    Chris> four different batches) as well as the interactions.
    Chris> I have tried the lme function, but R continues to
    Chris> crash everytime I run it...  My code is:

    Chris> model_lme(y~a,random=~b+a:b) 

    Chris> Is this way off base?  

well, yes: if this is your code and you don't get an error or a
strong warning, your version of R is so much outdated, that you
*really* should upgrade.

[The use of "_" as assignment operator is ``prehistoric''.
 It's still valid in some Jurassic implementations of the S language...]

    Chris> I just need to get F statistics to do hypothesis testing.

Martin




More information about the R-help mailing list