2. 実験の計画

ツールバーで空白シートアイコン   をクリックするか、“New Design” をクリックします。画面左のStandard Designs の下の、“Response Surface”をクリックしてツリーを展開し、Split-Plot セクションから “Central Composite”（CCD）を選択します。因子の数として “4” をドロップダウンリストから選択します。

4因子の分割法CCD計画

 

因子 a の名称に「Gap」と入力します。なお、この因子は既に変更困難な因子にラベル付けされていますが、これは分割法を選択したからです。因子 b をクリックして2番目の因子名として「Deflection」と入力します。Tab キーで Change 列に移動したら、「h」と入力し、これを Hard に変更します。ドロップダウンリストをクリックして変更することもできますが、直接入力するほうが簡単です。

ドロップダウンリストを使って因子bをHardに変更

 

Low と High の水準の因子は入力する必要はありません。この例の著者からは、コード化された（-1 から +1）の水準しか提供されていないからです（実測値はありません）。恐らく情報を保護する必要があったのでしょう。次に、因子 C をクリックして、因子名は「Angle」と入力します。さらに、因子 D の名前として「Reynolds」と入力します。これらは ETC  因子のため、Change 列はデフォルトの Easy のままにします。

4因子の名前すべてを入力

 

画面中央下にある “Options” ボタンをクリックして、Alpha の値を “Practical” に変更します。

Alpha値のオプションをPracticalに変更

 

この値は CCD 内の軸点（星印）をどのレベルで実行するかを決定します。Alpha は計画の中心からの（コード化された単位での）距離です。

※ 分割法 CCD の構成について:この CCD の 3因子の配置を図で示すと以下のようになります。核となる要因が立方体の頂点として構成されており、コード化された２つの構成単位（上記に記載した -1 から 1）が向かい合う面の距離となります。星印は軸点（axial point）を表します。中心からこの軸点をどれだけ遠ざけるべきかについては、統計学者間で大いに議論される問題です。この距離は、コード化された因子の水準を尺度とする “Alpha” で指定します。これから順次見て行きますが、Design-Expert にはアルファに関する様々なオプションが用意されています。アルファのデフォルト “Rotatable” は、予測内における標準誤差の分布（不確実性）を示すのに最も優れています。そのため、設定はデフォルトのままが好ましいですが、変更していけないわけではありません。回転可能なアルファには、中心からさらに遠くなるような星印のポイントが必要です（この場合は、コード化された単位で 2.0）。そのため、いつでも実行可能なわけではありません。星印のポイントの水準が中心から離れすぎていると、測定可能な結果が得られない場合があります。この理由から、アルファに Practical を使用することで星印を少し引き寄せて、実験を確実に試行することで有益な結果を得ることができます。Rotatable を指定すると因子が多くなったときに扱いにくくなるため、5 因子（k > 5）より多い場合は Alpha を Practical にすることをお勧めします。

3因子の中心複合計画

分割法を用いた中心複合計画においては、アルファ水準に加えてグループの複雑な関係が存在しています。ページの下部、Options ボタンのそばを見ると、この計画がどのようにグループ分けされたのか分かります。

グループに分解した分割法を用いたCCD計画

それぞれの種類のポイントで、グループの番号と各グループの試行番号が表示されます。これらは Options ダイアログボックスで変更することもできますが、この例ではデフォルトのままにしておきます。

 

“Next” をクリックします。R1 の Name に応答として「Lift」と入力します。

応答名としてLiftと入力

 

“Finish” をクリックして、design-building ウィザードを終了します。グループ間の因子の水準をリセットするように警告が表示されます。とりあえずこの警告は無視して、“OK” をクリックします。すると、Gap（因子 a）と Deflection（因子 b）にグループ分けされた無作為の試行順で、実験を達成するために必要とされる試行が表示されます。計画は無作為化されるため（制約なし）のため、ほとんどの場合、試行は各ユーザーによって異なります。また Design-Expert では、列のヘッダーで分かるように HTC が小文字で表示されるという決まりが採用されていることにも注意してください。こうして、HTC は ETC と区別され、これら ETC は完全無作為化された計画でいつも表示されるように大文字で示されます。

計画における最初の10試行（ユーザーにより異なる）

 

“whole-plot” を構成する HTC 因子の組み合わせを把握するために、Group 列が追加されています。グループは黒い横線で区切られている点に注目してください。どのグループの試行も HTC 因子の水準は一定に保たれています。例えば、Group 1 では、どの試行も Gap が 0 に設定され、Deflection が 1.41421 に設定されています。Group 2 の４つの試行については、すべて gap が mid-level（0）に設定され、因子 B も mid-level（0）に設定されています。

有効桁数を減らすには Deflection の列の上で右クリックし、“Edit Info” を選択します。Format ドロップダウンリストから “0.00” を選択すると小数点 2 桁を残すことができます。４因子すべてに対してこの操作を行って、計画をより見やすく、実行しやすくします。

上記因子のFactor設定を小数点2桁に設定

 

HTC 因子をグループ化するとさらに実験を試行しやすくなります。「はじめに」で述べたように Gap と Deflection を変更するには、密閉された試験機を試行ごとに開ける必要があります。分割法において、この中断が発生するのはグループ間でのみになります。この方法によって、航空エンジニアたちは試験機を開いて、再密閉し、最後に再加圧するといった時間のかかる一連の手順を踏むことなく、実験を4回行うことができます。しかし、HTC 因子の水準がたとえ変わらなくても、グループ間のすべての因子はリセットする必要があるということは忘れないでください。前に Yes をクリックしたときに、これを示唆するような警告が表示されたことを思い出してください。例えば、この計画で Groups 2 と 3 において因子 “a” と “b” は 0 で一定です。グループ間で実際に発生している変動をとらえるため、実験者は試験機を開けてこれらの因子の公称値を調整し、3番目のグループの試行が始まる前に 0 に戻す必要があります。こうして、これらの水準でリセットされて発生したノイズは正しく説明されます。