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
Mixed effects model help
3 messages · crevell@rice.edu, Fernando Henrique Ferraz P. da Rosa, Martin Maechler
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:
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
______________________________________________ R-help at stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
-- Fernando Henrique Ferraz P. da Rosa http://www.ime.usp.br/~feferraz
"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