Skip to content

Mixed effects model help

3 messages · crevell@rice.edu, Fernando Henrique Ferraz P. da Rosa, Martin Maechler

#
I am attempting to set up a mixed effects model, with one factor fixed (put 
in as.factor, with three factor levels) and one factor random(put in 
as.factor, with four different batches) as well as the interactions.  I have 
tried the lme function, but R continues to crash everytime I run it...
My code is:
model_lme(y~a,random=~b+a:b)
Is this way off base?
I just need to get F statistics to do hypothesis testing.

Thanks,
Chris
#
This has come up less than a month ago:
https://stat.ethz.ch/pipermail/r-help/2004-November/060038.html

crevell at rice.edu writes:
--
Fernando Henrique Ferraz P. da Rosa
http://www.ime.usp.br/~feferraz
#
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