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

以下に載っているマグカップの例題
https://www.amazon.co.jp/%E4%BE%8B%E9%A1%8C%E3%81%A8Excel%E6%BC%94%E7%BF%92%E3%81%A7%E5%AD%A6%E3%81%B6-%E5%AE%9F%E9%A8%93%E8%A8%88%E7%94%BB%E6%B3%95%E3%81%A8%E3%82%BF%E3%82%B0%E3%83%81%E3%83%A1%E3%82%BD%E3%83%83%E3%83%89-%E8%8F%85-%E6%B0%91%E9%83%8E/dp/4274219666

制御因子について（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比

平均値

結果が一致している。