[R] aov with Error and lme

Karl Knoblick karlknoblich at yahoo.de
Tue Nov 18 15:00:28 CET 2003


Thanks!

But I could not find it in MASS 2nd edition, also not
a similiar example to mine in Pinheiro Bates. I
thought at least for a balanced and quite simple
design the results should be equal. Is that wrong?

Karl.

 --- Prof Brian Ripley <ripley at stats.ox.ac.uk>
schrieb: > There are worked examples in MASS (at least
the last
> two editions) of
> getting `the same results with lme as with aov with
> Error()'.
> 
> Your over-long lines make it difficult to see, but
> this appears to be very 
> similar to examples both in MASS and in Pinheiro &
> Bates.  Have you 
> consulted either?
> 
> 
> On Tue, 18 Nov 2003, Karl Knoblick wrote:
> 
> > 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)
> 
> -- 
> Brian D. Ripley,                 
> ripley at stats.ox.ac.uk
> Professor of Applied Statistics, 
> http://www.stats.ox.ac.uk/~ripley/
> University of Oxford,             Tel:  +44 1865
> 272861 (self)
> 1 South Parks Road,                     +44 1865
> 272866 (PA)
> Oxford OX1 3TG, UK                Fax:  +44 1865
> 272595
>  


__________________________________________________________________

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




More information about the R-help mailing list