<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Ich bin neugierig, ob jemand weiß und kann mir erklären, wie JMP seinen Nicht-Zentralitätsparameter für benutzerdefinierte Designs berechnet. Ich kenne das allgemeine Formular für die Berechnung, aber ich kann meine Berechnungen nicht ganz auf die Ergebnisse von JMP beziehen. Der Grund dafür ist, dass der df Fehler bei der Berechnung der Leistung für Split-Plot-Designs sind in JMP 10 standardmäßig falsch. Während ich das df einstellen kann manuell wäre es schön, wenn ich ein Skript schreiben könnte, um die Berechnungen von Anfang an korrekt auszuführen, aber die Plattform für die Entwurfsbewertung ist in JMP 10 nicht skriptfähig (nicht dass ich es trotzdem herausfinden kann - auch wenn ich ein << Get Scriptable Object) probiere ein leeres Set zurückgegeben) Ich muss also in JSL die Berechnungen "von Hand" ausführen können.
Danke im Voraus!
Dr. Jones,
Danke für Ihre Antwort. Möchten Sie erklären, wie dies bei einer Split-Plot-Struktur funktioniert? Ich habe mit dem benutzerdefinierten Designer einen einfachen 5-Faktor (2 HTC-, 3 ETC-Faktoren) erstellt und ein Design mit 6 WP und 8 Replikaten in jedem WP erhalten, um ein perfekt repliziertes Design zu erstellen. Ich passe das Modell mit REML an und speicherte dann die Codierungstabelle, um die X-Matrix zu erhalten. Ich kann X'X jedoch nicht mit der Kodierung lösen, die für die 6 ganzen Plots verwendet wird, da sie (absichtlich) mit den HTC-Faktoren verwechselt werden. Wenn ich für die ETC-Faktoren wie angegeben vorkomme, erhalte ich dasselbe ncp wie JMP. Wenn ich die HTC-Faktoren auf nur die 6 unabhängigen Einstellungen reduziere, um sie so zu behandeln, als gäbe es nur 6 Läufe (was tatsächlich der Fall ist), bekomme ich sogar 6 für den ncp. JMP gibt mir jedoch 5.333 und 4.741 (die genauen Werte sind nicht wichtig). Es passiert also noch etwas anderes und ich bin mir nicht sicher, was es ist oder wie ich es erklären muss.
Ich freue mich, dass Sie sich die Zeit genommen haben, den oben genannten Code zur Verfügung zu stellen, um den NCP zu erhalten.
Bei den Leistungsberechnungen gibt es einen Unterschied zwischen JMP9 und JMP10. In JMP9 wurde die Leistung für einen Signal-zu-Rausch-Abstand berechnet, wobei der Regressionskoeffizient das Signal und die Fehlerstandardabweichung das Rauschen war. In JMP 10 ist das Signal für einen 2-Stufen-Faktor die Änderung der Antwort, wenn vom niedrigen Wert des Faktors auf den hohen Wert des Faktors übergegangen wird. Das Signal ist also 2-fach so groß wie der Regressionskoeffizient für einen Effekt von 1 Freiheitsgrad.
Hier ist die JSL, die Sie benötigen, um die Leistungsberechnungen in JMP 9 für einen Effekt von 1 Freiheitsgrad zu berechnen
Sei X die X-Matrix für das Modell.
Dann ist die Varianz der Koeffizienten
varBetaHat = vecdiag (invers (X` * X));
alpha = 0,05; // Standardwert der Signifikanz
NennerDF = 11; // df für ein 16-maliges orthogonales Design in 4 Faktoren
Fcrit = F Quantil (1-alpha, 1, NennerDF); // kritischer Wert der F-Verteilung
SNRatio = 1; // Standard-SNRatio
Nicht-Zentralität = (SNRatio ^ 2) / varBetaHat [2]; // Nichtzentralität für den 1. Haupteffekt unter der Annahme, dass der 1. Effekt der Intercept ist
Leistung = 1 - F Verteilung (Fcrit, 1, NennerDF, Nichtzentrale);
Hier ist die JSL, die Sie benötigen, um die Leistungsberechnungen in JMP 10 für einen Effekt mit 1 Freiheitsgrad unter Verwendung der gleichen X-Matrix zu berechnen.
varContrast = 4 * vecdiag (invers (X` * X)); // Da der Kontrast doppelt so groß ist wie der Koeffizient, ist seine Varianz viermal so groß.
alpha = 0,05; // Standardwert der Signifikanz
NennerDF = 11; // df für ein 16-maliges orthogonales Design in 4 Faktoren
Fcrit = F Quantil (1-alpha, 1, NennerDF); // kritischer Wert der F-Verteilung
SNRatio = 2; // default SNRatio ist doppelt so groß, dass es den gleichen Standardwert für die Leistung liefert.
Nicht-Zentralität = (SNRatio ^ 2) / varContrast [2]; // Nichtzentralität für den 1. Haupteffekt unter der Annahme, dass der 1. Effekt der Intercept ist
Leistung = 1 - F Verteilung (Fcrit, 1, NennerDF, Nichtzentrale);
