R 多元方差分析

当因变量（结果变量）不止一个时，可用多元方差分析（MANOVA）对它们同时进行分析。
以MASS包中的UScereal数据集为例（Venables，Ripley（1999）），我们将研究美国谷物中的卡路
里、脂肪和糖含量是否会因为储存架位置的不同而发生变化；其中1代表底层货架，2代表中层货
架，3代表顶层货架。卡路里、脂肪和糖含量是因变量，货架是三水平（1、2、3）的自变量。

> library(MASS) 
> attach(UScereal) 
> shelf <- factor(shelf) 
> y <- cbind(calories, fat, sugars) 
> aggregate(y, by=list(shelf), FUN=mean) 
 Group.1 calories fat sugars 
1 1 119 0.662 6.3 
2 2 130 1.341 12.5 
3 3 180 1.945 10.9 
> cov(y) 
 calories fat sugars 
calories 3895.2 60.67 180.38 
fat 60.7 2.71 4.00 
sugars 180.4 4.00 34.05 
> fit <- manova(y ~ shelf) 
> summary(fit) 
 Df Pillai approx F num Df den Df Pr(>F) 
shelf 2 0.402 5.12 6 122 1e-04 *** 
Residuals 62 
--- 
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
> summary.aov(fit)
Respon