Giuseppe GIven the information below it is difficult to provide help. But, you might use the lmer() function in the Matrix package in a fashion such as
lmer(A ~ B +(1|ID), data) where ID is a grouping variable. There are plently of papers on lme and lmer in the R news as well as a book by Pinhiero and Bates as well as multiple threads in the archives on this function as well. -----Original Message----- From: [EMAIL PROTECTED] on behalf of [EMAIL PROTECTED] Sent: Wed 11/30/2005 5:19 AM To: [email protected] Cc: Subject: [R] Random Effects for One-Way Anova Hello to All. I'd want to use a one-way ANOVA. This means that I have only one factor, with, lets say, 5 levels. I made a dataframe, called "DATA", with two Columns: A, that is my response, and it is "class numerical". B, that defines the different levels of my factor, and it is "class factor". If I want to use a fixed effect model, I know that the formula I have to use is: lm.1 <- lm(A~B,data=DATA) anova(lm.1) My questions is: which formula should I use if I want to use a random effects model? I think I should use "lme", but I don't know how. I hope to get a reply Best wishes to all Giuseppe ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html [[alternative HTML version deleted]] ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
