我试图在许多因变量上运行anova。我有一个自变量,这是我的分组变量(组)。我有大约25个DV - “TMTG,TMTF,CUE,CSE,TCUE,TCSE,WRS,WMAO,TWRS,TWMAO,JCP,JCPE ....等”。我使用以下代码作为前三个变量,我得到了所需的输出。如何调整代码以同时获取所有25个变量的输出,但是没有命名它们?我有另外一个100 DV的数据集 - 我不能写出来!
这是数据框
Group TMTG TMTF CUE CSE WRS
TN 27 33 35.12 13.56 0
TN 32 34 12.90 25.56 0
TN 14 78 11 14.78 0
TN 89 41 98 45.25 0
TL 65 11 18.5 23.89 0
TL 12 78 34.6 41.85 0
TL 11 20 35.5 45.5 0
TL 27 25 11.28 55.69 0
以下是代码:
mydataframe
manova_1 <-
manova(cbind(TMTG, TMTF, CUE) ~ as.factor(Group), data = mydataframe)
manova_1
summary.aov(manova_1)
这是输出
Response TMTG :
Df Sum Sq Mean Sq F value Pr(>F)
as.factor(Group) 1 0.535 0.5351 0.1683 0.6858
Residuals 21 66.769 3.1795
Response TMTF :
Df Sum Sq Mean Sq F value Pr(>F)
as.factor(Group) 1 0.02 0.016 5e-04 0.9831
Residuals 21 749.13 35.673
Response CUE :
Df Sum Sq Mean Sq F value Pr(>F)
as.factor(Group) 1 14.7 14.75 0.0372 0.8489
Residuals 21 8325.7 396.46
我想调整这一行:
manova(cbind(TMTG,TMTF,CUE)~as.factor(Group),data = mydataframe,
所以cbind可以接收所有列,而不必将其写出来。我试过cbind(2:24),但它不起作用!任何帮助将不胜感激!!!
答案 0 :(得分:1)
假设1)Group
是mydataframe
中的第一个变量,而2)你想做一个manova而不是一些单独的anovas,你可以替换这行:
manova(cbind(TMTG, TMTF, CUE) ~ as.factor(Group), data = mydataframe)
使用:
manova(as.matrix(mydataframe[, -1]) ~ as.factor(Group), data = mydataframe)