[R] aov with Error and lme

[Karl Knoblick]

Hi

I searched in the list and only found questions
without answers e.g.
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/19955.html
: Is there a way to get the same results with lme as
with aov with Error()?

Can anybody reproduce the following results with lme:
id<-c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5)
factA<-c("al","al","al","al","al","al","al","al","al","al","al","al","al","al","al","a2","a2","a2","a2","a2","a2","a2","a2","a2","a2","a2","a2","a2","a2","a2","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3")
factB<-c("B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3")
Y<-c(56,52,48,57,54,46,55,51,51,58,51,50,54,53,46,54,50,49,53,49,48,56,48,52,52,52,50,55,51,46,57,49,50,55,51,47,56,48,51,58,50,48,58,46,52)

df<-data.frame(id, factA, factB, Y)

df$id<-as.factor(df$id)
df.aov <- aov(Y ~ factA*factB + Error(factA:id),
data=df)
summary(df.aov)

Hope, somebody can help!

Karl





__________________________________________________________________

Gesendet von Yahoo! Mail - http://mail.yahoo.de
Logos und Klingeltöne fürs Handy bei http://sms.yahoo.de