Exemplo de Ajustar modelo de efeitos mistos

Um pesquisador testa a produção de seis variedades de alfafa em quatro campos selecionados aleatoriamente. A produção de cada variedade foi registrada para cada campo.

O pesquisador quer saber se a variedade de alfafa afeta o rendimento médio. O pesquisador tem 4 campos onde eles podem coletar dados. No entanto, o pesquisador quer ser capaz de modelar como as alfafas crescerão nos campos que não estão no experimento. Assim, o pesquisador faz com que o campo onde a alfafa cresce um fator aleatório. O pesquisador usa um modelo de efeitos mistos para avaliar os efeitos fixos e aleatórios juntos.

  1. Abra os dados das amostras Alfafa.MTW.
  2. Selecione Estat > ANOVA > Modelo de efeitos mistos > Ajustar modelo de efeitos mistos.
  3. Em Respostas, insira Rendimento.
  4. Em Fatores aleatórios (obrigatórios), insira Campo.
  5. Em Fatores fixos, insira Variedade.
  6. Clique em Gráficos.
  7. Em Resíduos para gráficos, selecione Padronizado condicional.
  8. Em Gráficos de resíduos, selecione Quatro em um.
  9. Clique em OK em cada caixa de diálogo.

Interpretar os resultados

Na tabela Componentes de Variância, o valor-p para Campo é 0,124. O teste de hipótese não mostra evidências de que o componente de variância é diferente de 0. O valor-p para o componente de variância de erro é de 0,003. Como o valor-p é menor do que o nível de significância de 0,05, o pesquisador pode concluir que o componente de variância para o erro não é 0.

O valor-p de aproximadamente 0 para o termo de fator fixo, Variedade, mostra que pelo menos um tipo de efeito sobre a produção de alfafa é significativamente diferente dos outros cinco tipos.

Os coeficientes para os efeitos principais representam a diferença entre cada nível médio e a média geral. Por exemplo, Variedade 1 está associada com um rendimento de alfafa que é de aproximadamente 0,385 unidades maior do que a média global. O valor-p de aproximadamente 0 para este coeficiente indica que o efeito da Variedade 1 sobre o rendimento é significativamente diferente de um outro efeito do nível do termo de Variedade. Para determinar quais os efeitos de nível são estatisticamente iguais, e quais efeitos de nível são estatisticamente diferentes, o pesquisador pretende fazer uma análise de comparação múltipla para o termo.

O valor R2 mostra que o modelo explica cerca de 92% da variação no rendimento. O valor de R-quadrado (adj) também é elevado, com um valor de aproximadamente 90,2%. O pesquisador usa esse valor para comparar os modelos que têm números diferentes de preditores.

As observações 1 e 5 são observações incomuns porque elas têm resíduos padronizados superiores a 2. O pesquisador examina os dados para se certificar de que os valores de resposta para essas observações estão corretos.

Os resíduos no gráfico de probabilidade normal se aproximam de uma reta, e os pontos parecem estar dispersos aleatoriamente em torno de 0 no gráfico de resíduos versus ajustes.

Método

Estimação de variânciaMáxima verossimilhança restrita
GL para efeitos fixosKenward-Roger

Informações dos Fatores

FatorTipoNíveisValores
CampoAleatório41; 2; 3; 4
VariedadeFixo61; 2; 3; 4; 5; 6

Componentes de Variância

FonteVar% do TotalEP de VarValor-ZValor-P
Campo0,07791972,93%0,0675801,1529960,124
Erro0,02892427,07%0,0105622,7386130,003
Total0,106843       
-2 Log-verossimilhança = 7,736012

Testes de Efeitos Fixos

TermoNum GLDen GLValor FValor-P
Variedade5,0015,0026,290,000

Sumário do Modelo

SR2R2(aj)AICcBIC
0,17007192,33%90,20%12,5413,52

Coeficientes

TermoCoefEP de CoefGLValor-TValor-P
Constante3,0945830,1438223,0021,5166920,000
Variedade         
  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

Valores Ajustados e Diagnósticos Marginais para Observações Atípicas

Obs.RendimentoAjusteResídResíd Pad
14,1000003,4800000,6200002,190221R
54,2200003,4900000,7300002,578808R
R  Resíduo grande

Valores Ajustados e Diagnósticos Condicionais para Observações Atípicas

Obs.RendimentoAjusteResídResíd Pad
54,2200003,8953390,3246612,400733R
R  Resíduo grande