Как выполнить тест 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