@マークベイリーと@フィル_ケイ、最終的な答えは何ですか？ まったく同じ DOE テーブルを再現できますか? 理解していない場合は申し訳ありません。
私は同様の問題を抱えています。 I, 1) 3 つの混合係数、1 つの中心点、合計 24 点、I-optimal を使用してカスタム計画を作成しました。 2）「テーブルの作成」をクリックしました。 3) テーブル ウィンドウで、[DOE ダイアログ] をクリックしました。それから私は別のデザインを得た. ランダムシードを確認したところ、初期デザインと同じでした。
DOE(
Custom Design,
{Add Response( Maximize, "Y", ., ., . ), Add Factor( Mixture, 0.25, 0.7, "A", 0 ),
Add Factor( Mixture, 0.25, 0.7, "B", 0 ), Add Factor( Mixture, 0.01, 0.04, "C", 0 ),
Set Random Seed( 1096779613 ), Number of Starts( 6896 ), Add Term( {1, 1} ),
Add Term( {2, 1} ), Add Term( {3, 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} ),
Add Term( {1, 1}, {2, 1} ), Add Term( {1, 1}, {3, 1} ), Add Term( {2, 1}, {3, 1} ),
Center Points( 1 ), Set Sample Size( 24 ),
Optimality Criterion( Name( "Make I-Optimal Design" ) ), Simulate Responses( 0 ),
Save X Matrix( 0 ), Make Design, Set Run Order( Randomize ), Make Table}
);
ランダム シード = 1096779613
開始回数 = 6896
必要に応じてテーブル スクリプト:
New Table( "Custom Design",
Add Rows( 24 ),
New Table Variable( "Design", "Custom Design" ),
New Table Variable( "Criterion", "I Optimal" ),
New Script(
"Model",
Fit Model(
Effects(
:A & Mixture,
:B & Mixture,
:C & Mixture,
:A * :B,
:A * :C,
:B * :C,
:A * :B * :C,
Scheffe Cubic( A, B ),
Scheffe Cubic( A, C ),
Scheffe Cubic( B, C )
),
Y( :Y ),
No Intercept( 1 )
)
),
New Script( "Evaluate Design", DOE( Evaluate Design, X( :A, :B, :C ) ) ),
New Script(
"DOE Dialog",
DOE(
Custom Design,
{Add Response( Maximize, "Y", ., ., . ), Add Factor( Mixture, 0.25, 0.7, "A", 0 ),
Add Factor( Mixture, 0.25, 0.7, "B", 0 ), Add Factor( Mixture, 0.01, 0.04, "C", 0 ),
Set Random Seed( 1096779613 ), Number of Starts( 6896 ), Add Term( {1, 1} ),
Add Term( {2, 1} ), Add Term( {3, 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} ), Add Term( {1, 1}, {2, 1} ),
Add Term( {1, 1}, {3, 1} ), Add Term( {2, 1}, {3, 1} ), Center Points( 1 ),
Set Sample Size( 24 ), Optimality Criterion( Name( "Make I-Optimal Design" ) ),
Simulate Responses( 0 ), Save X Matrix( 0 ), Make Design}
)
),
New Column( "A",
Numeric,
"Continuous",
Format( "Best", 12 ),
Set Property( "Mixture", {0.25, 0.7, 1, L PseudoComponent Coding} ),
Set Property( "Design Role", DesignRole( Mixture ) ),
Set Property( "Factor Changes", Easy ),
Set Values(
[0.312099999999998, 0.570700000000011, 0.3504, 0.437800000000002, 0.372, 0.7, 0.3504,
0.487411764705882, 0.615000000000001, 0.525199999999999, 0.535999999999999,
0.437799999999998, 0.7, 0.348, 0.3642, 0.612, 0.447400000000013, 0.330999999999991,
0.6563, 0.676, 0.5716, 0.331000000000001, 0.612599999999999, 0.699999999999997]
)
),
New Column( "B",
Numeric,
"Continuous",
Format( "Best", 12 ),
Set Property( "Mixture", {0.25, 0.7, 1, L PseudoComponent Coding} ),
Set Property( "Design Role", DesignRole( Mixture ) ),
Set Property( "Factor Changes", Easy ),
Set Values(
[0.656900000000003, 0.3983, 0.6126, 0.525200000000002, 0.618000000000005, 0.2666,
0.6126, 0.487411764705882, 0.36, 0.437800000000001, 0.454000000000005,
0.525200000000002, 0.29, 0.612, 0.6198, 0.348, 0.533599999999996, 0.659000000000019,
0.3067, 0.29, 0.4004, 0.659, 0.350400000000001, 0.26687487586833]
)
),
New Column( "C",
Numeric,
"Continuous",
Format( "Best", 12 ),
Set Property( "Mixture", {0.01, 0.04, 1, L PseudoComponent Coding} ),
Set Property( "Design Role", DesignRole( Mixture ) ),
Set Property( "Factor Changes", Easy ),
Set Values(
[0.031, 0.031, 0.037, 0.037, 0.01, 0.0333999999999994, 0.0370000000000005,
0.0251764705882343, 0.025, 0.0370000000000026, 0.01, 0.0370000000000022, 0.01, 0.04,
0.0159999999999976, 0.04, 0.019, 0.01, 0.0370000000000026, 0.034, 0.028, 0.01,
0.0370000000000035, 0.0331251241316703]
)
),
New Column( "Y",
Numeric,
"Continuous",
Format( "Best", 12 ),
Set Property(
"Response Limits",
{Goal( Maximize ), Lower( . ), Upper( . ), Importance( . )}
),
Set Values( [., ., ., ., ., ., ., ., ., ., ., ., ., ., ., ., ., ., ., ., ., ., ., .] )
)
)
