設計開発に役立つ実験計画手法～要因配置実験から応答曲面法まで～（株）日本科学技術研修所
数理部

犬伏秀生

１．設計開発と実験計画手法
設計開発においては，対象とするシステムの設計パラメータを決定したり，採用したアイデアやメカニズムを評価したりする作業が必要となる．
このとき，実験計画手法を利用することにより，設計パラメータの決定やメカニズムの評価などを効率的かつ的確に実施できることが期待できる．
実験計画手法に基づく実験は，大きく下図の流れで実行される．
実験の計画
特性・要因の決定
計画表作成

実験

効果の確認

解析
最適条件探索本資料では，JUSE-StatWorks/V4.0 を用いて実験計画法手法を適用する上で大切な「特性・要因の決定」および「解析」に関して，まずいくつかの実験計画手法の考え方を紹介し，次に JUSE-StatWorks/V4.0 の実験計画法関連機能を説明する．

２．実験計画法を成功するための特性・要因の決定
品質工学（タグチメソッド）では，品質を以下のように分類する．（［3］p.17）
分類

内容

下流の品質

市場において消費者が評価する品質．

品質に関わる活動の自動車のエンジンならば，燃費，騒音，振動などである．

目的

図面寸法やスペックなどで，工程や製品のコントロールに用いる品質特下流の品質を獲得す性．計測器で計測できる項目．

るための手段

中流の品質
上流の品質

特徴

製品が市場のいろいろな条件で目的とする機能（目的機能）をどの程度達成できるかを示す測度．
自動車のエンジンであれば，目的機能は燃焼エネルギーから機械出力（トルクなど）への変換である．

源流の品質

基本機能の機能性（機能のばらつきの程度）で，技術のロバストネスを示す測度．技術開発段階で用いる品質．
自動車のエンジンであれば，基本機能は化学反応（特に，酸素の化学反応）である．

実験を行うにあたっては，何を（どの品質を）評価するのかを決める必要がある．さらにそれは，どのような特性を計測し，解析特性（解析に使用する特性）としてどのような特性を用いるのかという問題につながる．こ
こでは，実験計画手法で使用される解析特性のいくつかを紹介する．
２－１．計測特性をそのまま解析特性として使用
計測したデータをそのまま解析特性として使用する．この場合は，品質特性（中流，下流の品質の特性）やその代用特性が計測される場合が多く，その特性は望大特性，もしくは，望小特性であることが多い．
２－２．計測特性から解析特性を計算
これには等分散性を確保するために対数変換を行う場合なども含まれるが，ここでは特に，「SN 比」，「感度」，「望ましさ」という解析特性について紹介する．

株式会社 日本科学技術研修所 数理部〒151-0051 東京都渋谷区千駄ヶ谷 5-10-9

TEL 03-5379-1465© 2010, The Institute of JUSE. All Rights Reserved.

２－２－１．SN 比，感度（基本機能，目的機能の評価）「SN 比」，「感度」は品質工学（タグチメソッド）で用いられる解析特性である．
品質工学（タグチメソッド）では，基本機能の機能性（機能のばらつきの程度）を評価することが推奨されており，基本機能を考えることにより，全ての品質問題を一度に解決できるとしている．（［1］p.79）
基本機能，目的機能は，一般的に入出力の関係で捉えられ（「動特性」と呼ばれる），その関係は次のような状態であると望ましいとされる；・入出力の関係が，設計者が理想的であるとする状態にある・誤差条件の影響を受けにくい
このような望ましい状態は，適切に設計パラメータ（制御因子の水準）を決めることにより達成されることが期待できる（達成されない場合は新たなアイデアや新たなメカニズムが必要となるかもしれない）．
例えば，基本機能（目的特性）の理想的な入出力の関係式が原点を通る直線（比例式）である場合，望ましい

77665544

出力

出力状態と望ましくない状態は下図のようになる（●：誤差条件 1 でのデータ，▲：誤差条件 2 でのデータ）．

33221100012340512345

入力

入力

望ましい状態

望ましくない状態（誤差条件の影響が少なく，直線上に乗っている）（誤差条件の影響が大きく，直線上に乗っていない）

品質工学（タグチメソッド）では，基本機能（目的機能）を望ましい状態とする設計パラメータ（制御因子の水準）の決定は，次の 2 つのステップで行う；
ステップ 1：

ばらつき（理想的な関係からのずれや誤差条件による影響）が小さくなるように設計パラメータを決める．

ステップ 2：

ばらつきにはあまり影響がないが出力の大きさに影響がある制御因子を用いて，目標への合わせ込みをおこなう．

ステップ 1 を実現するために使用される解析特性が「SN 比」であり，ステップ 2 を実現するためにしばしば用いられる解析特性が「感度」である．ただし，SN 比，感度は機能性評価のための解析特性の総称であり，具体的な定義は個別の問題毎に異なる．
例えば，基本機能（目的特性）の入出力の理想的な関係式が直線であれば，
SN 比 ＝ 出力の大きさ（Signal）とばらつき（Noise）との比
感度＝ 出力の大きさ

776665554443出力

7出力出力として定義される．このとき，入出力の状態と SN 比の大きさの関係は下図のとおりである．

3322211100012345001234501入力入力234

入力（１）傾き：大，ばらつき：少（２）傾き：小，ばらつき：少（３）傾き：中，ばらつき：大→SN 比が大きい→SN 比が小さい→SN 比が小さい© 2010, The Institute of JUSE. All Rights Reserved.

5

２－２－２．望ましさ（多特性の同時最適化）
設計パラメータ（制御因子の水準）が複数の特性と関係している場合，それらの特性が全てそれなりに望ましい状態となるように設計パラメータを決める必要がある．
この目的のための有力な方法の一つが，「望ましさ」（desirability）というスカラー量（単一の量）を解析特性として使用する方法である．このスカラー量である「望ましさ」は 0～1 の値をとり，その値が大きければ大きいほど，考えた特性全てが望ましい状態であることを表す．
「望ましさ」は，以下の手順で定義される；
手順１：特性毎に望ましさ関数を定義する．
望ましさ

望ましさ

1.0

1.0

0.00.0

特性の値
十分な値の下限許容される値の下限目標値
特性の値許容される値の上限許容される値の下限（1）望大特性に対する望ましさ関数（２）望目特性に対する望ましさ関数

手順２：特性毎の望ましさ関数を総合した望ましさ関数を定義する．
特性 y1，y2，‥，ym の望ましさ関数を d1，d2，‥，dm とするとき，それらを総合した望ましさ関数 D(x)を d1，d2，‥，dm の幾何平均として定義する（算術平均とする場合もある）．
D ( x) = m d1d 2 L d m
例えば，以下のような特性，因子を取り上げたものとする；（［2］p.202）
特性：転換量（望大特性），活性度（望目特性）
因子：反応時間，温度，触媒比率
そして，転換量，活性度の望ましさ関数を以下のように定義するものとする．
望ましさ

望ましさ

1.0

1.00.00.0
転換量80

活性度9755.057.560.0

このとき，それぞれの因子の水準に応じて，望ましさは以下のようになる；因子の水準
反応時間：x11
条件 1温度：x12

触媒比率：x13
反応時間：x21
条件 2温度：x22

触媒比率：x23
反応時間：x31
条件 3温度：x32

触媒比率：x33

特性の推定値→

転換量：81.10

推定式

活性度：59.85→

転換量：97.53

推定式

活性度：55.16→

転換量：93.86

推定式

活性度：58.59

→
望ましさ
関数
→
望ましさ
関数
→
望ましさ関数

望ましさ

望ましさ

（特性毎）

（総合）

転換量：0.0642→

活性度：0.0601

幾何平均

転換量：1.0000→

活性度：0.0640

幾何平均

転換量：0.8153→

活性度：0.5640

幾何平均© 2010, The Institute of JUSE. All Rights Reserved.

0.06210.25300.6781

これより，特定の特性値のみが完全に望ましい状態よりは，全ての特性値がそれなりに望ましい状態の方が「望ましさ」が大きくなる傾向があることが分かる．
ここで，ある特性を重視したい，ある特性はそれほど重視する必要はない，など特性にウェイト付けしたい場
合は，望ましさ関数の形状を以下のように定義する；（［4］p.250）
望ましさ

望ましさ

1.0

1.00.0

目標値

許容される値の下限

目標値0.0特性の値許容される値の上限許容される値の下限

特性の値許容される値の上限（１）重視する特性の望ましさ関数（２）それほど重視しない特性の望ましさ関数（望目特性の場合）（望目特性の場合）

望ましさ関数を（１）のように定義した特性については，その特性の値をかなり目標値に近づけないと「望ましさ」は大きくならない．逆に，（２）のように定義した特性については，その特性の値が目標値にあまり近づけなくても「望ましさ」はそれなりの大きさとなる．

３．応答曲面法による解析
どのような解析特性を使用するにせよ，解析特性の望ましい状態は設計パラメータ（制御因子の水準）を適切な値に設定することにより達成される．これには，解析特性に対する設計パラメータの因果的効果を利用する．
そして，解析特性と設計パラメータとの因果関係は，実験計画法に基づいて行われた実験で得られたデータ（実験データ）から推定することができる．
なぜならば，実験計画法での“実験”は，外的操作を加え観測する，いわゆる帰納的実験（特定の結果から一般的理論を構成する）であり（［3］p.53），また，実験計画法に基づいて実験を行うと，取り上げた要因同士の交絡（効果が分離できないこと）は実験の計画において避けられ，さらに，実験順序がランダマイズされている（無作為化が行われている）ことにより，取り上げた要因と取り上げていない要因との交絡も避けられている．
このことから，実験データから得られたモデル式は，（実験が行われた領域での）因果関係を表していると考えられる．
例えば，一元配置実験（特性：化学薬品の収率（望大特性），因子：反応温度）で下図のデータが得られたとす
る．
収率
9291908980℃

90℃

100℃

110℃

反応温度

このとき，「一元配置法」による分析では，以下の手順で分析を行う．
手順１：分散分析を行い，収率に対する反応温度の効果の有無を確認する（分散分析で有意となれば効果あり）
手順２：効果があれば，各反応温度での収率を推定し，最適条件（収率が最も高くなる反応温度）を求める．

© 2010, The Institute of JUSE. All Rights Reserved.

ここで，分散分析を行うということは，実験で取り上げた反応温度において収率の平均値に違いがあるかどうかを見ること，すなわち，
帰無仮設 H0：α1＝α2＝α3＝α4=0
帰無仮設 H1：α1＝α2＝α3＝α4=0 ではない（少なくとも一つの“=”が成立しない）の検定を行うことである．
収率9291

α3

α2

α4

α1
90

89

80℃

90℃

100℃

110℃

反応温度

もし分散分析で有意となれば，各反応温度における収率の平均値は等しくないと判断され，その違いは「反応温度を変えると収率の平均が変わる」という因果関係を示すものと判断される．
以下，因子が量的因子である場合の方法論である「応答曲面法」を紹介する．
３－１．応答曲面法の利点と方法
分散分析で効果を確認し，実験で取り上げた水準の組合せから最適条件を探すという解析方法（一元配置法，二元配置法，直交配列表など）はよく知られているが，取り上げた因子の水準の順序や間隔に意味があれば（すなわち，量的因子の場合），直線（1 次モデル）や曲線（2 次モデルなど）を当てはめるという解析方法も考えられる．
例えば，前述の一元配置実験（特性：化学薬品の収率（望大特性），因子：反応温度）のデータには二次式がよく当てはまり，下図のようになる；

92

収率

91

90

8980
回 帰式

:

Y=90100
反応温度

110

4 8 . 0 5 3 3 3 3 + 0 . 8 8 2 3 3 3 X - 0 . 0 0 4 5 0 0 X ^ 2

収率 = 88.6658 + 0.0273 × 反応温度 − 0.0045 × (反応温度 − 95)2

一般的に，直線（1 次モデル）や曲線（2 次モデルなど）を当てはめることの利点は以下のとおり；
・因子と特性との関係が把握し易くなる
・領域内の任意の水準が最適条件の候補となるため，より良い条件を得ることができる© 2010, The Institute of JUSE. All Rights Reserved.

このように，取り上げた因子が量的因子の場合に，直線（1 次モデル）や曲線（2 次モデルなど）を当てはめる方法論を「応答曲面法」（RSM：Response Surface Methodology）と呼ぶ．
ここで，１次モデルは１次項のみを含むモデル，２次モデルは１次項＋２次項＋交互作用項を含むモデルのことを指し，応答曲面法では通常１次モデルか２次モデルを当てはめる．
例えば，因子が 2 つの場合の 1 次モデル，２次モデルは以下のとおり；
１次モデル： y = b0 + b1x1 + b2 x2 ，２次モデル： y = b0 + b1x1 + b2 x2 + b11x12 + b22 x2 2 + b12 x1x2応答曲面法では，一元配置法や直交配列表などの方法論と解析方法が異なるだけではなく，使用される計画も異なり，１次モデルや 2 次モデルを当てはめることを念頭に置いた計画が使用される．
応答曲面法と従来の実験計画法との比較
方法論

主な計画

一元配置法

大まかな解析手順

要因配置計画，直交配列表

１．分散分析で効果がある要因を確認

二元配置法
直交配列表
応答曲面法

２．取り上げた水準の組み合わせの中から最など

適条件を探索＜１次モデルの当てはめ＞

１．１次モデル，２次モデルを当てはめる

2k 要因配置計画，2k 一部実施要因配置計

２．当てはまりの悪さ（Lack of Fit）に関す

画，プラケット・バーマン計画

る検定で，モデルの当てはまりを確認＜二次モデルの当てはめ＞

３．当てはめたモデルで最適条件を探索

ボックス－ベーンケン計画，中心複合計画（より詳細な解析手順は下図を参照）

First-ordermodelNo

Fit
Well

Add/Dropnewvariables

Yes

Yes
NoYou are still pretty
far away from theoptimal point

Find the direction
for improvement

You are probably
close to the
optimal setting

Second-ordermodelSet up a newcenter point

FitWellNo

Yes

Higherordermodel?

Optimization

応答曲面法の解析手順（［5］p.65）
ここで，１次モデル，２次モデルを当てはめる方法は，いわゆる回帰分析と同様である（通常，最小二乗法によりパラメータを推定する）．しかし，前述のとおり，得られた１次モデル，２次モデルは因果モデルを表すと考えられ，制御などに用いることができる．

© 2010, The Institute of JUSE. All Rights Reserved.

そして，当てはめられた１次モデル，２次モデルから，因子と特性（「応答」とも呼ばれる）との関係は通常，等高線図で表現される．この図は，２つの因子を横軸，縦軸にとり（ただし,ここでは他の因子の水準は固定される），それらの因子の水準の組合せに対する特性の様子を示す．等高線を概観することにより，曲面の形状や最適条件の位置が分かる（ただし，曲面の形状は，横軸，縦軸にとられていない因子（水準が固定されている因子）の水準に依存することに注意が必要である）．
x2
2.0x2
2.0x2
2.0

1.2

1.2

1.2

0.4

0.4

0.4

-0.4

-0.4

-0.4

-1.2

-1.2

-1.2

-2.0
-2.0

-1.2

-0.4

0.4

1.2

2.0

x1-2.0
-2.0

-1.2

-0.4

0.4

1.2

2.0

x1

-2.0
-2.0

-1.2

-0.4

0.4

1.2

2.06095

43
41

50

9039

40
85

37

30

35802033

210
1（２）２次モデル（最大値）

-2-10
10.5

-2-0.5-1-1.5002-2-1010.5-2-0.5-1-1.570

01.5-2-1.6-0.8-1.20-0.4

0.4

1.20.821.6（１）１次モデル

21225

1.22975

21.527-2-1.2-0.4
0.431（３）２次モデル（停留点）© 2010, The Institute of JUSE. All Rights Reserved.

x1

４．JUSE-StatWorks/V4.0 の実験計画法関連機能
最後に，JUSE-StatWorks/V4.0 の実験計画法関連機能について紹介する．
４－１．JUSE-StatWorks/V4.0 の実験計画法関連機能
要因配置実験

解析手法
要因配置実験のための計画

一元配置分散分析

二元配置分散分析

多元配置分散分析

直交表実験

直交表実験のための計画

直交配列表

応答曲面法

応答曲面法のための計画

１特性の最適化主な出力，機能の特徴≪主な出力≫
実験計画表≪機能の特徴≫・26 因子（水準数：2～24）まで対応
・実験順序のランダマイズ可能（分割実験にも対応）≪主な出力≫
データプロット，分散分析表，推定値プロット，推定値，差の推定，残差≪機能の特徴≫・64 水準まで解析≪主な出力≫
データプロット，分散分析表，推定値プロット，推定値，差の推定，残差≪機能の特徴≫・64 水準まで解析可能
・交互作用項を誤差へプーリング≪主な出力≫
分散分析表，推定値，推定値プロット，差の推定，残差≪機能の特徴≫・解析可能な実験種類：母数模型，混合模型，多段分割実験，枝分かれ実験，多方分割実験
・母数因子，変量因子，繰り返し，ブロック因子を解析可能
・6 因子まで解析可能（繰り返し，ブロック含む．また各因子 64 水準まで）
・分散分析表で指定した要因を誤差へプーリング
・推定式に取り込む要因を分散分析上で任意に指定可能
・測定の繰り返しがある場合の解析が可能≪主な出力≫
実験計画表≪機能の特徴≫・使用可能な直交表：L8，L16，L32，L64，L9，L27，L81
・最大 26 因子までわりつけ可能・線点図に基づくわりつけ支援機能を利用
・実験順序のランダマイズ可能（分割実験にも対応）≪主な出力≫
効果表，分散分析表，推定値，推定値プロット≪機能の特徴≫・解析可能な直交表：L8，L16，L32，L64，L9，L27，L81
・分割実験の解析が可能
・多水準作成法（2 水準系直交表に 4 水準因子のわりつけ）に対応
・測定の繰り返しがある場合の解析が可能
・分散分析表で指定した要因を誤差へプーリング
・分散分析表で寄与率を表示可能
・推定式に取り込む要因を分散分析上で任意に指定可能≪主な出力≫
実験計画表≪機能の特徴≫・使用可能な計画は以下のとおり；
分類
含まれる計画2 水準系
2^k 要因実験，2^(k-p)一部実施要因実験，プラケッ
ト・バーマン計画，L8，L16，L32，L64
3 水準系
3^k 要因実験，3^(k-p)一部実施要因実験，ボックス－ベーンケン計画，L9，L27，L18
中心複合計画中心複合計画
D－最適計画
D－最適計画≪主な出力≫
変数選択，応答曲面，停留点，残差の検討，選択履歴，回帰係数，カテゴリスコア，分散分析表，予測≪機能の特徴≫・指定可能な説明変数は 1～253（量質混在可）．
・「変数の指定」ダイアログで，量的因子の交互作用項，2 次項を作成可能
・中心化（交互作用項，2 次項を偏差で作成）可能© 2010, The Institute of JUSE. All Rights Reserved.

多特性の最適化

品質設計（ﾀｸﾞﾁﾒｿｯﾄﾞ）

パラメータ設計のための計画

パラメータ設計

許容差解析

効果プロット

マルチバリチャート・応答曲面を出力可能（横軸，縦軸にとっていない因子の水準を動的に変更することも可能）
・分散分析表に当てはまりの悪さ（Lack of Fit）を出力
・手動，もしくは，逐次変数選択法（変数増減法等）による変数選択可能≪主な出力≫
係数一覧，変数選択，選択履歴，最適化グラフ，回帰係数，分散分析表，予測，残差の検討≪機能の特徴≫・「変数の指定」ダイアログで，量的因子の交互作用項，2 次項を作成可能
・指定可能な特性は 1～6，因子は 1～8（量的因子のみ）
・最適化グラフで対話的に最適条件を探索
・分散分析表に当てはまりの悪さ（Lack of Fit）を出力・特性毎に変数選択，残差の検討が可能≪主な出力≫
実験計画表≪機能の特徴≫・使用可能な直交表：L8，L16，L32，L64，L9，L27，L81，L12，L18，L36(3^13)，L36(2^11×3^12)・最大 26 因子までわりつけ・外側因子の設定（信号因子，誤差因子，繰り返しの水準数等）が可能≪主な出力≫
要因効果図（SN 比，感度），分散分析表（SN 比，感度），推定値プロット（SN 比，感度），推定値（SN 比，感度）≪機能の特徴≫・解析可能な直交表：L8，L16，L32，L64，L9，L27，L81，L12，L18，L36(3^13)，L36(2^11×3^12)
・多水準作成法，ダミー法に対応
・解析可能な特性の種類は以下のとおり；
動特性：比例式，一次式，基準点比例式
静特性：望目特性，望大特性，望小特性，ゼロ望目特性
・既に SN 比，感度が計算済みの場合も解析可能
・SN 比，感度の要因効果図を一画面に表示比較
・SN 比，感度の利得を出力可能≪主な出力≫
分散分析表≪機能の特徴≫・解析可能な直交表：L12，L18，L36(3^13)・目標値の入力が可能
・ダミー法に対応・変動の分解が可能≪主な出力≫
交互作用図（因子の全ての組合せ），特性値グラフ（因子の主効果）≪機能の特徴≫・指定可能な因子は 1～255（量質混在可）
・一覧から注目するグラフを拡大表示することが可能≪主な出力≫
マルチバリチャート≪機能の特徴≫・指定可能な因子は 1～4（質的因子のみ）
・因子の順序の入れ替えなどが可能

４－２．今後の機能強化について
弊社では，引き続き JUSE-StatWorks の機能強化を行うが，特に実験計画法関連機能については今後以下の機能強化を検討している；
解析手法
直交配列表
１特性の最適化
多特性の最適化

パラメータ設計

機能強化内容・ダミー法への対応
・最適条件の探索機能強化（多特性の最適化も同様）
・質的変数への対応
・解析可能な因子数の拡張
・複数の応答曲面の重ね合わせ
・標準 SN 比への対応© 2010, The Institute of JUSE. All Rights Reserved.

５．参考文献［1］

田口玄一，（2005），ハード製品の研究開発における戦略

目的機能と基本機能，標準化と品質管理，

58，72-83［2］

山田秀，（2004），実験計画法－方法編－，日科技連［3］

宮川雅巳，（2000），品質を獲得する技術，日科技連［4］Myers, R. H. and Montgomery, D. C.，（1995），Response Surface Methodology，John Wiley & Sons.

［5］

DENNIS K. J. LIN，（1999），Discussion，Journal of Quality Technology，31，61-66［6］

日本科学技術研修所，（2004），JUSE-StatWork/V4.0 ユーザーズマニュアル© 2010, The Institute of JUSE. All Rights Reserved.

掲載されている著作物の著作権については，制作した当事者に帰属します．
著作者の許可なく営利・⾮営利・イントラネットを問わず，本著作物の複製・転⽤・販売等を禁⽌します．
所属および役職等は，公開当時のものです．

■公開資料ページ
弊社ウェブページで各種資料をご覧いただけます http://www.i-juse.co.jp/statistics/jirei/■お問い合わせ先(株)⽇科技研 数理事業部 パッケージサポート係 http:/www.i-juse.co.jp/statistics/support/contact.html© 2010, The Institute of JUSE. All Rights Reserved.