我试图找到Group
中的独立变量与我的Level1
数据集中的两个因变量Level2
和data
样本之间的任何互动(缩小尺寸)。
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
有非常基本的了解,所以也许我完全错了一般来说)。
答案 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