r - 如何在 R 中执行 MANOVA 测试？

Question

我试图在我的数据集中找到自变量和Group2 个因变量样本之间的任何交互（减小了大小）。Level1Level2data

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 测试。但是，当我尝试运行时，虽然是一个列表manova(data[,2:3] ~ as.factor(Group), data=data)，但我得到了invalid type (list) for variable 'data[,2:3]'错误。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