https://seinzumtode.hatenadiary.jp/entry/2025/06/25/231509

www.youtube.com

データ
Rice Fish Meat Km distance
200 100 150 2.50
200 125 180 2.85
200 180 195 3.00
250 100 180 3.20
250 125 195 3.70
250 180 150 3.90
300 100 195 2.50
300 125 150 2.60
300 180 180 2.57

効果要因図

RのqualityToolsと比較

Replicate=1(反復なし)だとSN plotが表示されなかった。
全く同じ値で反復すると値がInfになってしまうので、小さい値（０.001）を加えて2反復とした。

# install.packages("devtools")
library(devtools)
# install_github("cran/qualityTools")
if(!require(tidyr)){install.packages("tidyr")}
if(!require(dplyr)){install.packages("dplyr")}
library(qualityTools)

tdo = taguchiDesign(design="L9",randomize=F,replicates=2)
write.csv(tdo,"L9.csv",row.names = F)

#ここでyの値を手動で入力する

d =read.table("L9.csv",sep=",",header = TRUE)
y=d$y
response(tdo)=y

snPlot(tdo,
       #type="larger",
       # type="smaller",
       type="nominal",
       factors=c('A','B','C'),
       # points = TRUE, 
       col = "black", 
       single=TRUE,
       pch = 16,
       lwd=5, 
       lty = 1,
       ylab = "Mean of S/N ratio"
)
par(mfrow = c(1,3)) #1行にする

effectPlot(tdo,
           single=TRUE,
           factors=c('A','B','C'),
           # points = TRUE,
           col='black',
           pch=16,
           lwd=5, 
           lty=1,
           ylab = "mean of y"
)
par(mfrow = c(1,3)) #1行にする

定性的にはあっているのだが、特にSN比の値はMinitabとかけ離れている。
（どのくらい誤差を与えるかに依存してしまう）