[R] ANOVA procedure on the sufficient statistics

Michael Friendly friendly at yorku.ca
Thu Feb 19 14:45:01 CET 2004

>You can apply the trick that Prof. Dalgaard recently posted in response to a
>similar question (for one-way ANOVA).  For each cell, generate data as:
>y <- cell.mean + cell.sd * scale(rnorm(cell.count))
>Then generate the data frame to feed to aov.
That will generate only approximately the same means and MSE, however.

Larsen's procedure generates a weighted data set that gives the same 
ANOVA table
as the raw data, but requires that the weight= argument be used in aov().
(weights are not supported in all related functions)


Michael Friendly     Email: friendly at yorku.ca 
Professor, Psychology Dept.
York University      Voice: 416 736-5115 x66249 Fax: 416 736-5814
4700 Keele Street    http://www.math.yorku.ca/SCS/friendly.html
Toronto, ONT  M3J 1P3 CANADA

More information about the R-help mailing list