Как выполнить тест MANOVA в R?

Question

Я пытаюсь найти любое взаимодействие между независимыми переменными в Group и 2 образца зависимых переменных 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(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