63200_es-ES

liliconcom

Opción Guardar columnas con modelo mixto

Hola todos,

 

Cuando ejecutamos un modelo mixto utilizando la plataforma Fit Model, podemos guardar valores residuales y valores predichos, pero también tenemos la posibilidad de guardar valores residuales condicionales y valores predecibles condicionales (que tienen en cuenta los efectos fijos y aleatorios). Sin embargo, no tenemos opciones similares para los residuos estudiados y las columnas StdErr Pred Formula.

¿Cuál es la mejor manera de obtener o tal vez volver a calcular la "verificación condicional" de estas dos columnas?

 

Tenga en cuenta que no estoy trabajando con JMP Pro.

 

¡Gracias de antemano por tu ayuda!

 

 

Re: Guardar opción de columnas con modelo mixto

Sí, eso es correcto, voy a hacer eso!

Re: Guardar opción de columnas con modelo mixto

Puede sugerir esta función en la parte de la Lista de deseos de la comunidad JMP.

Re: Guardar opción de columnas con modelo mixto

Hola GM

Sí, tienes razón, podemos extraer los errores estándar de efectos aleatorios como ese, pero lo que estoy buscando son

Los errores estándar de residuos y de los valores previstos. No creo que sea la misma información aquí.

Re: Guardar opción de columnas con modelo mixto

Gracias Cameron por esta respuesta detallada, la probaré!

Pero, de hecho, no es una forma fácil de obtener lo que queremos. Me pregunto por qué no está disponible directamente desde el menú Guardar columnas.

Re: Guardar opción de columnas con modelo mixto

Si estima los efectos aleatorios utilizando la plataforma Fit Model con Personality = mínimos cuadrados estándar y Method = REML, puede extraer los errores estándar de los efectos aleatorios haciendo una tabla de datos fuera de la tabla de Predicción de efectos aleatorios; se incluyen los errores estándar.   La personalidad del modelo mixto no genera los errores estándar con los coeficientes de efectos aleatorios en el informe.

 

Re: Guardar opción de columnas con modelo mixto

Tienes razón.  Duh  No tengo acceso a JMP Pro, por lo que no puedo verificarlo, pero no creo que pueda obtener todos los componentes que necesitaría para calcular los residuos condicionales estandarizados directamente del informe.  La columna PredSE proporciona la matriz X como el segundo argumento de Vec Quadratic (), y el informe le proporciona el MSE.  Tendrías que construir G y Z. 

 

Podrías obtener Z de algo como:

Design Nom( :RandomEffectCol << Get Values, {list of levels})

 Si tiene múltiples efectos aleatorios, podría simplemente concatenarlos juntos.

Para obtener G, eso dependería de la estructura de correlación elegida.  Si tiene una estructura de correlación simple con todos los 0 en las diagonales fuera de línea, puede construir G multiplicando su vector de componentes de varianza con una matriz de identidad.

Una vez que tenga X, G, Z y el MSE, podría calcular la varianza de la matriz de error de la siguiente manera:

n = N Row(dt);

V = Z*G*Z` + MSE*Identity(n);
K = Identity(n) - Z*G*Z`*Inv(V);
Q = X*Inv(X`*Inv(V)*X)*X`;

Var_e = K*(V-Q)*K`;

A partir de ahí, tomaría el cuadrado de los elementos diagonales de Var_e para obtener los errores estándar condicionales para sus residuales condicionales.  No creo que JMP tenga una función para devolver solo los elementos diagonales de la matriz, pero aquí hay uno simple que hará el truco:

getDiag = function({X},
 n = ncol(X);
 diags = {};
 for(i=1, i<=N Col(X),i++,
  Insert Into(diags,X[i,i])
 )
);

¡Eso es mucho esfuerzo! Lo siento, me gustaría saber de una manera más fácil.

Re: Guardar opción de columnas con modelo mixto

Gracias por tu respuesta Cameron!

 

Sin embargo, si guardó el Std Err de los residuos, creo que está relacionado con aquellos que no incluyen el efecto aleatorio y no con los residuos condicionales. Parece que no es posible guardar el error estándar condicional de los residuos. Creo que podemos rehacer el cálculo en función de los valores de los sombreros, pero el problema es el mismo: probablemente estén vinculados a los "residuos básicos".

Re: Guardar opción de columnas con modelo mixto

Puedes estudiar los residuos condicionales dividiéndolos por sus errores estándar.  Guarde ambas columnas en la tabla y luego calcule los valores estudiantiles con una fórmula de columna.  No estoy seguro de cómo ayudarlo con la Fórmula de Pred.