Liebe Community,
Ich bin Professor und integriere (langsam) neue Inhalte von JMP Pro in meine Chemievorlesungen. Ein Thema, das ich behandle, ist der Tukey-Test mit einem Freiheitsgrad auf Nichtadditivität bei nicht replizierten faktoriellen Versuchsreihen. Gibt es Alternativen, um diesen Test mit JMP Pro durchzuführen (ich habe online nichts dazu gefunden)?
Danke
Rita
Ich glaube nicht, dass es eine integrierte Möglichkeit gibt, diesen Test zu erhalten, aber Sie können JSL verwenden, um ihn zu erstellen. Die Schritte sind: Passen Sie Ihr 2-Weg-ANOVA-Modell an. Speichern Sie die vorhergesagten Werte. Quadratiere sie. Passen Sie nun das 2-Wege-Modell wieder an, aber fügen Sie die quadrierten Werte als Spalte hinzu. Dieser Test ist Ihr nicht-additiver Test. Ich habe einen Beispieldatensatz mit einem Skript beigefügt, das dies tut. Hier ist das Skript, das ich verwendet habe und das für dieses Beispiel fest codiert ist.
Names Default to Here(1);
dt = Current Data Table();
// Launch platform: Fit Least Squares
dt = Current Data Table();
Mod1 = dt << Fit Model(
Y( :Impurity ),
Effects( :Temp, :Pressure ),
Personality( "Standard Least Squares" ),
Emphasis( "Effect Leverage" ),
Run(
:Impurity << {Summary of Fit( 1 ), Analysis of Variance( 1 ),
Parameter Estimates( 1 ), Lack of Fit( 0 ), Scaled Estimates( 0 ),
Plot Actual by Predicted( 1 ), Plot Regression( 0 ),
Plot Residual by Predicted( 1 ), Plot Studentized Residuals( 0 ),
Plot Effect Leverage( 1 ), Plot Residual by Normal Quantiles( 0 ),
Box Cox Y Transformation( 0 )}
), Invisible
);
mod1 << Prediction Formula;
mod1 << Close Window();
// New formula column: Pred Formula Impurity^2
dt << New Formula Column(
Operation( Category( "Transform" ), "Square" ),
Columns( :Pred Formula Impurity )
);
// Change column name: Pred Formula Impurity^2 → Nonadd
dt:"Pred Formula Impurity^2"n << Set Name( "Non-Additive Test" );
// Launch platform: Fit Least Squares
dt << Fit Model(
Y( :Impurity ),
Effects( :Temp, :Pressure, :"Non-Additive Test"n ),
Personality( "Standard Least Squares" ),
Emphasis( "Effect Leverage" ),
Run(
:Impurity << {Summary of Fit( 0 ), Analysis of Variance( 0 ),
Parameter Estimates( 1 ), Sequential Tests( 1 ), Lack of Fit( 0 ), Scaled Estimates( 0 ),
Effect Tests( 0 ),
Plot Actual by Predicted( 0 ), Plot Regression( 0 ),
Plot Residual by Predicted( 0 ), Plot Studentized Residuals( 0 ),
Plot Effect Leverage( 0 ), Plot Residual by Normal Quantiles( 0 ),
Box Cox Y Transformation( 0 )}
)
);