anova - 同时选择多个DV

时间:2017-11-22 21:00:55

标签: r anova

我试图在许多因变量上运行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),但它不起作用!任何帮助将不胜感激!!!

1 个答案:

答案 0 :(得分:1)

假设1)Groupmydataframe中的第一个变量,而2)你想做一个manova而不是一些单独的anovas,你可以替换这行:

manova(cbind(TMTG, TMTF, CUE) ~ as.factor(Group), data = mydataframe)

使用:

manova(as.matrix(mydataframe[, -1]) ~ as.factor(Group), data = mydataframe)