実験計画(DOE) > カスタム計画の例 > 応答曲面計画 > 制約とカテゴリカルな因子を含む応答曲面計画
公開日: 04/01/2021

制約とカテゴリカルな因子を含む応答曲面計画

この例では、化学反応の収率を最適化する計画を作成します。実験には、3水準のカテゴリカルな因子が1つ含まれます。このうち2水準で制約があります。この例では、「フィルタで許可しない組み合わせを指定」というオプションを使って、この制約を入力します。

応答は「収率」です。因子は次のとおりです。

時間: 興味がある範囲は500~560秒です。

温度: 興味がある範囲は350~750度(ケルビン)です。

触媒: A、B、Cの3種類の触媒を試験する必要があります。

「触媒」の2水準について、次の制約を含める必要があります。

触媒がBのときは、温度は400度を超えるものとする。

触媒がCのときは、温度は650度未満であるものとする。

応答と因子の定義

1. [実験計画(DOE)]>[カスタム計画]を選択します。

2. 「応答」アウトラインで[Y]をダブルクリックし、「収率」に変更します。

目標は「収率」を最大化することなので、「目標」の設定は[最大化]のままにします。

3. 「因子」アウトラインで、「N個の因子を追加」の横のボックスに「2」と入力します。

4. [因子の追加]>[連続変数]を選択します。

5. [因子の追加]>[カテゴリカル]>[3水準]を選択します。

6. 因子の名前を、「時間」「温度」「触媒」に変更します。

7. 「時間」「値」を「500」と「560」に変更します。

8. 「温度」「値」を「350」と「750」に変更します。

9. 「触媒」「値」を「A」、「B」、「C」に変更します。

図5.46 因子設定 

Image shown here

10. [続行]をクリックします。

制約の定義

1. 「因子の制約を定義」アウトラインで、[フィルタで許可しない組み合わせを指定]を選択します。

2. 「フィルタ因子の追加」のリストから、「温度」「触媒」を選択します。

3. [追加]をクリックします。

4. 温度のスライダの上に表示されている式の「750」をクリックし、「400」に変更します。Enterキーを押します(または、どこか別の場所をクリックします)。

これにより、「温度」が400度以下の値になるような因子の設定は許可されなくなります。

5. 「触媒」の下の「B」のブロックをクリックします。

「温度」の制約と合わせて、「触媒」がBで「温度」が400度以下になるような因子の設定が許可されなくなります。

6. [OR]をクリックします。

2つ目の制約を定義できます。

7. 「フィルタ因子の追加」のリストから、「温度」「触媒」を選択します。(直前の選択が維持されている場合もあります。)

8. [追加]をクリックします。

9. 「OR」の下に表示されるパネルで、温度のスライダの上に表示されている式の「350」をクリックし、「650」に変更します。Enterキーを押します(または、どこか別の場所をクリックします)。

これにより、「温度」が650度以上の値になるような因子の設定は許可されなくなります。

10. 「触媒」の下の「C」のブロックをクリックします。

「温度」の制約と合わせて、「触媒」がCで「温度」が650度以上になるような因子の設定が許可されなくなります。

図5.47 制約の定義 

Image shown here

応答曲面の項の追加と計画の作成

1. 「モデル」アウトラインで[RSM]を選択します。

JMPの警告画面が開き、連続尺度の因子から構成される2次項のみをモデルに追加するというメッセージが表示されます。

2. [OK]をクリックして、警告を閉じます。

JMPは、適切な応答曲面の項だけをモデルに追加します。

注: 乱数シード値(ステップ3)と開始点の数(ステップ4)を設定すると、以下の数値例と同じ結果が得られます。同じ実験設定でなくても良い場合は、これらの手順は不要です。

3. (オプション)「カスタム計画」の赤い三角ボタンをクリックし、[乱数シード値の設定]を選択します。次に、「654321」と入力して[OK]をクリックします。

4. (オプション)「カスタム計画」の赤い三角ボタンをクリックし、[開始点の数]を選択して「1000」に設定します。[OK]をクリックします。

5. [計画の作成]をクリックします。

図5.48 制約を満たす計画 

Image shown here

「モデル」アウトラインでRSM(応答曲面の項)を追加したので、これはI-最適計画となります。この計画は、以下の制約を満たします。

触媒がBのときは、温度は400度を超えるものとする。

触媒がCのときは、温度は650度未満であるものとする。

計画の表示

はじめに、この計画の範囲として、「時間」は500~560秒の間、「温度」は350~750度の間と設定しました。このような制約をもつ計画の計画点をグラフ上にプロットしてみましょう。そのためには、次の操作を行います。

1. [テーブルの作成]をクリックします。

2. [グラフ]>[グラフビルダー]を選択します。

3. 「時間」「温度」を選択し、グラフの中央にドラッグします。

4. 平滑線のアイコンの選択を解除します。このアイコンは、グラフの上に並んでいるアイコンのうち左から2番目にあります。

5. 「触媒」を選択し、グラフ上側の「グループX」ゾーンにドラッグします。

図5.49 3水準の触媒の計画点 

Image shown here

6. 「温度」の軸をダブルクリックします。

7. 「参照線」で、次の操作を行います。

「値」に「350」と入力します。

[追加]をクリックします。

「値」に「750」と入力します。

[追加]をクリックします。

8. [OK]をクリックします。

図5.50 3水準の触媒の温度の範囲 

Image shown here

9. 「時間」の軸をダブルクリックします。

10. 「参照線」で、次の操作を行います。

「値」に「500」と入力します。

[追加]をクリックします。

「値」に「560」と入力します。

[追加]をクリックします。

11. [OK]をクリックします。

12. [終了]をクリックします。

注: 図5.51のプロットは、縦方向のサイズを変更したものです。

図5.51 3水準の触媒の計画の範囲 

Image shown here

触媒Aにおいては、6つの点が、計画の範囲の境界線上に位置しています。触媒Bにおいては、設定した制約によって、すべての点の温度が400度以上になっています。触媒Cにおいては、設定した制約によって、すべての点の温度が650度以下になっています。

より詳細な情報が必要な場合や、質問があるときは、JMPユーザーコミュニティで答えを見つけましょう (community.jmp.com).