Hola, equipo,
Uno de los profesores que enseña DOE usando JMP está usando el libro Experiments with Mixtures: Designs, Models and Analysis of Mixture data. de la publicación de Wiley. Está intentando replicar el resultado y el análisis del libro.
He adjuntado el conjunto de datos aquí.
Por favor aconséjame
gracias
Murali
Supongo que está preguntando por las correlaciones entre las estimaciones de los parámetros en la regresión. Agregué un script de tabla para el cuadro de diálogo de inicio Ajustar modelo y la plataforma Ajustar mínimos cuadrados. Haga clic en el triángulo rojo en la parte superior y seleccione Estimaciones > Correlación de estimaciones.
Normalmente no preguntamos sobre las correlaciones entre los componentes.
Vea la tabla de datos adjunta que modifiqué.
Tenga en cuenta que también puede abrir la tabla de datos del experimento y luego seleccionar DOE > Diagnóstico de diseño > Evaluar diseño. Este comando producirá los diagnósticos de diseño, incluidas las correlaciones entre estimaciones.
(10/11/23) Me comuniqué con Ryan Lekivetz para verificar mi explicación. Proporcionó un guión que demuestra la diferencia entre las correlaciones reportadas en la sección Evaluación de diseño en Diseño personalizado y la tabla Correlación de estimaciones en Ajustar mínimos cuadrados. El primero informa las correlaciones por pares entre las columnas de la matriz del modelo. Este último tiene en cuenta la procedimiento de regresión .
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 ) );
Muchas gracias @Mark_Bailey .
El profesor estaba interesado en la matriz (X'X)-1 sigma^2 real. Creo que JMP (el profesor también piensa lo mismo) tiene disposiciones para la matriz X'X que, sin embargo, debe invertirse.
Por favor aconséjame
gracias
murali
