如何在R中进行MANOVA测试？

Question

我试图在Group中的自变量与Level1和Level2数据集中的两个因变量样本之间找到任何交互作用(缩小大小)。

Group Level1 Level2
a     1      0
a     2      3
a     4      3
b     2      4
b     1      3
b     3      2
c     2      4
c     3      2
c     1      3

为此，我想做一次马诺瓦测试。然而，当我试图运行manova(data[,2:3] ~ as.factor(Group), data=data)时，虽然data[,2:3]是一个列表，但会得到invalid type (list) for variable 'data[,2:3]'错误。

(显然，这是我第一次尝试使用manova和基本R来完成这项任务，而且我对R有非常基本的了解，所以在一般情况下，我可能完全错了。)

Answer 1

您需要使用cbind像这样运行它

数据

df <- read.table(header=T, text='Group Level1 Level2
a     1      0
a     2      3
a     4      3
b     2      4
b     1      3
b     3      2
c     2      4
c     3      2
c     1      3')

解决方案：

> manova( cbind(Level1,Level2) ~ Group, data=df)
Call:
   manova(cbind(Level1, Level2) ~ Group, data = df)

Terms:
                   Group Residuals
resp 1          0.222222  8.666667
resp 2                 2        10
Deg. of Freedom        2         6

Residual standard errors: 1.20185 1.290994
Estimated effects may be unbalanced