無法建立具有大量約束的 DOE<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />我的客戶有以下情況來設置 DOE:
混合物的 3 個成分(S、X、M,可以視為混合物或連續)加上填料成分 (W),以及成分 S 和 X 各自的兩個供應商,稱為 OS 和 OX(每個分類有兩個水平)。
限制帶來的困難,首先是范圍:
範圍 S 0 -0.04
範圍 X 0-0.006
範圍 M 0-0.07
範圍 W 0.884-1
操作系統 L1、L2
牛 L1, L2
或者是連續的,其中範圍表示百分比並且不再需要填料,因為他們對填料本身的影響根本不感興趣
範圍 S 0 -4 [%]
範圍 X 0-0.6 [%]
範圍 M 0-7 [%]
操作系統 L1、L2
牛 L1, L2
此外還有以下要求:
- S+X >0,所以S和X不能同時為0
- 如果 X==0,S 必須 >0.5%(但如果 X!=0,S 可以為 0)
- 如果 S ==0,X 必須 >0.05%(但如果 S!=0,X 可以為 0)
- 如果 S==0,則不允許操作系統/沒有意義
- 如果 X==0,OX 是不允許的/沒有意義
問題是除了範圍或簡單的線性約束之外,不能有混合因子的約束來確保 S+X>0。 在連續情況下,我們無法設法以滿足所有約束的方式設置約束。 不是線性組合,也不是不允許的組合,兩者的混合也是不可能的。
有什麼想法可以使用上述設置創建設計嗎?
附上我到目前為止嘗試過的一些東西。 在一個中,我創建了一個接近要求的設計,但認為它不夠平衡,因為沒有單獨出現 S 或 X 的運行。 請參閱三個腳本來說明行為,尤其是具有 3 級分類的腳本,因為其中包含多種設計。 為方便起見,因子保存在表格中
感謝您的支持,
馬丁
回复:不能設置有很多約束的 DOE好的,非常感謝瑞恩!
回复:不能設置有很多約束的 DOE嗨,馬丁,
沒有什麼能引起我的注意 - 僅僅因為空間有限,看到不允許的組合看起來很糟糕的診斷並不罕見。
至於內在層次的強制,如果你添加一些必要的權力,這是一種選擇,或者你可以嘗試盡可能保留至少一些權力。
如果您使用更簡單的模型(具有必要的權力)設計作為比較設計中的參考設計,您可以獲得更直接的比較,以了解可能首選的設計。
乾杯,
瑞安
回复:不能設置有很多約束的 DOE嗨瑞安,
可能你不需要回答我上面的回复。
我與客戶一起嘗試了更多的東西,並通過移除效果 S*OS、X*OX 和 OS*OX 提出了一個非 splot-plot 設計:
美國能源部(
定制設計,
{添加響應(最大化, “是” , . , . , . ),
添加因子(離散數字,{ 0 , 1 , 2 , 3 , 4 }, “S” , 0 ),
添加因子(離散數字,{ 0 , 0.15 , 0.3 , 0.45 , 0.6 }, “X” , 0 ),
添加因子(離散數字,{ 0 , 1.75 , 3.5 , 5.25 , 7 }, “米” , 0 ),
添加因子(分類,{ “L1” , “L2” , “沒有” }, “操作系統” , 0 ),
添加因子(分類,{ “L1” , “L2” , “沒有” }, “牛” , 0 ),
設置隨機種子( 695079790 ), 開始次數( 1000 ), 添加項( { 1 , 0 }),
添加項({ 1 , 1 } ), 添加項( { 2 , 1 } ), 添加項( { 3 , 1 } ), 添加項( { 4 , 1 }),
添加項({ 5 , 1 } ), 添加項( { 1 , 1 }, { 2 , 1 } ), 添加項( { 1 , 1 }, { 3 , 1 }),
添加項({ 1 , 1 }, { 5 , 1 } ), 添加項( { 2 , 1 }, { 3 , 1 }),
添加項({ 2 , 1 }, { 4 , 1 } ), 添加項( { 3 , 1 }, { 4 , 1 }),
添加項({ 3 , 1 }, { 5 , 1 } ), 設置樣本大小( 36 ),
不允許的組合(
小號 == 0 & X == 0 | 0 < S <= 0.5 | S > 0 & 操作系統 == 3 |小號 == 0 & (操作系統 == 1 |操作系統
== 2 ) | X > 0 & 牛== 3 | X == 0 & (OX == 1 |牛== 2 )
)}
);
這似乎是相當好的預測方差,FDS 圖看起來也不錯。對於大多數效果,功率都非常好,相關性也很好。 當然,S*OS 和 X*OX 的別名效應更大。
你會選擇那個設計,還是有什麼你會說讓你感到緊張的事情?
最後一個問題:有沒有辦法強制在設計中也有 S 和 X 的內部級別?有時濃度的下半部分可能遵循與較高濃度不同的模式
謝謝,
馬丁
回复:不能設置有很多約束的 DOE嗨瑞安,
請先閱讀第二個答复,因為該答复可能會使此答復過時。
所以這是我的最終腳本,如果添加更多運行(例如 36),它具有不錯的功能和預測方差:
// 24 次運行
能源部(
定制設計,
{添加響應(最大化,“Y”,.,.,.),
添加因子(離散數值,{0, 1, 2, 3, 4}, "S", 0 ),
添加因子(離散數值,{0, 0.15, 0.3, 0.45, 0.6}, "X", 0),
添加因子(離散數值,{0, 1.75, 3.5, 5.25, 7}, "M", 0),
添加因子(分類,{“L1”,“L2”,“無”},“OS”,0),
添加因子(分類,{“L1”,“L2”},“OX”,1),設置隨機種子(984492505),
開始次數(1000),添加項({1, 0}),添加項({1, 1}),
添加項({2, 1}),添加項({3, 1}),添加項({4, 1}),添加項({5, 1}),
添加項({1, 1}, {2, 1}),添加項({1, 1}, {3, 1}),
添加項({1, 1}, {5, 1}),添加項({2, 1}, {3, 1}),
添加項({2, 1}, {4, 1}),添加項({2, 1}, {5, 1}),
添加項({3, 1}, {4, 1}),添加項({3, 1}, {5, 1}),
添加項({4, 1}, {5, 1}),優化整區數,
設置樣本量(24),不允許的組合(
0 < S <= 0.5 | S > 0 & 操作系統 == 3 | S == 0 & (OS == 1 | OS == 2) | X < 0.05
)}
);
// 36 次運行,具有更大的功率和更少的預測方差
能源部(
定制設計,
{添加響應(最大化,“Y”,.,.,.),
添加因子(離散數值,{0, 1, 2, 3, 4}, "S", 0 ),
添加因子(離散數值,{0, 0.15, 0.3, 0.45, 0.6}, "X", 0),
添加因子(離散數值,{0, 1.75, 3.5, 5.25, 7}, "M", 0),
添加因子(分類,{“L1”,“L2”,“無”},“OS”,0),
添加因子(分類,{“L1”,“L2”},“OX”,1),設置隨機種子(1173606031),
開始次數(1000),添加項({1, 0}),添加項({1, 1}),
添加項({2, 1}),添加項({3, 1}),添加項({4, 1}),添加項({5, 1}),
添加項({1, 1}, {2, 1}),添加項({1, 1}, {3, 1}),
添加項({1, 1}, {5, 1}),添加項({2, 1}, {3, 1}),
添加項({2, 1}, {4, 1}),添加項({2, 1}, {5, 1}),
添加項({3, 1}, {4, 1}),添加項({3, 1}, {5, 1}),
添加項({4, 1}, {5, 1}),添加別名項({1, 1}, {4, 1}),
設置 N 個整區(4),設置樣本量(36),
不允許的組合(
0 < S <= 0.5 | S > 0 & 操作系統 == 3 | S == 0 & (OS == 1 | OS == 2) | X < 0.05
)}
);
正如您所建議的,我手動刪除了高階項以及 S*OS 交互。 感謝那。
我有兩個附加問題:
- 1)我剛收到回复。設置分裂圖後,我們必須始終讓 X 進入。 但他們也希望獲得 X == 0 但 S > 0.5 的估計值。 我試圖用 OX 建立一個新的設計,它有 3 個級別。 我在 X == 0 時強制 OX "None" 遇到麻煩,如果 X > 0 強制 OX != "None" (3) 並且仍然得到功率計算,所以我刪除了交互效果 X*OX 和 OS *牛。 這是我最終想出的設計:
美國能源部(
定制設計,
{添加響應(最大化, “是” , . , . , . ),
添加因子(離散數字,{ 0 , 1 , 2 , 3 , 4 }, “S” , 0 ),
添加因子(離散數字,{ 0 , 0.15 , 0.3 , 0.45 , 0.6 }, “X” , 0 ),
添加因子(離散數字,{ 0 , 1.75 , 3.5 , 5.25 , 7 }, “米” , 0 ),
添加因子(分類,{ “L1” , “L2” , “沒有” }, “操作系統” , 0 ),
添加因子(分類,{ “L1” , “L2” , “沒有” }, “牛” , 1 ),
設置隨機種子( 1317349002 ), 開始次數( 1000 ), 添加項( { 1 , 0 }),
添加項({ 1 , 1 } ), 添加項( { 2 , 1 } ), 添加項( { 3 , 1 } ), 添加項( { 4 , 1 }),
添加項({ 5 , 1 } ), 添加項( { 1 , 1 }, { 2 , 1 } ), 添加項( { 1 , 1 }, { 3 , 1 }),
添加項({ 1 , 1 }, { 5 , 1 } ), 添加項( { 2 , 1 }, { 3 , 1 }),
添加項({ 2 , 1 }, { 4 , 1 } ), 添加項( { 3 , 1 }, { 4 , 1 }),
添加項({ 3 , 1 }, { 5 , 1 } ), 添加別名項( { 1 , 1 }, { 4 , 1 }),
設置 N 個整區( 6 ), 設置樣本量( 36 ),
不允許的組合(
0 < S <= 0.5 | S > 0 & 操作系統 == 3 |小號 == 0 & (操作系統 == 1 |操作系統 == 2 ) | X > 0 &牛
== 3 | X == 0 & (OX == 1 |牛== 2 )
)}
);
我應該避免使用這種設計嗎? 除了設計空間情節的一部分似乎都還可以。 我認為即使是 FDS 的情節也不錯。 你會改變什麼。
2) 保留最初的裂區設計並為 S/OS 而不是 X/OX 做另一個以獲取 X==0 選項的估計值是否有意義? 如果整體上不同組合的供應商存在顯著差異,我認為將兩個結果結合起來回答這個問題不是一個好主意,你同意嗎?
回复:不能設置有很多約束的 DOE哎呀,我之前看的不夠仔細…… 我認為這裡發生了幾件事 - 首先,如果可能,將添加更高階的離散數字項(這些必須手動刪除)。 更大的問題是 S*OS 交互項,因為它不能用不允許的組合來估計。 設計師應該已經抓住了這一點,但沒有,我將不得不進一步調查。 運行以下內容並刪除那些可能的術語,事情應該看起來更好。 我會看看我們是否可以在未來的版本中解決這些問題。
您可以嘗試刪除該交互並返回到連續路線......但您可能想要 S < 0.01 之類的東西,而不是使用 S==0。 (注意,由於存儲雙精度數,這也可能是離散數字的問題,所以即使那樣我也需要使用 abs(S)<.0001 而不是 S==0)。
乾杯,
瑞安
美國能源部(定制設計,
{添加響應(最大化, “是” , . , . , . ),
添加因子(離散數字,{ 0 , 1 , 2 , 3 , 4 }, “S” , 0 ),
添加因子(離散數字,{ 0 , 0.15 , 0.3 , 0.45 , 0.6 }, “X” , 0 ),
添加因子(離散數字,{ 0 , 1.75 , 3.5 , 5.25 , 7 }, “米” , 0 ),
添加因子(分類,{ “L1” , “L2” , “沒有” }, “操作系統” , 0 ),
添加因子(分類,{ “L1” , “L2” }, “牛” , 1 ), 設置隨機種子( 891554152 ),
開始次數( 1000 ), 添加項( { 1 , 0 } ), 添加項( { 1 , 1 }),
添加項({ 2 , 1 } ), 添加項( { 3 , 1 } ), 添加項( { 4 , 1 } ), 添加項( { 5 , 1 }),
添加項({ 1 , 1 }, { 2 , 1 } ), 添加項( { 1 , 1 }, { 3 , 1 }),
添加項({ 1 , 1 }, { 5 , 1 } ), 添加項( { 2 , 1 }, { 3 , 1 }),
添加項({ 2 , 1 }, { 4 , 1 } ), 添加項( { 2 , 1 }, { 5 , 1 }),
添加項({ 3 , 1 }, { 4 , 1 } ), 添加項( { 3 , 1 }, { 5 , 1 }),
添加項({ 4 , 1 }, { 5 , 1 } ), 添加別名項( { 1 , 1 }, { 4 , 1 }),
設置 N 個整區( 4 ), 設置樣本量( 28 ),
0 < S <= 0.5 | S > 0 & 操作系統 == 3 |小號 == 0 & (操作系統 == 1 |操作系統 == 2 ))});
回复:不能設置有很多約束的 DOE嗨瑞安,
不,我沒有。 它似乎有效,我將不得不與客戶核實這是否足夠好並且可以使用離散值。
我看到設計評估和功率分析的一件事,我無法估計功率,但我想這是由於離散數字,對吧? 顏色圖對我來說看起來不錯:
設計空間圖的分數在 0,6 之後看起來很奇怪:

那是因為約束集嗎?
謝謝,
馬丁
回复:不能設置有很多約束的 DOE我認為您在使用不允許的連續因子等式組合時遇到了問題。 特別是,當使用 X==0 或 S==0 時,最終可能會使用真正接近 0(甚至顯示為 0)的東西。 您是否嘗試過將離散數字用於連續? 使用您最後的腳本之一,如下所示:
美國能源部(
{添加響應(最大化, “是” , . , . , . ),
//添加因子(Continuous, 0, 4, "S", 0 ),
//添加因子( Continuous, 0, 0.6, "X", 0 ), 添加因子( Continuous, 0, 7, "M", 0 ),
添加因子(離散數字,{ 0 , 1 , 2 , 3 , 4 }, “S” , 0 ),
添加因子(離散數字,{ 0 , .15 , .3 , .45 , .6 }, “X” , 0 ),
添加因子(離散數字,{ 0 , 1.75 , 3.5 , 5.25 , 7 }, “米” , 0 ),
添加因子(分類,{ “L1” , “L2” , “沒有” }, “操作系統” , 0 ),
添加因子(分類,{ “L1” , “L2” }, “牛” , 1 ), 設置隨機種子( 513965188 ),
開始次數( 1000 ), 添加項( { 1 , 0 } ), 添加項( { 1 , 1 }),
添加項({ 2 , 1 } ), 添加項( { 3 , 1 } ), 添加項( { 4 , 1 } ), 添加項( { 5 , 1 }),
添加項({ 1 , 1 }, { 2 , 1 } ), 添加項( { 1 , 1 }, { 3 , 1 }),
添加項({ 1 , 1 }, { 4 , 1 } ), 添加項( { 1 , 1 }, { 5 , 1 }),
添加項({ 2 , 1 }, { 3 , 1 } ), 添加項( { 2 , 1 }, { 4 , 1 }),
添加項({ 2 , 1 }, { 5 , 1 } ), 添加項( { 3 , 1 }, { 4 , 1 }),
添加項({ 3 , 1 }, { 5 , 1 } ), 添加項( { 4 , 1 }, { 5 , 1 } ), 設置 N 個整區 ( 4 ),
設置樣本量( 28 ), 不允許的組合(
X < 0.05 | 0 < S <= 0.5 | S > 0 & 操作系統 == 3 |小號 == 0 & (操作系統 == 1 |操作系統 == 2 )
Re: 不能設置有很多約束的 DOE嗨,馬丁,
您是否嘗試過使用不允許的組合腳本? 這為您定義約束提供了更大的靈活性。
http://www.jmp.com/support/help/Define_Factor_Constraints.shtml#205355
問候,
菲爾