myscratchpad

#aovの結果
Error: a:s
          Df  Sum Sq Mean Sq F value    Pr(>F)   
a          3 217.500  72.500 21.4437 1.347e-06 ***
s          7  77.000  11.000  3.2535   0.01682 * 
Residuals 21  71.000   3.381 

# Anovaの結果
                SS num Df Error SS den Df         F       Pr(>F)
(Intercept) 3444.5      1       77      7 313.13636 4.537639e-07
afact        217.5      3       71     21  21.44366 1.347309e-06

Anova関数の被験者要因は
平方和: (Intercept) のError SS = 77
平均平方和: (Intercept) のError SS/ (Intercept) のden DF 77/7 = 11
誤差平方和: 要因のError SS/要因のden Df = 71/21 = 3.380952
F値: 平均平方和/誤差平方和 = 11/3.380952 = 3.253521
要因自由度: Interceptのden Df = 7
誤差自由度: 要因のden Df = 21


## aov結果のとりだし
aovrestb <- summary(aov(result ~ a + s + Error(a:s), RB))
## aovrestb[[1]][[1]]["Sum Sq"]