r - R-偏乘方平方，用于重复测量方差分析(汽车包装)

我有一个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

我知道如何对汽车包装进行重复测量方差分析(3型SS在我的领域中是标准配置，尽管我知道这会导致逻辑错误。.如果有人想像我5岁那样向我解释这一点，我会很喜欢了解它):

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

我认为我想做的是将摘要的这一部分打印出来:

单变量III型重复测量方差分析假设球形度

                     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 Error)。这是正确的方法吗？如果是这样，我该怎么办？我不知道如何从摘要打印出来引用值。

最佳答案

可以使用etasq包中的heplots函数来计算部分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")。

关于r - R-偏乘方平方，用于重复测量方差分析(汽车包装)，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/17961349/