Question

我想将为aov（）编写的代码转换为car-package中的Anova（）函数。

anovadata3 <- within(anovadata3, {
  subject <- factor(subject)
  time <- factor(time)
  gender <- factor(gender)
  group <- factor(group)
  groupgender <- factor(groupgender)
})

anovadata3.aov <- aov(values ~ time*group*gender + Error(subject),
                      data = anovadata3)
summary(anovadata3.aov)

此代码为我提供以下输出：

Error: subject
              Df Sum Sq Mean Sq F value  Pr(>F)   
group          1  32220   32220   8.632 0.00365 **
gender         1     30      30   0.008 0.92819   
group:gender   1     15      15   0.004 0.94952   
Residuals    221 824913    3733                   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Error: Within
                   Df Sum Sq Mean Sq F value   Pr(>F)    
time                3  21160    7053   9.223 5.53e-06 ***
time:group          3  18338    6113   7.993 3.06e-05 ***
time:gender         3   1916     639   0.835  0.47486    
time:group:gender   3  11679    3893   5.091  0.00172 ** 
Residuals         663 507012     765                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

然而，当我尝试使用car-package中的Anova（）函数时，我会这样做：

library(car)
mlm <- lm(values ~ time*group*gender, data = anovadata3)
rfactor <- factor(c("time", "gender","group","groupgender","subject"))
anovadata3.aov <- Anova(mlm, idata = dataframe(rfactor), idesign = ~ rfactor, type ="III")
summary(anovadata3.aov)

这给了我这个输出。

     Sum Sq                Df         F value              Pr(>F)       
 Min.   :     58.3   Min.   :  1   Min.   :  0.03871   Min.   :0.00000  
 1st Qu.:   1324.9   1st Qu.:  1   1st Qu.:  0.54230   1st Qu.:0.04997  
 Median :  10281.5   Median :  3   Median :  1.57697   Median :0.21357  
 Mean   : 196286.3   Mean   :100   Mean   : 34.03103   Mean   :0.31053  
 3rd Qu.:  12290.4   3rd Qu.:  3   3rd Qu.:  2.61758   3rd Qu.:0.50989  
 Max.   :1331924.5   Max.   :884   Max.   :262.67095   Max.   :0.84408  
                                   NA's   :1           NA's   :1

有谁知道我如何重新制作我用于aov（）以适合Anova（）的代码。我试着按照以下教程： https://gribblelab.wordpress.com/2009/03/09/repeated-measures-anova-using-r/

尝试让Anova（）正确无误。但它没有给出看起来相似的输出。我也从网页上看到它假设给Mauchlys和Greenhouse，我没有得到。也有人知道如何在anova结果中获得eta平方吗？或者是否有必要使用单独的函数来计算eta（etaSquared（））。

以下数据用于测试，我试图测试时间，性别和群体之间的“值”与因素之间的交互效应是否存在显着差异。

          values  testperiod subject gender group groupgender time
1   118.82660110     Pretest       1      2     2   BSTfemale    1
2    61.07615138     Pretest       2      2     2   BSTfemale    1
3    57.51022740     Pretest       3      2     2   BSTfemale    1
4    70.73637347     Pretest       4      2     2   BSTfemale    1
5     9.86907880     Pretest       5      2     2   BSTfemale    1
6    64.51579546     Pretest       6      2     2   BSTfemale    1
7    63.25669342     Pretest       7      2     2   BSTfemale    1
8   109.09354856     Pretest       8      2     2   BSTfemale    1
9   140.69340502     Pretest       9      2     2   BSTfemale    1
10   93.94269807     Pretest      10      2     2   BSTfemale    1
11   43.76802256     Pretest      11      2     2   BSTfemale    1

...
898  60.85271722 FU_12_month     223      1     2     BSTmale    4
899  82.75598576 FU_12_month     224      1     2     BSTmale    4
900 -32.38497309 FU_12_month     225      1     2     BSTmale    4

Answer 1

对于所期望的问题，这个问题看起来有点模糊，但这肯定不是第一次有人因为R多元函数而难以得到他们想要的东西。我希望“多变量测试：＆＃34;＆＃34;＆＃34;＆＃34;＆＃34;＆＃34;＆＃34;在eta-squared与summary.Anova.mlm的结果中，返回Pillai，Wilks，Hotelling-Lawley和Roy测试，然后分开（multivariate=TRUE参数的值）的部分（＆＃34; Mauchly测试球形度＆＃34; ＆＃34; Greenhouse-Geisser和Huynh-Feldt对球形离开的修正＆＃34;。我不知道你是否只想从控制台会话中删除结果，或者你是否想要对结果进行进一步处理。（而后一种希望在过去曾对Rhelp提出过多次请求。）

如果需要，可以使用car :: Anova帮助页面上的代码构建一个示例以供进一步讨论（因为您没有包含可以复制的数据示例）：

multivariate

如何将aov（）重新制作汽车包Anova（）以获得Mauchly对球形度，Greenhouse-Geisser和eta-squared的测试？

1 个答案: