> 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版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请联系我们举报，一经查实，本站将立刻删除。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/211413.html原文链接：https://javaforall.net
