704722_ja-JP

liliconcom

ブロックされた複製を使用した 2^3 階乗

こんにちは

 

1 回複製される 2^3 要因計画 (合計実行 = 16) を構築しようとしていました。 最後に複製をブロックする必要があるため、各ブロックには元の 2^3 要因計画が含まれている必要があります。 簡単な方法はありますか？ カスタム デザイン ツールを使用した私の方法は、各ブロックで半分階乗が複製されるという結果になりました。

 

ありがとう ：）

 

 

 

 

 

 

実験計画法Re: ブロックされた複製を使用した 2^3 階乗

こんにちは@Luca_Indrizzi 、

 

提供されたスクリーンショットからは、すべてが正常に期待どおりの結果に達しているように見えます。
ブロッキング係数はどのように設定しましたか? ブロックごとに 8 回の実行として設定しましたか (各ブロックで 2^3 要因計画を使用するため)。 :

 

「はい」の場合、モデルに 16 の実行と主効果を指定すると、期待どおりの結果が得られるはずです。

同じデザインを生成する場合のスクリプトは次のとおりです。

DOE(
 Custom Design,
 {Add Response( Maximize, "Y", ., ., . ),
 Add Factor( Continuous, -1, 1, "X1", 0 ),
 Add Factor( Continuous, -1, 1, "X2", 0 ),
 Add Factor( Continuous, -1, 1, "X3", 0 ), Add Factor( Blocking, 8, "X4" ),
 Set Random Seed( 1874557131 ), Number of Starts( 1 ), Add Term( {1, 0} ),
 Add Term( {1, 1} ), Add Term( {2, 1} ), Add Term( {3, 1} ), Add Term( {4, 1} ),
 Add Alias Term( {1, 1}, {2, 1} ), Add Alias Term( {1, 1}, {3, 1} ),
 Add Alias Term( {2, 1}, {3, 1} ), Set Sample Size( 16 ), Simulate Responses( 0 ),
 Save X Matrix( 0 ), Make Design}
)

見ていただくためにデータテーブル (色付き!) を添付しました。

 

完全な要因計画が必要な場合は、モデル内で 2 因子交互作用を指定する必要があります。同じ手順で、期待どおりに 2 つのブロックに 1 回完全な要因計画が複製されることになります。

これを生成するスクリプトは次のとおりです。

DOE(
 Custom Design,
 {Add Response( Maximize, "Y", ., ., . ),
 Add Factor( Continuous, -1, 1, "X1", 0 ),
 Add Factor( Continuous, -1, 1, "X2", 0 ),
 Add Factor( Continuous, -1, 1, "X3", 0 ), Add Factor( Blocking, 8, "X4" ),
 Set Random Seed( 1421844496 ), Number of Starts( 2 ), Add Term( {1, 0} ),
 Add Term( {1, 1} ), Add Term( {2, 1} ), Add Term( {3, 1} ), Add Term( {4, 1} ),
 Add Term( {1, 1}, {2, 1} ), Add Term( {1, 1}, {3, 1} ),
 Add Term( {2, 1}, {3, 1} ), Set Sample Size( 16 ),
 Optimality Criterion( "Make A-Optimal Design" ),
 "A-Optimality Parameter Weights"n( [1 1 1 1 1 1 1 1] ), Simulate Responses( 0 ),
 Save X Matrix( 0 ), Make Design, Set Run Order( Randomize within Blocks ),
 Make Table}
)

データテーブルも添付されています。

 

別のオプションは、8 実行要因計画を作成し、プラットフォームを使用することです。オーグメントデザイン" 実行を 1 回複製します。 ただし、最終的なデータテーブルにデザインのブロック/部分を示す列を導入し、それに応じて (手動で) 列のデータ型 (文字)、モデリング タイプ (公称)、および必要なプロパティを設定する必要があります。ブロック 1 と 2 の場合 ブロック 2 の場合)、RunsperBlock (8)、設計の役割 (ブロック)、および因子の変更 (簡単)

 

これがお役に立てば幸いです。

Re: ブロックされた複製を使用した 2^3 階乗

こんにちは@Victor_G

ご返信ありがとうございます。

正確には、1ブロックあたり8ランを設定しました。

しかし、あなたと私が得たのは、各ブロックに複製された 2^3-1 設計です。 1,1,1 はブロック 1 と 2 に 2 回現れます。 -1、-1、-1 はまったく表示されません。 私が望んでいるのは、ブロック 1 と 2 にそれぞれ 2^3 要因計画全体が含まれていることです。

よろしくお願いします

 

Re: ブロックされた複製を使用した 2^3 階乗

はい、最初の応答の後、あなたが半分の階乗ではなく完全な階乗を望んでいることがわかりましたので、修正しました:)

完全な要因計画を取得するには、モデルで 2 因子交互作用を指定する必要があります。

改訂された回答では、両方の設計のオプションがあり、2 つのデータテーブルとスクリプトが提供されています (1 つは半分完全階乗複製用、もう 1 つは完全階乗複製用)。

これがお役に立てば幸いです。

Re: ブロックされた複製を使用した 2^3 階乗

なるほど、Blocking Factor X4 の導入によりレプリケートを追加する必要がなくなりました、ありがとう@Victor_G

Re: ブロックされた複製を使用した 2^3 階乗

カスタム デザインの結果は、指定したモデルによって異なります。 この場合、すべての 2 方向および 3 方向のインタラクションを指定しましたか。 私の結果はあなたが求めたものです。

DOE(
 Custom Design,
 {Add Response( Maximize, "Y", ., ., . ),
 Add Factor( Continuous, -1, 1, "X1", 0 ),
 Add Factor( Continuous, -1, 1, "X2", 0 ),
 Add Factor( Continuous, -1, 1, "X3", 0 ), Add Factor( Blocking, 8, "X4" ),
 Set Random Seed( 21226362 ), Number of Starts( 9 ), Add Term( {1, 0} ),
 Add Term( {1, 1} ), Add Term( {2, 1} ), Add Term( {3, 1} ), Add Term( {4, 1} ),
 Add Term( {1, 1}, {2, 1} ), Add Term( {1, 1}, {3, 1} ),
 Add Term( {2, 1}, {3, 1} ), Add Term( {1, 1}, {2, 1}, {3, 1} ),
 Set Sample Size( 16 ), Simulate Responses( 0 ), Save X Matrix( 0 ), Make Design,
 Set Run Order( Randomize within Blocks ), Make Table}
);