R- Частичная эта квадрат для повторных измерений ANOVA (автомобильный пакет)

Question

У меня есть двухэтапный дизайн повторных измерений (3 x 2), и я хотел бы получить представление о том, как рассчитать размеры эффекта (частичное значение в квадрате).

У меня есть матрица с данными в нем (так называемый), как это (повторные измерения)

         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 Error для этой строки) для каждого эффекта. Это правильный путь? Если да, то как мне это сделать? Я не могу понять, как ссылаться на значения из краткой распечатки.

Спасибо

Answer 1

Частичная эта-квадрат может быть рассчитана с помощью etasq функция в heplots пакет

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). Они идентичны в одностороннем дизайне.

Если вы хотите проверить код для расчетов η ^ 2 для модели с class как в примере, вы можете использовать getS3method(f = "etasq", class = "Anova.mlm"),