直交表L8を使ったパラメータ設計がわかる
「直交表を使ったロバストパラメータ設計ができない」などと困っていませんか？
こういう疑問に答えます。
本記事のテーマ
おさえておきたいポイント
- ①パラメータ設計にはどの直交表を使っても良い
- ➁実験計画法の復習
- ➂直交表L8を使ったパラメータ設計事例
- ➃SN比と感度の計算
- ➄最適条件の選定
ロバストパラメータ設計
タグチメソッド
手法に溺れるな！
数式と理論で理解しよう！
品質工学、ロバストパラメータ設計、タグチメソッド
結局わからない！
⇒ＱＣプラネッツが解決！
①パラメータ設計にはどの直交表を使っても良い
品質工学、ロバストパラメータ設計、タグチメソッドで扱う直交表について疑問に思うことが2つあります。
- 「品質工学＝混合系直交表」じゃないとダメなのか？
- 「品質工学≠実験計画法」は正しいのか？
「品質工学＝混合系直交表」は正しいのか？
品質工学、ロバストパラメータ設計、タグチメソッドって、最初から、特殊な直交表L12,L18が出て来ます。その理由は
でも、これが意味わからないんですよ！
データが決めるはずで、
我々ではないから
データは様々な要因が絡み合って数値化されます。であれば、ナチュラルに主効果、交互作用を含めて検討しても良いと考えます。
「品質工学≠実験計画法」は正しいのか？
同じ直交表を使うのに、
これもピンと来ません。同じで良いではないかと！
実験計画法の延長として品質工学を検討したり
一般の直交表を使ってパラメータ設計してみましょう。
で、おかしい結果になるか確かめてみましょう。
解析結果は
データが決めるはずで、
我々ではないよ！
品質工学に一般の直交表を使ってみる
ＱＣプラネッツは長年にわたりこの疑問を抱いていましたので、実際に
●L8(本記事)
●L16
●L27
と混合系の
●L12
●L18
を取り上げてみます。
➁実験計画法の復習
ＱＣプラネッツは
品質工学＝実験計画法 という考えなので、
まず実験計画法を復習しましょう。
関連記事で、確認ください。
 |
究める！実験計画法 QCプラネッツが解説する究める実験計画法。多くの教科書がある中、勉強してもどうしても分からない、苦労している難解な箇所をすべて解説します。多元配置実験、乱塊法、分割法、直交表などなど多くの手法を個別に公式暗記せず、データの構造式をみればすべて導出できる新しい実験計画法を解説します。 |
➂直交表L8を使ったパラメータ設計事例
事例
次の問いを考えます。
あるデータから最適条件を直交表L8を使って求めたい。
(1) 分散分析表を作れ。
(2) 各因子の、各水準におけるＳＮ比と感度を計算し、要因効果図を作れ。
(3) 最適条件を選び、その条件におけるＳＮ比と感度を計算せよ。
直交表L8
各因子の平方和と分散分析を解析
直交表L8は2水準系なので、各列の平方和を計算する公式があります。関連記事で解説しています。
 |
【本記事限定】直交表の各列の平方和の式は自力で導出できる【必見】 直交表の各列の平方和を導出する方法を知っていますか？公式暗記で済ませていませんか？本記事では、実験計画法の直交表の各列の平方和を導出する方法を詳しく解説します。本記事しか書いていない、直交表の知見を広げたい方は必見です。 |
これをもとに各列の平方和を計算すると、下表になります。
●分散分析表
| L8 | 効果 | 平方和S | 自由度φ | 平均平方V |
| 1 | A | 60.5 | 1 | 60.5 |
| 2 | B | 24.5 | 1 | 24.5 |
| 3 | A×B | 128 | 1 | 128 |
| 4 | C | 32 | 1 | 32 |
| 5 | A×C | 24.5 | 1 | 24.5 |
| 6 | B×C | 12.5 | 1 | 12.5 |
| 7 | e | 18 | 1 | 18 |
| 15 | e | 18 | 1 | 18 |
| 合計 | T | 300 | 7 | – |
データの構造式
分散分析を扱うための最も重要なデータの構造式を定義します。今回は全列を成分に合わせた効果とするので、
\(x_{ijk}\)=\(μ\)+\(a_i\)+\(b_j\)+\(c_k\)
+\(ab_{ij}\)+\(ac_{ik}\)+\(bc_{jk}\)
+\(e_{ijk}\)
(\(i,j,k=1,2\))
➃SN比と感度の計算
SN比と感度Sの公式
関連記事にも公式導出過程を解説しています。
 |
品質工学のSN比が導出できる 品質工学のSN比の式 η=10log (Sm-Ve)/Veがちゃんと導出できますか？ 本記事はSN比を導出します。公式暗記に頼らず、式変形から意味を理解して、式を使うようにしましょう。 |
●SN比ηは
η=10\(log \frac{μ^2}{σ^2}\)=10\(log \frac{\frac{1}{n}(S_m-V_e)}{V_e}\)●感度Sは
S=10\(log μ^2\)=10\(log \frac{1}{n}(S_m-V_e)\)
ですが、今回簡略化のため、
●SN比ηは
η=\(\frac{μ^2}{σ^2}\)=\(\frac{\bar{x}^2}{s^2}\) (\(s\)は標準偏差)
●感度Sは
S=\(μ^2\)=\(\bar{x}^2\)
で計算します。
各効果の各水準における平均\(\bar{x}\)と標準偏差\(s\)
各効果の水準１，２に属するデータの平均と標準偏差を計算すると下表になります。
| – | 水準1 | 水準2 | |||
| No | 列 | 平均\(\bar{x}\) | 標準偏差\(s\) | 平均\(\bar{x}\) | 標準偏差\(s\) |
| 1 | A | 17.88 | 8.54 | 20.63 | 2.63 |
| 2 | B | 18.38 | 6.34 | 20.13 | 7.18 |
| 3 | A×B | 18.25 | 5.97 | 20.25 | 7.33 |
| 4 | C | 21.25 | 5.45 | 17.25 | 5.26 |
| 5 | A×C | 20.13 | 5.5 | 18.38 | 7.85 |
| 6 | B×C | 18.63 | 5.56 | 19.88 | 8.06 |
| 7 | e | 18.5 | 5.16 | 20 | 8.21 |
例えば、因子Aにおいて、
●水準１：35,29,48,31
●水準２：44,39,43,39
なので、それぞれ4個の平均と標準偏差を計算します。
それを因子Aから誤差eまでの7列分を計算します。
直交表L8の各列のSN比と感度
●SN比ηは
η=\(\frac{μ^2}{σ^2}\)=\(\frac{\bar{x}^2}{s^2}\) (\(s\)は標準偏差)
●感度Sは
S=\(μ^2\)=\(\bar{x}^2\)
で計算すると、下表になります。
| – | 水準1 | 水準2 | |||
| No | 列 | SN比 η | 感度S | SN比 η | 感度S |
| 1 | A | 4.38 | 319.52 | 61.5 | 425.39 |
| 2 | B | 8.39 | 337.64 | 7.85 | 405.02 |
| 3 | C | 9.34 | 333.06 | 7.64 | 410.06 |
| 4 | A×B | 15.22 | 451.56 | 10.76 | 297.56 |
| 5 | A×C | 13.39 | 405.02 | 5.48 | 337.64 |
| 6 | B×C | 11.22 | 346.89 | 6.08 | 395.02 |
| 7 | e | 12.83 | 342.25 | 5.94 | 400 |
よく 対数を取ってSN比や感度Sの値を計算しますが、別になくてもOKなので、対数にしていません。
要因効果図があると見やすいですが、数表からも確認できるので、割愛します。
➄最適条件の選定
ここで、因子A,Cの水準の高い方を選択します。
\(μ(A_i C_k)\の式を先に作ります。
関連記事で解説しています。
 |
【簡単】データの構造式で実験計画法がわかる(必読) 実験計画法が難しい、多元配置実験、乱塊法、分割法、などたくさんの手法を学ぶのが大変など困っていませんか？本記事では、データの構造式さえ理解すれば実験計画法がすぐマスタできるように、わかりやすく解説します。早く実験計画法をマスターした方は必見です。 |
\(μ(A_i C_k)\)
=\(μ\)+\(a_i\)+\(c_k\)+\(ac_{ik}\)
=\(\bar{\bar{x}}\)+(\(\bar{x_{i・・}}-\bar{\bar{x}}\))+(\(\bar{x_{・・k}}-\bar{\bar{x}}\))
+\((\bar{x_{i・k}}-\bar{x_{i・・}}-\bar{x_{・・k}}+\bar{\bar{x}})\)
=\((\bar{x_{i・k}}-\bar{\bar{x}}\))
=(\(\widehat{μ+a_i+c_k}\))-\(μ\)
とすっきりした式になります。
暗記不要で、データの構造式からどんな組み合わせパターンも式が作れます！
●SN比において、
A,Cで値のSN比が大きい水準をみると
AC_1なので、
\(μ(A_ i C_k)\)= (\(\widehat{μ+ac_1}\))-\(μ\)
=13.39-12.85=0.53
●感度において、
A,Cで値の感度が大きい水準をみると
AC_1なので、
\(μ(A_ i C_k)\)= (\(\widehat{μ+ac_1}\))-\(μ\)
=405.02-371.90=33.11
と計算できました。
直交表L8を使って、SN比、感度の計算を実施しました。
まとめ
「直交表L8を使ったパラメータ設計がわかる」を解説しました。
- ①パラメータ設計にはどの直交表を使っても良い
- ➁実験計画法の復習
- ➂直交表L8を使ったパラメータ設計事例
- ➃SN比と感度の計算
- ➄最適条件の選定
Warning: count(): Parameter must be an array or an object that implements Countable in /home/qcplanets/qcplanets.com/public_html/wp-content/themes/m_theme/sns.php on line 119