２要因ANOVA混合計画 by R

大学院統計で２要因ANOVA（一要因に対応があり、一要因に対応なし）の例をSPSSの反復測定を使わずに試してもらった。以下の例は、RでSPSSの例（下記）と比較できるようにしたもの。データも共通。

cf.
1つの被験者間要因と 1つの被験者内要因 （1変量）（星野祐司氏 ( 立命館大学））
http://www.psy.ritsumei.ac.jp/~hoshino/spss/anova11x.html


> library(foreign)
> x<-read.spss("anova11x.sav")
> x
$CONDTION
[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
[49] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
[97] 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3

$INTERVAL
[1] 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6
[49] 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6
[97] 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6

$SUBJECT
[1] 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 5 5 5 5 5 5 6 6
[33] 6 6 6 6 7 7 7 7 7 7 8 8 8 8 8 8 9 9 9 9 9 9 10 10 10 10 10 10 11 11 11 11
[65] 11 11 12 12 12 12 12 12 13 13 13 13 13 13 14 14 14 14 14 14 15 15 15 15 15 15 16 16 16 16 16 16
[97] 17 17 17 17 17 17 18 18 18 18 18 18 19 19 19 19 19 19 20 20 20 20 20 20 21 21 21 21 21 21 22 22
[129] 22 22 22 22 23 23 23 23 23 23 24 24 24 24 24 24

$ACTIVITY
[1] 150 44 71 59 132 74 335 270 156 160 118 230 149 52 91 115 43 154 159 31 127 212 71 224
[25] 159 0 35 75 71 34 292 125 184 246 225 170 297 187 66 96 209 74 170 37 42 66 114 81
[49] 346 175 177 192 239 140 426 329 236 76 102 232 359 238 183 123 183 30 272 60 82 85 101 98
[73] 200 271 263 216 241 227 366 291 263 144 220 180 371 364 270 308 219 267 497 402 294 216 284 255
[97] 282 186 225 134 189 169 317 31 85 120 131 205 362 104 144 114 115 127 338 132 91 77 108 169
[121] 263 94 141 142 120 195 138 38 16 95 39 55 329 62 62 6 93 67 292 139 104 184 193 122

> attach(x)
> c<- factor(CONDTION)
> i<- factor(INTERVAL)
> s<- factor(SUBJECT)

> aov.ex=aov(ACTIVITY~c*i+Error(s:c+s:i:c),x)
> summary(aov.ex)

Error: s:c
Df Sum Sq Mean Sq F value Pr(>F)
c 2 285815 142908 7.8006 0.002928 **
Residuals 21 384722 18320
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Error: s:c:i
Df Sum Sq Mean Sq F value Pr(>F)
i 5 399737 79947 29.8524 <2.2e-16 *** c*i 10 80820 8082 3.0178 0.002164 **
Residuals 105 281199 2678
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

なお、SPSSで反復測定を用いた例は以下。
http://www.psy.ritsumei.ac.jp/~hoshino/spss/anova11.html