[R] AOV() may misslabel random effects.

peter dalgaard pdalgd at gmail.com
Mon Mar 14 20:58:17 CET 2011


On Mar 14, 2011, at 17:57 , Giles Crane wrote:

> Greetings,
> 
> The aov() function may mislabel
> the random effects as in the example below:
> Has anybody else noticed this?

What's "mislabeled" about it??? Looks like you nave an unbalanced design (in which case, aov() may be the wrong tool.)

-pd

> 
> Cordially,
> Giles Crane, MPH, ASA, NJPHA
> gilescrane at verizon.net
> 
> > m2
> 
> Call:
> aov(formula = y ~ ap + pe + Error(ju), data = d)
> 
> Grand Mean: 77.50667
> 
> Stratum 1: ju
> 
> Terms:
>                      ap
> Sum of Squares  4322.538
> Deg. of Freedom       12
> 
> 13 out of 25 effects not estimable
> Estimated effects may be unbalanced
> 
> Stratum 2: Within
> 
> Terms:
>                      ap       pe Residuals
> Sum of Squares  7047.885  255.034  2981.290
> Deg. of Freedom       25        2        35
> 
> Residual standard error: 9.229285
> Estimated effects may be unbalanced
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

-- 
Peter Dalgaard
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: pd.mes at cbs.dk  Priv: PDalgd at gmail.com



More information about the R-help mailing list