Ciao,
Ho i seguenti fattori e livelli per i quali sto cercando di creare un design D-ottimale:
Fattore A - Categoria due livelli (Jolt, Triad)
Fattore B - Categorico a due livelli (Blunt, ventosa)
Fattore C - Continuo (10,25)
Fattore D - Categorico a due livelli (in piedi, seduto)
Fattore E - Categorico tre livelli (Stazionario, Schivare, Correre)
Fattore F - Categoria due livelli (Principiante, Esperto)
Sto tentando di creare un esempio che gli studenti possano utilizzare per insegnare loro come generare progetti ottimali con combinazioni non consentite. La combinazione non consentita che sto cercando di creare è Factor D (Seated) e Factor F (Expert). Sto anche cercando di costruire un progetto con effetti principali E interazioni a due fattori. Tuttavia, continuo a ricevere il messaggio di errore (il designer ottimale non è riuscito a convergere). Ho provato tutto il possibile per generare un design che converge con combinazioni non consentite tra fattori categorici senza alcun risultato. Ho provato ad aggiungere fattori, aggiungendo più livelli a questi fattori. Posso tuttavia convincere il progettista a convergere quando aggiunge una combinazione non consentita tra un fattore categorico e un fattore continuo. Il secondo che aggiungo in una combinazione non consentita tra fattori categorici l'algoritmo non convergerà. Qualsiasi aiuto è apprezzato!
Ciao @FaceStatistics6 ,
Riproduco le stesse impostazioni e riscontro un messaggio di errore simile.
Il problema potrebbe risiedere nella tua combinazione di combinazioni e termini non consentiti nel tuo modello: se blocchi la combinazione di "Seated" da Factor D (categoria a 2 livelli) con "Expert" da Factor F (categoria a 2 livelli), ciò significa non puoi stimare il termine di interazione di questi due fattori poiché un vincolo ti impedisce di testare tutte le possibili combinazioni di livelli per questi due fattori.
Puoi:
- Rimuovere il termine di interazione a 2 fattori "Factor D*Factor F" e utilizzare questo script per generare la tabella DoE (o utilizzare la tabella dati allegata "Custom Design_Remove2FI"):
DOE(
Custom Design,
{Add Response( Maximize, "Y", ., ., . ),
Add Factor( Categorical, {"Jolt", "Triad"}, "Factor A", 0 ),
Add Factor( Categorical, {"Blunt", "Suction Cup"}, "Factor B", 0 ),
Add Factor( Categorical, {"Standing", "Seated"}, "Factor D", 0 ),
Add Factor( Categorical, {"Novice", "Expert"}, "Factor F", 0 ),
Add Factor( Continuous, 10, 25, "Factor C", 0 ),
Add Factor( Categorical, {"Stationary", "Evading", "Running"}, "Factor E", 0 ),
Set Random Seed( 80459262 ), Number of Starts( 2161 ), Add Term( {1, 0} ),
Add Term( {1, 1} ), Add Term( {2, 1} ), Add Term( {3, 1} ), Add Term( {4, 1} ),
Add Term( {5, 1} ), Add Term( {6, 1} ), Add Term( {1, 1}, {2, 1} ),
Add Term( {1, 1}, {3, 1} ), Add Term( {1, 1}, {4, 1} ),
Add Term( {1, 1}, {5, 1} ), Add Term( {1, 1}, {6, 1} ),
Add Term( {2, 1}, {3, 1} ), Add Term( {2, 1}, {4, 1} ),
Add Term( {2, 1}, {5, 1} ), Add Term( {2, 1}, {6, 1} ),
Add Term( {3, 1}, {5, 1} ), Add Term( {3, 1}, {6, 1} ),
Add Term( {4, 1}, {5, 1} ), Add Term( {4, 1}, {6, 1} ),
Add Term( {5, 1}, {6, 1} ), Set Sample Size( 36 ),
Disallowed Combinations( Factor D == "Seated" & Factor F == "Expert" ),
Simulate Responses( 0 ), Save X Matrix( 0 ), Make Design}
);
- Oppure imposta la stimabilità dell'interazione dei 2 fattori "Factor D*Factor F" come "If Possible" e usa questo script per generare la tabella DoE (o usa la tabella dati allegata "Custom Design_IfPossible"):
DOE(
Custom Design,
{Add Response( Maximize, "Y", ., ., . ),
Add Factor( Categorical, {"Jolt", "Triad"}, "Factor A", 0 ),
Add Factor( Categorical, {"Blunt", "Suction Cup"}, "Factor B", 0 ),
Add Factor( Categorical, {"Standing", "Seated"}, "Factor D", 0 ),
Add Factor( Categorical, {"Novice", "Expert"}, "Factor F", 0 ),
Add Factor( Continuous, 10, 25, "Factor C", 0 ),
Add Factor( Categorical, {"Stationary", "Evading", "Running"}, "Factor E", 0 ),
Set Random Seed( 349202301 ), Number of Starts( 2106 ), Add Term( {1, 0} ),
Add Term( {1, 1} ), Add Term( {2, 1} ), Add Term( {3, 1} ), Add Term( {4, 1} ),
Add Term( {5, 1} ), Add Term( {6, 1} ), Add Term( {1, 1}, {2, 1} ),
Add Term( {1, 1}, {3, 1} ), Add Term( {1, 1}, {4, 1} ),
Add Term( {1, 1}, {5, 1} ), Add Term( {1, 1}, {6, 1} ),
Add Term( {2, 1}, {3, 1} ), Add Term( {2, 1}, {4, 1} ),
Add Term( {2, 1}, {5, 1} ), Add Term( {2, 1}, {6, 1} ),
Add Term( {3, 1}, {5, 1} ), Add Term( {3, 1}, {6, 1} ),
Add Term( {4, 1}, {5, 1} ), Add Term( {4, 1}, {6, 1} ),
Add Term( {5, 1}, {6, 1} ), Add Potential Term( {3, 1}, {4, 1} ),
Set Sample Size( 36 ), Disallowed Combinations(
Factor D == "Seated" & Factor F == "Expert"
), Simulate Responses( 0 ), Save X Matrix( 0 ), Make Design}
);
Spero che questa risposta ti possa aiutare,