DOE mit vielen Einschränkungen kann nicht eingerichtet werden<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Mein Kunde hat die folgende Situation, um eine DOE einzurichten:
3 Komponenten einer Mischung (S, X, M, kann als Mischung oder kontinuierlich behandelt werden) plus eine Füllstoffkomponente (W), sowie zwei Lieferanten für jede der Zutaten S und X, OS und OX genannt (jeweils kategorial mit zwei Ebenen).
Die Schwierigkeiten steigen bei den Einschränkungen, zuerst die Bereiche:
Bereich S 0 -0,04
Bereich X 0-0,006
Bereich M 0-0,07
Bereich W 0,884-1
Betriebssystem L1, L2
Ochse L1, L2
Oder als kontinuierlich, wo die Reichweite den Prozentsatz ausdrückt und der Füllstoff nicht mehr benötigt wird, da sie sich überhaupt nicht für die Wirkung des Füllstoffs selbst interessieren
Bereich S 0 -4 [%]
Bereich X 0-0,6 [%]
Bereich M 0-7 [%]
Betriebssystem L1, L2
Ochse L1, L2
Darüber hinaus gibt es folgende Anforderungen:
- S+X >0, sodass S und X nicht gleichzeitig 0 sein können
- Wenn X==0, muss S >0,5 % sein (aber wenn X!=0, S kann 0 sein)
- Wenn S ==0, muss X >0,05 % sein (aber wenn S!=0, kann X 0 sein)
- Wenn S==0, ist OS nicht erlaubt/sinnlos
- Bei X==0 ist OX nicht erlaubt/sinnlos
Das Problem ist, dass es keine anderen Beschränkungen für die Mischungsfaktoren als die Bereiche oder eine einfache lineare Beschränkung geben kann, um sicherzustellen, dass S+X > 0 ist. Und im kontinuierlichen Fall konnten wir es nicht schaffen, die Bedingungen so zu setzen, dass alle Bedingungen erfüllt sind. Weder bei linearen, noch bei unzulässigen Kombinationen, und eine Mischung aus beidem ist auch nicht möglich.
Irgendwelche Ideen, wie ich ein Design mit den obigen Einstellungen erstellen könnte?
Anbei einige Dinge, die ich bisher versucht habe. In einem habe ich ein Design erstellt, das den Anforderungen nahe kommt, aber ich fand es nicht ausgewogen genug, da es keinen Lauf gibt, in dem S oder X allein vorkommen. Bitte sehen Sie sich die drei Skripte an, um die Verhaltensweisen zu veranschaulichen, insbesondere das mit der 3-stufigen Kategorie, da mehrere Designs enthalten sind. Faktoren werden bequem in Tabellen gespeichert
Danke für deine Unterstützung,
Martin
Betreff: DOE mit vielen Einschränkungen kann nicht eingerichtet werdenOK, vielen Dank Ryan!
Betreff: DOE mit vielen Einschränkungen kann nicht eingerichtet werdenHallo Martin,
Mir fällt nichts auf - und es ist nicht ungewöhnlich, Diagnosen zu sehen, die bei unzulässigen Kombinationen schlecht erscheinen können, einfach wegen des begrenzten Platzes.
Was das Erzwingen der inneren Ebenen betrifft, so ist das eine Option, wenn Sie einige der Kräfte nach Bedarf hinzufügen, oder Sie könnten versuchen, zumindest einige der Kräfte nach Möglichkeit beizubehalten.
Wenn Sie das einfachere Modelldesign (mit erforderlichen Potenzen) als Referenzdesign in Designs vergleichen verwenden, können Sie einen direkteren Vergleich erhalten, um festzustellen, welches Design bevorzugt werden könnte.
Prost,
Ryan
Betreff: DOE mit vielen Einschränkungen kann nicht eingerichtet werdenHallo Ryan,
Vielleicht brauchen Sie meine obige Antwort nicht zu beantworten.
Ich habe noch ein paar Dinge mit dem Kunden ausprobiert und bin zu einem Nicht-Splot-Plot-Design gekommen, indem ich einfach die Effekte S*OS, X*OX und OS*OX entfernt habe:
DAMHIRSCHKUH (
Maßgefertigtes Design,
{Antwort hinzufügen( maximieren, "Y" , . , . , . ),
Faktor hinzufügen( diskret numerisch, { 0 , 1 , 2 , 3 , 4 }, "S" , 0 ),
Faktor hinzufügen( diskret numerisch, { 0 , 0,15 , 0,3 , 0,45 , 0,6 }, "X" , 0 ),
Faktor hinzufügen( diskret numerisch, { 0 , 1,75 , 3.5 , 5.25 , 7 }, "M" , 0 ),
Faktor hinzufügen( Kategorial, { "L1" , "L2" , "Keiner" }, "Betriebssystem" , 0 ),
Faktor hinzufügen( Kategorial, { "L1" , "L2" , "Keiner" }, "OCHSE" , 0 ),
Zufallsstartwert festlegen ( 695079790 ), Anzahl der Starts ( 1000 ), Begriff hinzufügen ( { 1 , 0 } ),
Begriff hinzufügen ( { 1 , 1 } ), Begriff hinzufügen ( { 2 , 1 } ), Begriff hinzufügen ( { 3 , 1 } ), Begriff hinzufügen ( { 4 , 1 } ),
Begriff hinzufügen ( { 5 , 1 } ), Begriff hinzufügen ( { 1 , 1 }, { 2 , 1 } ), Begriff hinzufügen ( { 1 , 1 }, { 3 , 1 } ),
Begriff hinzufügen ( { 1 , 1 }, { 5 , 1 } ), Begriff hinzufügen ( { 2 , 1 }, { 3 , 1 } ),
Begriff hinzufügen ( { 2 , 1 }, { 4 , 1 } ), Begriff hinzufügen ( { 3 , 1 }, { 4 , 1 } ),
Begriff hinzufügen ( { 3 , 1 }, { 5 , 1 } ), Stichprobengröße festlegen ( 36 ),
Unzulässige Kombinationen(
s == 0 & X == 0 | 0 < S <= 0,5 | S > 0 & Betriebssystem == 3 | S == 0 & (Betriebssystem == 1 | Betriebssystem
== 2 ) | X > 0 & OX == 3 | X == 0 & (OX == 1 | Ochse == 2 )
)}
);
Dies scheint eine ziemlich gute Vorhersagevarianz zu sein, und der FDS-Plot sieht auch in Ordnung aus. Die Leistung ist für die meisten Effekte ziemlich gut und die Korrelation ist auch in Ordnung. Natürlich ist der Alias-Effekt bei S*OS und X*OX größer.
Würden Sie sich für dieses Design entscheiden oder gibt es irgendetwas, was Sie nervös machen würde?
Und eine letzte Frage: Gibt es eine Möglichkeit, die inneren Ebenen für S und X auch im Design zu haben? manchmal kann die untere Hälfte der Konzentrationen einem anderen Muster folgen als die höheren
Danke,
Martin
Betreff: DOE mit vielen Einschränkungen kann nicht eingerichtet werdenHallo Ryan,
Bitte lesen Sie zuerst die zweite Antwort, da diese wahrscheinlich diese obsolet macht.
Hier ist also mein endgültiges Skript, das eine anständige Leistung und Vorhersagevarianz hat, wenn man weitere Läufe hinzufügt (z. B. 36):
// 24 Läufe
DAMHIRSCHKUH(
Maßgefertigtes Design,
{Antwort hinzufügen (Maximieren, "Y", ., ., . ),
Faktor hinzufügen (diskret numerisch, {0, 1, 2, 3, 4}, "S", 0 ),
Faktor hinzufügen (diskret numerisch, {0, 0,15, 0,3, 0,45, 0,6}, "X", 0 ),
Faktor hinzufügen (diskret numerisch, {0, 1,75, 3,5, 5,25, 7}, "M", 0 ),
Add Factor( Categorical, {"L1", "L2", "None"}, "OS", 0 ),
Faktor hinzufügen (Kategorial, {"L1", "L2"}, "OX", 1 ), Zufallsstartwert festlegen ( 984492505 ),
Anzahl der Starts ( 1000 ), Term hinzufügen ( {1, 0} ), Term hinzufügen ( {1, 1} ),
Begriff hinzufügen( {2, 1} ), Begriff hinzufügen( {3, 1} ), Begriff hinzufügen( {4, 1} ), Begriff hinzufügen( {5, 1} ),
Begriff hinzufügen( {1, 1}, {2, 1} ), Begriff hinzufügen( {1, 1}, {3, 1} ),
Begriff hinzufügen( {1, 1}, {5, 1} ), Begriff hinzufügen( {2, 1}, {3, 1} ),
Begriff hinzufügen( {2, 1}, {4, 1} ), Begriff hinzufügen( {2, 1}, {5, 1} ),
Begriff hinzufügen( {3, 1}, {4, 1} ), Begriff hinzufügen( {3, 1}, {5, 1} ),
Term hinzufügen ( {4, 1}, {5, 1} ), Anzahl ganzer Plots optimieren,
Stichprobenumfang festlegen( 24 ), Unzulässige Kombinationen(
0 < S <= 0,5 | S > 0 & OS == 3 | S == 0 & (OS == 1 | OS == 2) | X < 0,05
)}
);
// 36 Läufe, hat viel mehr Leistung und weniger Vorhersagevarianz
DAMHIRSCHKUH(
Maßgefertigtes Design,
{Antwort hinzufügen (Maximieren, "Y", ., ., . ),
Faktor hinzufügen (diskret numerisch, {0, 1, 2, 3, 4}, "S", 0 ),
Faktor hinzufügen (diskret numerisch, {0, 0,15, 0,3, 0,45, 0,6}, "X", 0 ),
Faktor hinzufügen (diskret numerisch, {0, 1,75, 3,5, 5,25, 7}, "M", 0 ),
Add Factor( Categorical, {"L1", "L2", "None"}, "OS", 0 ),
Faktor hinzufügen (kategorisch, {"L1", "L2"}, "OX", 1 ), Zufallsstartwert festlegen ( 1173606031 ),
Anzahl der Starts ( 1000 ), Term hinzufügen ( {1, 0} ), Term hinzufügen ( {1, 1} ),
Begriff hinzufügen( {2, 1} ), Begriff hinzufügen( {3, 1} ), Begriff hinzufügen( {4, 1} ), Begriff hinzufügen( {5, 1} ),
Begriff hinzufügen( {1, 1}, {2, 1} ), Begriff hinzufügen( {1, 1}, {3, 1} ),
Begriff hinzufügen( {1, 1}, {5, 1} ), Begriff hinzufügen( {2, 1}, {3, 1} ),
Begriff hinzufügen( {2, 1}, {4, 1} ), Begriff hinzufügen( {2, 1}, {5, 1} ),
Begriff hinzufügen( {3, 1}, {4, 1} ), Begriff hinzufügen( {3, 1}, {5, 1} ),
Begriff hinzufügen( {4, 1}, {5, 1} ), Alias-Begriff hinzufügen ( {1, 1}, {4, 1} ),
Setze N ganze Plots( 4 ), Setze Stichprobengröße( 36 ),
Unzulässige Kombinationen(
0 < S <= 0,5 | S > 0 & OS == 3 | S == 0 & (OS == 1 | OS == 2) | X < 0,05
)}
);
Wie Sie vorgeschlagen haben, habe ich die Terme höherer Ordnung sowie die S * OS-Interaktion manuell entfernt. Dank dafür.
Ich habe zwei Zusatzfragen:
- 1) Ich habe gerade eine Antwort bekommen. Mit dem so eingerichteten Split Plot sind wir gezwungen, X immer drin zu haben. Aber sie hätten gerne auch die Schätzung für X == 0, aber S >0,5. Ich habe versucht, ein neues Design mit OX mit 3 Ebenen einzurichten. Ich habe Probleme, OX "None" zu erzwingen, wenn X == 0 , und wenn X > 0 OX != "None" (3) erzwingt und trotzdem Leistungsberechnungen erhalte, habe ich den Interaktionseffekt X*OX und OS entfernt *OCHSE. Hier ist mein Design, auf das ich endlich gekommen bin:
DAMHIRSCHKUH (
Maßgefertigtes Design,
{Antwort hinzufügen( Maximieren, "Y" , . , . , . ),
Faktor hinzufügen( diskret numerisch, { 0 , 1 , 2 , 3 , 4 }, "S" , 0 ),
Faktor hinzufügen( diskret numerisch, { 0 , 0,15 , 0,3 , 0,45 , 0,6 }, "X" , 0 ),
Faktor hinzufügen( diskret numerisch, { 0 , 1,75 , 3.5 , 5.25 , 7 }, "M" , 0 ),
Faktor hinzufügen( Kategorial, { "L1" , "L2" , "Keiner" }, "Betriebssystem" , 0 ),
Faktor hinzufügen( Kategorial, { "L1" , "L2" , "Keiner" }, "OCHSE" , 1 ),
Zufallsstartwert festlegen ( 1317349002 ), Anzahl der Starts ( 1000 ), Begriff hinzufügen ( { 1 , 0 } ),
Begriff hinzufügen ( { 1 , 1 } ), Begriff hinzufügen ( { 2 , 1 } ), Begriff hinzufügen ( { 3 , 1 } ), Begriff hinzufügen ( { 4 , 1 } ),
Begriff hinzufügen ( { 5 , 1 } ), Begriff hinzufügen ( { 1 , 1 }, { 2 , 1 } ), Begriff hinzufügen ( { 1 , 1 }, { 3 , 1 } ),
Begriff hinzufügen ( { 1 , 1 }, { 5 , 1 } ), Begriff hinzufügen ( { 2 , 1 }, { 3 , 1 } ),
Begriff hinzufügen ( { 2 , 1 }, { 4 , 1 } ), Begriff hinzufügen ( { 3 , 1 }, { 4 , 1 } ),
Begriff hinzufügen ( { 3 , 1 }, { 5 , 1 } ), Alias-Begriff hinzufügen ( { 1 , 1 }, { 4 , 1 } ),
Setze N ganze Plots( 6 ), Stichprobengröße festlegen ( 36 ),
Unzulässige Kombinationen(
0 < S <= 0,5 | S > 0 & Betriebssystem == 3 | S == 0 & (Betriebssystem == 1 | Betriebssystem == 2 ) | X > 0 & OX
== 3 | X == 0 & (OX == 1 | Ochse == 2 )
)}
);
Gibt es irgendetwas, das ich von der Verwendung dieses Designs fernhalten sollte? Bis auf den Bruchteil des Designspace-Plots scheint alles in Ordnung zu sein. Und selbst der FDS-Plot ist gar nicht so schlecht, finde ich. Würdest du etwas ändern.
2) Wäre es sinnvoll, das anfängliche Split-Plot-Design beizubehalten und ein weiteres für S/OS anstelle von X/OX zu erstellen, um auch für die Option X==0 eine Schätzung zu erhalten? Ich halte es nicht für eine gute Idee, beide Ergebnisse zu kombinieren, um die Frage zu beantworten, ob es signifikante Unterschiede bei den Anbietern für die verschiedenen Kombinationen insgesamt gibt, stimmen Sie zu?
Betreff: DOE mit vielen Einschränkungen kann nicht eingerichtet werdenUps, ich hatte vorher nicht genau genug hingeschaut... Ich denke, hier passiert ein paar Dinge - erstens werden die diskreten numerischen Terme für höhere Ordnungen nach Möglichkeit hinzugefügt (diese müssen manuell entfernt werden). Das größere Problem ist der S*OS-Interaktionsterm, da er mit den unzulässigen Kombinationen nicht geschätzt werden kann. Der Designer hätte das bemerken sollen, hat es aber nicht, ich muss es mir genauer ansehen. Wenn Sie das Folgende ausführen und diese, wenn möglich, Begriffe entfernen, sollten die Dinge besser aussehen. Ich werde sehen, ob wir diese Probleme in einer zukünftigen Version beheben können.
Sie könnten versuchen, diese Interaktion zu entfernen und zur kontinuierlichen Route zurückzukehren ... aber anstatt S==0 zu verwenden, möchten Sie vielleicht etwas wie S < 0,01. (Hinweis: Dies kann aufgrund der Speicherung von Doubles auch ein Problem mit diskreten Zahlen sein, sodass ich selbst dann abs (S) <.0001 anstelle von S == 0 verwenden musste).
Prost,
Ryan
DAMHIRSCHKUH (Maßgefertigtes Design,
{Antwort hinzufügen( maximieren, "Y" , . , . , . ),
Faktor hinzufügen( diskret numerisch, { 0 , 1 , 2 , 3 , 4 }, "S" , 0 ),
Faktor hinzufügen( diskret numerisch, { 0 , 0,15 , 0,3 , 0,45 , 0,6 }, "X" , 0 ),
Faktor hinzufügen( diskret numerisch, { 0 , 1,75 , 3.5 , 5.25 , 7 }, "M" , 0 ),
Faktor hinzufügen( Kategorial, { "L1" , "L2" , "Keiner" }, "Betriebssystem" , 0 ),
Faktor hinzufügen( Kategorial, { "L1" , "L2" }, "OCHSE" , 1 ), Zufallsstartwert festlegen ( 891554152 ),
Anzahl der Starts ( 1000 ), Begriff hinzufügen ( { 1 , 0 } ), Begriff hinzufügen ( { 1 , 1 } ),
Begriff hinzufügen ( { 2 , 1 } ), Begriff hinzufügen ( { 3 , 1 } ), Begriff hinzufügen ( { 4 , 1 } ), Begriff hinzufügen ( { 5 , 1 } ),
Begriff hinzufügen ( { 1 , 1 }, { 2 , 1 } ), Begriff hinzufügen ( { 1 , 1 }, { 3 , 1 } ),
Begriff hinzufügen ( { 1 , 1 }, { 5 , 1 } ), Begriff hinzufügen ( { 2 , 1 }, { 3 , 1 } ),
Begriff hinzufügen ( { 2 , 1 }, { 4 , 1 } ), Begriff hinzufügen ( { 2 , 1 }, { 5 , 1 } ),
Begriff hinzufügen ( { 3 , 1 }, { 4 , 1 } ), Begriff hinzufügen ( { 3 , 1 }, { 5 , 1 } ),
Begriff hinzufügen ( { 4 , 1 }, { 5 , 1 } ), Alias-Begriff hinzufügen ( { 1 , 1 }, { 4 , 1 } ),
Setze N ganze Plots( 4 ), Stichprobengröße festlegen ( 28 ),
0 < S <= 0,5 | S > 0 & Betriebssystem == 3 | S == 0 & (Betriebssystem == 1 | Betriebssystem == 2 ))});
Betreff: DOE mit vielen Einschränkungen kann nicht eingerichtet werdenHallo Ryan,
Nein, habe ich nicht. Es scheint zu funktionieren, ich muss mich beim Kunden erkundigen, ob das gut genug und in Ordnung ist, um die diskreten Werte zu verwenden.
Eine Sache, die mir bei der Designbewertung und Leistungsanalyse aufgefallen ist, ich kann die Leistung nicht abschätzen, aber ich denke, das liegt an der diskreten Numerik, oder? Die Farbkarte sieht für mich ganz ok aus:
Der Bruchteil des Designraums sieht nach 0,6 seltsam aus:
Liegt das an den eingestellten Einschränkungen?
Danke,
Martin
Betreff: DOE mit vielen Einschränkungen kann nicht eingerichtet werdenIch denke, Sie stoßen auf Probleme bei der Verwendung der nicht zulässigen Kombinationen mit Gleichheiten bei kontinuierlichen Faktoren. Insbesondere wenn X==0 oder S==0 verwendet wird, kann etwas verwendet werden, das wirklich nahe an 0 liegt (sogar als 0 angezeigt). Haben Sie es mit diskreten Zahlen für die kontinuierliche versucht? Verwenden Sie eines Ihrer letzten Skripte, etwa so:
DAMHIRSCHKUH (
{Antwort hinzufügen( maximieren, "Y" , . , . , . ),
//Faktor hinzufügen (kontinuierlich, 0, 4, "S", 0 ),
//Faktor hinzufügen (kontinuierlich, 0, 0,6, "X", 0 ), Faktor hinzufügen (kontinuierlich, 0, 7, "M", 0 ),
Faktor hinzufügen( diskret numerisch, { 0 , 1 , 2 , 3 , 4 }, "S" , 0 ),
Faktor hinzufügen( diskret numerisch, { 0 , .fünfzehn , .3 , .45 , .6 }, "X" , 0 ),
Faktor hinzufügen( diskret numerisch, { 0 , 1,75 , 3.5 , 5.25 , 7 }, "M" , 0 ),
Faktor hinzufügen( Kategorial, { "L1" , "L2" , "Keiner" }, "Betriebssystem" , 0 ),
Faktor hinzufügen( Kategorial, { "L1" , "L2" }, "OCHSE" , 1 ), Zufallsstartwert festlegen ( 513965188 ),
Anzahl der Starts ( 1000 ), Begriff hinzufügen ( { 1 , 0 } ), Begriff hinzufügen ( { 1 , 1 } ),
Begriff hinzufügen ( { 2 , 1 } ), Begriff hinzufügen ( { 3 , 1 } ), Begriff hinzufügen ( { 4 , 1 } ), Begriff hinzufügen ( { 5 , 1 } ),
Begriff hinzufügen ( { 1 , 1 }, { 2 , 1 } ), Begriff hinzufügen ( { 1 , 1 }, { 3 , 1 } ),
Begriff hinzufügen ( { 1 , 1 }, { 4 , 1 } ), Begriff hinzufügen ( { 1 , 1 }, { 5 , 1 } ),
Begriff hinzufügen ( { 2 , 1 }, { 3 , 1 } ), Begriff hinzufügen ( { 2 , 1 }, { 4 , 1 } ),
Begriff hinzufügen ( { 2 , 1 }, { 5 , 1 } ), Begriff hinzufügen ( { 3 , 1 }, { 4 , 1 } ),
Begriff hinzufügen ( { 3 , 1 }, { 5 , 1 } ), Begriff hinzufügen ( { 4 , 1 }, { 5 , 1 } ), setze N ganze Plots( 4 ),
Stichprobengröße festlegen ( 28 ), unzulässige Kombinationen (
X< 0,05 | 0 < S <= 0,5 | S > 0 & Betriebssystem == 3 | S == 0 & (Betriebssystem == 1 | Betriebssystem == 2 )
Betreff: DOE mit vielen Einschränkungen kann nicht eingerichtet werdenHallo Martin,
Haben Sie versucht, das Skript für nicht zulässige Kombinationen zu verwenden? Dies gibt Ihnen mehr Flexibilität beim Definieren von Einschränkungen.
http://www.jmp.com/support/help/Define_Factor_Constraints.shtml#205355
Grüße,
Phil
