Displaying 1 result from an estimated 1 matches for "factaid".
Did you mean:
aacraid
2003 Dec 02
1
lme: reproducing example
...51,46,57,49,50,55,51,47)
Y[37:45]<-c(56,48,51,58,50,48,58,46,52)
df<-data.frame(id, factA, factB, Y)
df.aov <- aov(Y ~ factA*factB + Error(factA:id),
data=df)
summary(df.aov)
Is there a way to get the same results with lme as
with aov with Error()? HOW???
One idea was the following:
df$factAid=factor(paste(as.character(df$factA),":",as.character(df$id),sep=""))
df.lme <-
lme(Y~factA*factB,df,random=~1|factAid,method="REML")
The degrees of freedom look right, but the F values
don't match aov.
Hope somebody can help! Thanks!!
Karl