R-部分eta平方重復測量ANOVA（汽車包裝）

Question

我有一個雙向重復測量設計（3 x 2），我想弄清楚如何計算效果大小（部分eta平方）。

我有一個矩陣，其中包含數據（稱為a），如此（重復測量）

         A.a          A.b           B.a        B.b          C.a           C.b
1        514.0479     483.4246      541.1342   516.4149     595.5404      588.8000
2        569.0741     550.0809      569.7574   599.1509     621.4725      656.8136
3        738.2037     660.3058      812.2970   735.8543     767.0683      738.7920
4        627.1101     638.1338      641.2478   682.7028     694.3569      761.6241
5        599.3417     637.2846      599.4951   632.5684     626.4102      677.2634
6        655.1394     600.9598      729.3096   669.4189     728.8995      716.4605

idata =

    Caps    Lower
       A       a
       A       b
       B       a
       B       b
       C       a
       C       b

我知道如何使用汽車包重復測量ANOVA（3型SS在我的領域是標准的，雖然我知道它會導致邏輯錯誤..如果有人想向我解釋這就像我5我會喜歡了解它）：

summary(Anova(lm(a ~ 1),
          idata=idata,type=3, 
          idesign=~Caps*Lower)),
    multivariate=FALSE)

我想我想要做的是將這部分摘要打印出來：

單變量III型重復測量ANOVA假設球形度

                     SS num Df Error SS den Df        F    Pr(>F)    
(Intercept)     14920141     1   153687      5 485.4072 3.577e-06 ***
Caps            33782        2     8770     10  19.2589  0.000372 ***
Lower           195          1    13887      5   0.0703  0.801451    
Caps:Lower      2481         2      907     10  13.6740  0.001376 ** 

並用它來計算部分ETA平方。 所以，如果我沒有弄錯的話，我需要從第一列中取出SS並將其除以（本身+該行的SS錯誤）每個效果。 這是正確的方法嗎？ 如果是這樣，我該怎么辦？ 我無法弄清楚如何從摘要打印中引用值。

Answer 1

可以使用heplots包中的etasq函數計算部分eta平方

library(car)
mod <- Anova(lm(a ~ 1),
idata = idata,
type = 3,
idesign = ~Caps*Lower)

mod

library(heplots)
etasq(mod, anova = TRUE)

既然你在詢問計算：
來自?etasq ：'對於單變量線性模型，經典η^ 2 = SSH / SST和部分η^ 2 = SSH /（SSH + SSE）。 這些在單向設計中是相同的。'

如果您希望檢查具有class的模型的η^ 2的計算代碼，則可以使用getS3method(f = "etasq", class = "Anova.mlm") 。