Salut l'équipe,
L'un des professeurs qui enseigne le DOE à l'aide de JMP utilise le livre Experiments with Mixtures: Designs, Models and Analysis of Mixture data. de la publication Wiley. Il essaie de reproduire le résultat et l’analyse du livre.
J'ai joint l'ensemble de données ici.
S'il vous plaît, des conseils
merci
Peinture murale
Je suppose qu'il demande les corrélations entre les estimations des paramètres dans la régression. J'ai ajouté un script de table pour la boîte de dialogue de lancement Ajuster le modèle et la plateforme Ajuster les moindres carrés. Cliquez sur le triangle rouge en haut et sélectionnez Estimations > Corrélation des estimations.
Nous ne posons généralement pas de questions sur les corrélations entre les composants.
Voir le tableau de données ci-joint que j'ai modifié.
Notez que vous pouvez également ouvrir la table de données de l'expérience, puis sélectionner DOE > Design Diagnostics > Evaluate Design. Cette commande produira les diagnostics de conception, y compris les corrélations entre les estimations.
(10/11/23) J'ai contacté Ryan Lekivetz pour vérifier mon explication. Il a fourni un script démontrant la différence entre les corrélations rapportées dans la section Évaluation de la conception dans Conception personnalisée et le tableau Corrélation des estimations dans Ajuster les moindres carrés. Le premier rapporte les corrélations par paires entre les colonnes de la matrice du modèle. Cette dernière prend en compte le procédure de régression .
Names Default To Here( 1 );
// use the open data table for example
xmat = Current Data Table() << Get As Matrix;
// make another empty, conforming matrix
xmatint = Matrix( N Row( xmat ), 6 );
// copy factor levels
xmatint[0, 1 :: 3] = xmat[0, 1 :: 3];
// create cross terms
xmatint[0, 4] = xmat[0, 1] :*</img> xmat[0, 2];
xmatint[0, 5] = xmat[0, 1] :*</img> xmat[0, 3];
xmatint[0, 6] = xmat[0, 2] :*</img> xmat[0, 3];
//evaluate design version (Color Map on Correlations)
Show( Correlation( xmatint ) );
//Fit Model version
V = Sqrt( Diag( Inverse( Transpose( xmatint ) * xmatint ) ) );
Show( Inverse( V ) * Inverse( Transpose( xmatint ) * xmatint ) * Inverse( V ) );
Merci beaucoup @Mark_Bailey .
Le professeur s'intéressait à la matrice réelle (X'X)-1 sigma^2. Je pense que JMP (le professeur pense aussi la même chose) a des dispositions pour la matrice X'X qui doit cependant être inversée.
S'il vous plaît, des conseils
merci
murali
