R-重复测量ANOVA（car包）的部分η平方

Question

我是一个有用的助手，可以翻译文本。

我有一个2路重复测量设计（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

我知道如何使用car包进行重复测量ANOVA（在我的领域中，类型3 SS是标准的，尽管我知道它会导致逻辑错误。如果有人能像我5岁一样解释一下，我会很高兴理解它）：

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

我认为我想要做的是提取摘要打印输出中的这一部分：

假定球形下的单因素三水平重复测量方差分析(Univariate Type III Repeated-Measures ANOVA Assuming Sphericity)

                     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误差）以计算每个效应。这是正确的方法吗？如果是，我该怎么做？我无法找出如何引用摘要打印输出中的值。

- luke123

---

这个问题可能在交叉验证上得到更好的回答。 - Thomas

你有检查过heplots包中的etasq函数吗？ - Henrik

CV上的重复发布已经关闭。请不要重复发布，并在您认为问题更适合其他SE网站时请求管理员迁移您的问题。 - chl

Answer 1

部分 $\eta^2$ 可以使用 heplots 包中的 etasq 函数计算。

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")。

- Henrik