我有两组人,GRP0 和 GRP1,我测量了三个连续变量:VAR1、VAR2 和 VAR3。
我想在 R 中使用 Mancova: - VAR1、VAR2 和 VAR3 作为结果变量 - GRP={0,1} 作为预测变量 - 年龄和性别作为协变量
在 R 中制定这个模型的正确方法是什么?
此外,测量是在 100 个实例(它们是串行相关的)上进行的,因此任何关于如何在顶部应用基于排列的多重比较校正的帮助都是理想的。
非常感谢
我有两组人,GRP0 和 GRP1,我测量了三个连续变量:VAR1、VAR2 和 VAR3。
我想在 R 中使用 Mancova: - VAR1、VAR2 和 VAR3 作为结果变量 - GRP={0,1} 作为预测变量 - 年龄和性别作为协变量
在 R 中制定这个模型的正确方法是什么?
此外,测量是在 100 个实例(它们是串行相关的)上进行的,因此任何关于如何在顶部应用基于排列的多重比较校正的帮助都是理想的。
非常感谢
见: http: //www.statmethods.net/stats/anova.html
Y <- cbind(VAR1, VAR2, VAR3)
fit <- manova(Y ~ GRP+age+gender)
summary(fit, test="Pillai")
其他测试选项是“Wilks”、“Hotelling-Lawley”和“Roy”。
summary.aov(fit) # for univariate statistics.
以下链接可能有助于基于排列的多重比较:http ://www.ncbi.nlm.nih.gov/pmc/articles/PMC2611984/pdf/1751-0473-3-15.pdf
此页面上的 R 代码也很有帮助:http ://biostatistics1014.blogspot.in/2013/04/one-way-anova-permutation-test-and.html