Wie kann ich ein Testausführungsmodell für virtuelles DoE simulieren?Könnten Sie mir bitte bei diesem Problem helfen?
Ich möchte einen Computeremulator für Lernzwecke erstellen. Es emuliert ein binäres Ergebnis und verwendet mehrere Ja / Nein-Faktoren als Eingabe (beispielsweise ABC).
Die Ausgabe erfolgt in% der Erfolge.
Die Ausgabe, wenn alle Faktoren auf "-" stehen, beträgt 5% der Erfolge.
Wann
A bei "+" 5,5% (0,5% Gewinn)
B bei "+" 6,5% (1,5% Gewinn)
C bei "+" 5% (0% Gewinn)
A und B bei "+" (AB-Wechselwirkung) 7,5% (2,5% Gewinn)
Daher möchte ich es in einer Software als Testausführungsmodell mit der Liste der zufälligen Versuche und Ausgaben implementieren (abhängig von den Faktoren der ABC-Faktoren).
Beispiel für das Verhalten der Software:
Reihenfolge der Testversion von ABC
1. Versuch + + + Fehlschlagen
2. Versuch + - + Fail
3. Versuch + - - Erfolg
.....
10000 Versuch + - - Fehler
Ich möchte, dass es realistisch aussieht, so dass die Kombination von A + B- C- nicht immer zu "Erfolg" führt, weil der Prozess nicht physisch ist, sondern eine menschliche Reaktion.
Wie sieht der allgemeine Aufbau eines solchen Algorithmus aus?
Kann es in JMP gemacht werden?
Darüber hinaus erwäge ich Werkzeuge wie Excel, R, Simulationssoftware wie Anylogic, Promodel.
Vielen Dank!
Re: Wie kann ich ein Testausführungsmodell für virtuelles DoE simulieren?Sorry, vielleicht sollte ich ein neues Thema erstellen.
Ist es jedoch möglich, mein Ziel mithilfe des in JMP integrierten Simulate Response-Skripts zu erreichen?
Re: Wie kann ich ein Testausführungsmodell für virtuelles DoE simulieren?Vielen Dank für Ihre Hilfe bei der Simulation!
Könnten Sie mir bitte helfen, dieses Skript zu ändern oder einen neuen Ansatz für diese Aufgabe vorzuschlagen?
Jetzt möchte ich ein definitives Screening-Design simulieren. Es gibt 18 Faktoren.
Effektgrößen:
| Faktor | % der Erfolge |
| EIN | 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% |
Wenn alle Faktoren auf niedrigem "0" -Niveau liegen, beträgt die Erfolgsquote 1,2%.
Hier ist das genaue Design, das ich simulieren möchte:
| 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 |
Re: Wie kann ich ein Testausführungsmodell für virtuelles DoE simulieren?Sie passen Modelle an, die komplexer sind, als dies mit dem Rechner für die Stichprobengröße der Fall war. Es werden zusätzliche Freiheitsgrade verwendet, um andere Begriffe zu schätzen.
Ich würde vorschlagen, dass Sie zu den Formeln für jedes der 8 verschiedenen Szenarien (Kombinationen von A, B und C) zurückkehren, um sicherzustellen, dass Sie die richtigen Wahrscheinlichkeiten verwenden. Sie hatten nur vier der acht Situationen angegeben, daher habe ich einige Annahmen über die anderen vier Fälle gemacht. Sie können diese überprüfen, um sicherzustellen, dass alles so ist, wie Sie es erwartet haben. Ich vermute, Sie erwarten einen Interaktionseffekt von 2,5%, aber der Interaktionstermin beträgt NICHT 2,5%. Sie erhalten nur 2,5% Gewinn, wenn sowohl A als auch B wie angegeben positiv sind. Die Interaktion beinhaltet auch, was für die anderen 3 Kombinationen von A und B geschieht.
Als Beispiel: A und B beide positiv: Wahrscheinlichkeit ist 7,5%
A und B beide Null: Die Wahrscheinlichkeit beträgt 5%. Wenn sie also "dasselbe Zeichen" sind, ergibt sich eine durchschnittliche Wahrscheinlichkeit von 6,25%.
A null, B positiv: Wahrscheinlichkeit 6,5%
A positiv, B null: Die Wahrscheinlichkeit beträgt 5,5% bei einer durchschnittlichen Wahrscheinlichkeit von 6%.
Daher ist der Interaktionseffekt tatsächlich 0,25%, was sehr gering ist (und ziemlich nahe an Ihrer A * B-Interaktionsschätzung in Ihrer Analyse liegt).
Re: Wie kann ich ein Testausführungsmodell für virtuelles DoE simulieren?Vielen Dank für die Klarstellung!
Ich habe die Berechnung der minimalen Stichprobengröße durchgeführt, damit der Unterschied von 0,5% erkannt werden konnte. Bei der Analyse der simulierten Ergebnisse wird die AB-Wechselwirkung jedoch manchmal bei 64000 Studien nicht als signifikant erkannt.
Habe ich die Berechnungen richtig gemacht?
Re: Wie kann ich ein Testausführungsmodell für virtuelles DoE simulieren?Natürlich ist jede Simulation anders, aber Sie suchen nach sehr kleinen Effektgrößen. Bei 10.000 Versuchen haben Sie ungefähr 5000 Versuche mit A bei 1 und 5000 Versuche mit A bei 0. Von diesen 5000 bei 0 erwarten Sie eine Erfolgsquote von 5%, etwa 250 Beobachtungen. Bei 5000 bei 1 erwarten Sie eine Erfolgsquote von 5,5%, was etwa 275 Beobachtungen entspricht. Nur ein Unterschied von 25 Beobachtungen von insgesamt 10.000. Das ist ein kleiner Unterschied. Fügen Sie das zufällige Rauschen hinzu und Sie können sehen, dass der Unterschied tatsächlich noch geringer sein kann.
Kategoriale Daten erfordern normalerweise größere Stichprobengrößen. Das, verbunden mit der Suche nach sehr kleinen Effekten, führt zur Nicht-Signifikanz. In statistischen Worten ausgedrückt, haben Sie eine sehr geringe Leistung, um die gewünschten Effektgrößen zu ermitteln.
Re: Wie kann ich ein Testausführungsmodell für virtuelles DoE simulieren?Vielen Dank für dieses Skript! Dafür habe ich gesucht.
Könnten Sie mir bitte den Grund nennen, warum die generelle logistische Regression bei der Generierung von 10000 Zeilen nicht zeigt, dass A und AB signifikante Auswirkungen haben?
Re: Wie kann ich ein Testausführungsmodell für virtuelles DoE simulieren?Dies könnte wahrscheinlich eleganter erfolgen, aber ein schneller "Brute-Force-Ansatz" wäre die Eingabe der Zufallsformeln in eine JMP-Datentabelle. Ich habe die A-, B- und C-Spalten zu einer zufälligen ganzzahligen Spalte zwischen 0 und 1 gemacht (ich habe zur Vereinfachung 0 anstelle von -1 gewählt). Dann habe ich für die Ergebnisspalte gerade eine Folge von "If" -Bedingungen mit den verschiedenen Wahrscheinlichkeiten aus einer Binomialverteilung mit der entsprechenden Erfolgswahrscheinlichkeit eingegeben. Ziemlich einfach in diesem Fall, da es nur 8 mögliche Kombinationen von A, B und C gibt. Für ein größeres Problem würde ich versuchen, eleganter zu sein. Die Anzahl der Zeilen in Ihrer Tabelle ist die Anzahl der Stichproben, die Sie wünschen. Für diesen einfachen Fall hatte ich gerade 10 Reihen.
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 ) ) ) ) )
