2.今のところの結果
2.1.直交表とデータ
データは,鷲尾泰俊さんの「実験計画法入門 改訂版」,1997,ISBN:4-542-50330-5, P193 例10.3. より取っています.
2.2.ANOVA表
微妙に「実験計画法入門 改訂版」と結果がずれています.
2.3 要因効果図
要因効果図も出来ます.
3.上記の結果の出力プログラム(直交表のところだけ)
2.1.直交表とデータ
データは,鷲尾泰俊さんの「実験計画法入門 改訂版」,1997,ISBN:4-542-50330-5, P193 例10.3. より取っています.
f01 | f02 | f03 | f04 | f05 | f06 | f07 | f08 | f09 | f10 | f11 | f12 | f13 | y | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 64.20 |
2 | 1 | 1 | 1 | 1 | 2 | 2 | 2 | 2 | 2 | 2 | 2 | 2 | 2 | 64.60 |
3 | 1 | 1 | 1 | 1 | 3 | 3 | 3 | 3 | 3 | 3 | 3 | 3 | 3 | 65.10 |
4 | 1 | 2 | 2 | 2 | 1 | 1 | 1 | 2 | 2 | 2 | 3 | 3 | 3 | 64.50 |
5 | 1 | 2 | 2 | 2 | 2 | 2 | 2 | 3 | 3 | 3 | 1 | 1 | 1 | 65.00 |
6 | 1 | 2 | 2 | 2 | 3 | 3 | 3 | 1 | 1 | 1 | 2 | 2 | 2 | 63.90 |
7 | 1 | 3 | 3 | 3 | 1 | 1 | 1 | 3 | 3 | 3 | 2 | 2 | 2 | 65.10 |
8 | 1 | 3 | 3 | 3 | 2 | 2 | 2 | 1 | 1 | 1 | 3 | 3 | 3 | 64.10 |
9 | 1 | 3 | 3 | 3 | 3 | 3 | 3 | 2 | 2 | 2 | 1 | 1 | 1 | 65.00 |
10 | 2 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 64.70 |
11 | 2 | 1 | 2 | 3 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 64.30 |
12 | 2 | 1 | 2 | 3 | 3 | 1 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 64.60 |
13 | 2 | 2 | 3 | 1 | 1 | 2 | 3 | 2 | 3 | 1 | 3 | 1 | 2 | 64.30 |
14 | 2 | 2 | 3 | 1 | 2 | 3 | 1 | 3 | 1 | 2 | 1 | 2 | 3 | 64.50 |
15 | 2 | 2 | 3 | 1 | 3 | 1 | 2 | 1 | 2 | 3 | 2 | 3 | 1 | 64.30 |
16 | 2 | 3 | 1 | 2 | 1 | 2 | 3 | 3 | 1 | 2 | 2 | 3 | 1 | 63.90 |
17 | 2 | 3 | 1 | 2 | 2 | 3 | 1 | 1 | 2 | 3 | 3 | 1 | 2 | 64.90 |
18 | 2 | 3 | 1 | 2 | 3 | 1 | 2 | 2 | 3 | 1 | 1 | 2 | 3 | 64.10 |
19 | 3 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 65.00 |
20 | 3 | 1 | 3 | 2 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 62.90 |
21 | 3 | 1 | 3 | 2 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 2 | 1 | 64.80 |
22 | 3 | 2 | 1 | 3 | 1 | 3 | 2 | 2 | 1 | 3 | 3 | 2 | 1 | 64.20 |
23 | 3 | 2 | 1 | 3 | 2 | 1 | 3 | 3 | 2 | 1 | 1 | 3 | 2 | 65.30 |
24 | 3 | 2 | 1 | 3 | 3 | 2 | 1 | 1 | 3 | 2 | 2 | 1 | 3 | 64.00 |
25 | 3 | 3 | 2 | 1 | 1 | 3 | 2 | 3 | 2 | 1 | 2 | 1 | 3 | 64.80 |
26 | 3 | 3 | 2 | 1 | 2 | 1 | 3 | 1 | 3 | 2 | 3 | 2 | 1 | 64.80 |
27 | 3 | 3 | 2 | 1 | 3 | 2 | 1 | 2 | 1 | 3 | 1 | 3 | 2 | 64.00 |
2.2.ANOVA表
微妙に「実験計画法入門 改訂版」と結果がずれています.
Df | Sum Sq | Mean Sq | F value | Pr(> F) | |
---|---|---|---|---|---|
f01 | 2 | 0.24 | 0.12 | 1.93 | 0.2071 |
f02 | 2 | 0.03 | 0.01 | 0.23 | 0.7995 |
f05 | 2 | 0.05 | 0.02 | 0.37 | 0.7009 |
f09 | 2 | 2.75 | 1.37 | 21.88 | 0.0006 |
f13 | 2 | 0.51 | 0.25 | 4.04 | 0.0614 |
f01:f02 | 4 | 0.33 | 0.08 | 1.31 | 0.3444 |
f02:f05 | 4 | 2.58 | 0.65 | 10.29 | 0.0030 |
Residuals | 8 | 0.50 | 0.06 |
2.3 要因効果図
要因効果図も出来ます.
3.上記の結果の出力プログラム(直交表のところだけ)
# Test2.1 for latin Table # Cite from # Washio, Yasutoshi,実験計画法入門 改訂版, # "Introduction of Eexperimentation Design", 2nd. Ed.,1997, # ISBN:4-542-50330-5, P193 ex. 10.3. # library("xtable") y<-c(64.2,64.6,65.1,64.5,65.0,63.9,65.1,64.1,65.0,64.7,64.3,64.6,64.3,64.5,
64.3,63.9,64.9,64.1,65.0,62.9,64.8,64.2,65.3,64.0,64.8,64.8,64.0) res1<-latin.table.L27(y)
res2<-summary(aov(data=res1,formula=y ~ f01+f02+f05+f09+f13+f01:f02+f02:f05))
res1<-xtable(res1)
print(res1,type="html") res2<-xtable(res2)
print(res2,type="html")