Come posso simulare il modello di esecuzione di prova per DoE virtuale?Potresti per favore aiutarmi con questo problema?
Voglio creare un emulatore di computer per scopi didattici. Emulerà un risultato binario e userà diversi fattori Sì / No come input (diciamo ABC).
L'output sarà la% dei successi.
L'output quando tutti i fattori sono a "-" è il 5% dei successi.
quando
A a "+" 5,5% (guadagno dello 0,5%)
B a "+" 6,5% (guadagno 1,5%)
C a "+" 5% (guadagno 0%)
A e B a "+" (interazione AB) 7,5% (guadagno del 2,5%)
Quindi voglio implementarlo in un software come un modello di esecuzione di prova con l'elenco di prove e uscite casuali (a seconda dei livelli di fattori ABC).
Esempio di comportamento del software:
Ordine di prova ABC Output
1a prova + + + fallimento
2 ° prova + - + Fail
3 ° prova + - - Successo
.....
10000 prova + - - Fail
Voglio che appaia realistico in modo che non sempre la combinazione di A + B- C- possa risultare in "Successo" perché il processo non è fisico, è una risposta umana.
Qual è la procedura generale per costruire un tale algoritmo?
Può essere fatto in JMP?
Inoltre sto prendendo in considerazione strumenti come Excel, R, software di simulazione come Anylogic, Promodel è necessario.
Grazie tante!
Ri: Come posso simulare il modello di esecuzione di prova per DoE virtuale?Scusa, forse dovrei creare un nuovo argomento.
Tuttavia, è possibile raggiungere il mio obiettivo utilizzando lo script Simulate Response costruito in JMP?
Ri: Come posso simulare il modello di esecuzione di prova per DoE virtuale?Grazie mille per il vostro aiuto con la simulazione!
Potresti aiutarmi per favore con la modifica di questo script o suggerendo un nuovo approccio per questo compito?
Ora voglio simulare un progetto di screening definitivo. Ci sono 18 fattori.
Dimensioni degli effetti:
| Fattore | % di successi |
| UN | 1,31% |
| B | 1,42% |
| D | 1,58% |
| E | 1,64% |
| F | 1,46% |
| H | 1,39% |
| K | 1,34% |
| M | 1,37% |
| | |
| DE | 1,67% |
| FK | 1,54% |
Quando tutti i fattori con un tasso di successo a livello "0" basso sono pari all'1,2%.
Ecco il design esatto che voglio simulare:
| X1 | X2 | X3 | X4 | X5 | X6 | X7 | X8 | X9 | X10 | X11 | X12 | X13 | X14 | X15 | X16 | X17 | X18 |
| 0 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 1 |
| 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 |
| 0 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 |
| 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 |
| 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 |
| 1 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 |
| 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 |
| 0 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 |
| 1 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 |
| 0 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 |
| 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| 1 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 |
| 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 |
| 0 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 1 |
| 0 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 1 |
| 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 |
| 0 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 |
| 1 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 |
| 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 |
| 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| 0 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 |
| 0 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 |
| 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 |
| 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 |
| 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 |
| 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
| 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 |
| 1 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 |
| 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 |
| 1 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 |
| 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 |
| 0 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 1 | 0 |
| 1 | 0 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 |
| 1 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 1 | 0 | 0 |
Ri: Come posso simulare il modello di esecuzione di prova per DoE virtuale?Si stanno adattando modelli che sono più complessi di quanto questo calcolatore di dimensioni del campione intendesse fare. Ci sono ulteriori gradi di libertà usati per stimare altri termini.
Ti suggerirei di tornare alle formule per ciascuno degli 8 diversi scenari (combinazioni di A, B e C) per assicurarti di utilizzare le probabilità appropriate. Hai fornito solo quattro delle otto situazioni, quindi ho fatto delle ipotesi sugli altri quattro casi. Potresti voler controllare quelli per accertarti che tutto sia come ti aspettavi. Immagino che ti aspetti che l'effetto di interazione sia un guadagno del 2,5%, ma il termine di interazione NON è un guadagno del 2,5%. Ottieni un guadagno del 2,5% solo quando sia A sia B sono positivi come hai indicato. L'interazione include anche ciò che accade per le altre 3 combinazioni di A e B.
Ad esempio: A e B entrambi positivi: la probabilità è del 7,5%
A e B entrambi zero: la probabilità è del 5%. Quindi quando sono "lo stesso segno" danno una probabilità media del 6,25%
A zero, B positivo: la probabilità è del 6,5%
A positivo, B zero: la probabilità è del 5,5% per una probabilità media del 6%.
Pertanto, l'effetto di interazione è in realtà dello 0,25% che è molto piccolo (e molto vicino alla stima dell'interazione A * B nell'analisi).
Ri: Come posso simulare il modello di esecuzione di prova per DoE virtuale?Molte grazie per i chiarimenti!
Ho eseguito il calcolo della dimensione minima del campione in modo che potesse rilevare la differenza dello 0,5%. Tuttavia, a volte l'analisi dei risultati simulati non rileva l'interazione AB come significativa con 64000 studi.
Ho fatto i calcoli, giusto?
Ri: Come posso simulare il modello di esecuzione di prova per DoE virtuale?Naturalmente ogni simulazione è diversa, ma stai cercando dimensioni di effetti molto piccole. Con 10.000 prove, avrai circa 5000 prove con A a 1 e 5000 prove con A a 0. Di quelli 5000 a 0, ci si aspetta un tasso di "successo" del 5% che è di circa 250 osservazioni. Con 5000 a 1, si prevede una percentuale di successo del 5,5%, ovvero circa 275 osservazioni. Solo 25 differenze di osservazione su 10.000 totali. Questa è una piccola differenza. Aggiungi il rumore casuale a questo e puoi vedere che la differenza potrebbe effettivamente essere ancora più piccola.
I dati categoriali richiedono in genere dimensioni del campione maggiori. Ciò, unito alla ricerca di effetti molto piccoli, porta a non significanza. Indicato in termini più statistici, si ha una potenza molto bassa per rilevare le dimensioni dell'effetto desiderato.
Ri: Come posso simulare il modello di esecuzione di prova per DoE virtuale?Grazie mille per questo script! Questo è quello per cui stavo cercando.
Potresti dirmi per favore il motivo per cui quando ho generato 10000 righe la regressione Logistica Nominale non mostra che A e AB sono effetti significativi?
Ri: Come posso simulare il modello di esecuzione di prova per DoE virtuale?Questo potrebbe probabilmente essere fatto in modo più elegante, ma un rapido "approccio forza bruta" sarebbe quello di inserire le formule casuali in una tabella di dati JMP. Ho reso le colonne A, B e C una colonna di numeri casuali tra 0 e 1 (ho scelto 0 invece di -1 per facilità). Quindi per la colonna Risultato, ho appena inserito una stringa di condizioni "If" con le varie probabilità da una distribuzione binomiale con la probabilità di successo appropriata. Abbastanza semplice in questo caso poiché ci sono solo 8 combinazioni possibili di A, B e C. Per un problema più grande, proverei ad essere più elegante. Il numero di righe nella tua tabella sarà il numero di prove casuali che desideri avere. Per questo semplice caso ho appena avuto 10 righe.
New Table( "Untitled",
Add Rows( 10 ),
New Column( "A", Numeric, "Nominal", Format( "Best", 12 ),
Formula( Random Integer( 0, 1 ) ) ),
New Column( "B", Numeric, "Nominal", Format( "Best", 12 ),
Formula( Random Integer( 0, 1 ) ) ),
New Column( "C", Numeric, "Nominal", Format( "Best", 12 ),
Formula( Random Integer( 0, 1 ) ) ),
New Column( "Success?", Numeric, "Nominal", Format( "Best", 12 ),
Formula(
If(
:A == 0 & :B == 0 & :C == 0, Random Binomial( 1, 0.05 ),
:A == 0 & :B == 0 & :C == 1, Random Binomial( 1, 0.05 ),
:A == 0 & :B == 1 & :C == 0, Random Binomial( 1, 0.065 ),
:A == 0 & :B == 1 & :C == 1, Random Binomial( 1, 0.065 ),
:A == 1 & :B == 0 & :C == 0, Random Binomial( 1, 0.055 ),
:A == 1 & :B == 0 & :C == 1, Random Binomial( 1, 0.055 ),
:A == 1 & :B == 1 & :C == 0, Random Binomial( 1, 0.075 ),
:A == 1 & :B == 1 & :C == 1, Random Binomial( 1, 0.075 ),
Random Binomial( 1, 0.05 ) ) ) ) )
