如何在R中进行MANOVA测试?

时间:2015-02-27 15:27:35

标签: r statistics manova

我试图找到Group中的独立变量与我的Level1数据集中的两个因变量Level2data样本之间的任何互动(缩小尺寸)。

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]错误。

(显然这是我第一次尝试使用manovaR来完成这项任务,而且我对R有非常基本的了解,所以也许我完全错了一般来说)。

1 个答案:

答案 0 :(得分: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