ロバストパラメータ設計（タグチメソッド）をRで実装する　その２　

以下に載っているマグカップの例題

制御因子について（6因子、２水準）

誤差因子
ノイズは3レベル（衝突、温度負荷、持ち運び性）、入力は2レベル（それぞれのノイズについて２つずつ）

直交表はL8

評価結果は以下

結果だけ抜き出す

6.8	7.8	7.4	7.1
4.9	5.2	5.5	5.7
6.8	7.1	7	6.9
5.3	6.1	6.2	6.3
6	6.3	6.3	6.2
5.1	5.3	5.4	5.9
5.2	5.5	5.5	5.8
6.1	6.7	6.7	6.4

あとはRで実装する。

library(qualityTools)

tdo = taguchiDesign(design="L8",randomize=F, replicates=4)
tdo
write.csv(tdo,"mug.csv",row.names = F)
# ここでmug.csvを開いて、コピペする
d =read.table("mug.csv",sep=",",header = TRUE)

response(tdo) <- d$y
snPlot(tdo,
       # type="larger",
       # type="smaller",
       type="nominal",
       # factors=c('B','C','D','E','F','G','H'),
       # points = TRUE, 
       col = "black",
       single=TRUE,
       pch = 16,
       lwd=5, 
       lty = 1,
       ylab = "Mean of S/N ratio"
)

par(mfrow = c(1,8)) #1行にする

effectPlot(tdo,
           single=TRUE,
           # points = TRUE,
           col='black',
           pch=16,
           lwd=5, 
           lty=1,
           ylab = "mean of y"
)

上記のコピペのやり方

mug.csv
NAの列に、結果をコピペしていく。
４回コピペすれば良い。

SN比

平均値

結果が一致している。