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:	r-help-bounces@stat.math.ethz.ch on behalf of Giuseppe.Palermo@bo.infn.it
Sent:	Wed 11/30/2005 5:19 AM
To:	r-help@stat.math.ethz.ch
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

______________________________________________
R-help@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




	[[alternative HTML version deleted]]


