Ejemplo de Ajustar modelo de efectos mixtos

Un investigador evalúa el rendimiento de seis variedades de alfalfa en cuatro campos seleccionados de forma aleatoria. El rendimiento de cada variedad se registró para cada campo.

El investigador desea saber si la variedad de alfalfa afecta el rendimiento medio. El investigador tiene 4 campos donde puede recopilar datos. Sin embargo, quiere poder modelar cómo crecerán las alfalfas en campos que no están incluidos en el experimento. Por lo tanto, el investigador hace que el campo donde crece la alfalfa sea un factor aleatorio. El investigador utiliza un modelo de efectos mixtos para evaluar los efectos fijos y aleatorios de forma conjunta.

  1. Abra los datos de muestra Alfalfa.MTW.
  2. Elija Estadísticas > ANOVA > Modelo de efectos mixtos > Ajustar modelo de efectos mixtos.
  3. En Respuestas, ingrese Rendimiento.
  4. En Factores aleatorios (obligatorio), ingrese Campo.
  5. En Factores fijos, ingrese Variedad.
  6. Haga clic en Gráficas.
  7. En Residuos para las gráficas, seleccione Condicionales estandarizados.
  8. En Gráficas de residuos, seleccione Cuatro en una.
  9. Haga clic en Aceptar en cada cuadro de diálogo.

Interpretar los resultados

En la tabla Componentes de la varianza, el valor p de Campo es 0.124. La prueba de hipótesis no muestra evidencia de que el componente de la varianza sea diferente de 0. El valor p para el componente de la varianza del error es 0.003. Puesto que el valor p es menor que el nivel de significancia de 0.05, el investigador puede concluir que el componente de la varianza del error no es 0.

El valor p de aproximadamente 0 para el término de factor fijo, Variedad, muestra que al menos el efecto de un tipo de alfalfa sobre el rendimiento es significativamente diferente de los otros cinco tipos.

Los coeficientes de los efectos principales representan la diferencia entre cada media de nivel y la media general. Por ejemplo, la Variedad 1 está asociada a un rendimiento de alfalfa que es aproximadamente de 0.385 unidades mayor que la media general. El valor p de aproximadamente 0 de este coeficiente indica que el efecto de la Variedad 1 en el rendimiento es significativamente diferente de otro efecto de nivel del término Variedad. Para determinar qué efectos de nivel son estadísticamente iguales y qué efectos de nivel son estadísticamente diferentes, el investigador planea hacer un análisis de comparación múltiple para el término.

El valor de R2 muestra que el modelo explica aproximadamente el 92% de la variación en el rendimiento. El valor de R-cuad. (ajust.) también es alto, con un valor de aproximadamente 90.2%. El investigador utiliza este valor para comparar modelos que tienen diferente número de predictores.

Las observaciones 1 y 5 son observaciones poco comunes, porque tienen residuos estandarizados mayores que 2. El investigador examina los datos para asegurarse de que los valores de respuesta de esas observaciones sean correctos.

Los residuos en la gráfica de probabilidad normal se aproximan a una línea recta, mientras que los puntos parecen estar dispersos alrededor del 0 en la gráfica de residuos vs. ajustes.

Método

Estimación de la varianzaMáxima verosimilitud restringida
GL para efectos fijosKenward-Roger

Información del factor

FactorTipoNivelesValores
CampoAleatorio41, 2, 3, 4
VariedadFijo61, 2, 3, 4, 5, 6

Componentes de la varianza

FuenteVar% del totalEE de la
var.
Valor ZValor p
Campo0.07791972.93%0.0675801.1529960.124
Error0.02892427.07%0.0105622.7386130.003
Total0.106843       
-2 log likelihood = 7.736012

Pruebas de efectos fijos

TérminoGL NumGL DenValor FValor p
Variedad5.0015.0026.290.000

Resumen del modelo

SR-cuadradoR-cuadrado(ajustado)AICcBIC
0.17007192.33%90.20%12.5413.52

Coeficientes

TérminoCoefEE del coef.GLValor TValor p
Constante3.0945830.1438223.0021.5166920.000
Variedad         
  10.3854170.07762615.004.9650160.000
  20.1454170.07762615.001.8732870.081
  30.1079170.07762615.001.3902050.185
  4-0.3195830.07762615.00-4.1169380.001
  50.3954170.07762615.005.0938380.000

Ajustes y diagnósticos marginales para observaciones poco comunes

ObsRendimientoAjusteResidResid est.
14.1000003.4800000.6200002.190221R
54.2200003.4900000.7300002.578808R
Residuo grande R

Ajustes y diagnósticos condicionales para observaciones poco comunes

ObsRendimientoAjusteResidResid est.
54.2200003.8953390.3246612.400733R
Residuo grande R