Histórico: Nota: Meu conjunto de dados e código-r estão incluídos abaixo do texto
Desejo usar o AIC para comparar dois modelos de efeitos mistos gerados usando o pacote lme4 em R. Cada modelo tem um efeito fixo e um efeito aleatório. O efeito fixo difere entre os modelos, mas o efeito aleatório permanece o mesmo entre os modelos. Descobri que se eu usar REML = T, o modelo2 tem a pontuação mais baixa na AIC, mas se eu usar REML = F, o modelo1 tem a pontuação mais baixa na AIC.
Suporte para o uso de ML:
Zuur et al. (2009; PÁGINA 122) sugerem que "Para comparar modelos com efeitos fixos aninhados (mas com a mesma estrutura aleatória), a estimativa de ML deve ser usada e não a REML". Isso indica para mim que devo usar o ML, pois meus efeitos aleatórios são os mesmos nos dois modelos, mas meus efeitos fixos diferem. [Zuur et al. 2009. Modelos de efeitos mistos e extensões em ecologia com R. Springer.]
Suporte para usar REML:
No entanto, percebo que, quando uso ML, a variação residual associada aos efeitos aleatórios difere entre os dois modelos (modelo1 = 136,3; modelo2 = 112,9), mas quando uso REML, é a mesma entre os modelos (modelo1 = modelo2 = 151,5). Isso implica para mim que eu deveria usar REML para que a variação residual aleatória permaneça a mesma entre modelos com a mesma variável aleatória.
Questão:
Não faz mais sentido usar REML do que ML para comparações de modelos em que os efeitos fixos mudam e os efeitos aleatórios permanecem os mesmos? Caso contrário, você pode explicar por que ou me apontar para outra literatura que explica mais?
# Model2 "wins" if REML=T:
REMLmodel1 = lmer(Response ~ Fixed1 + (1|Random1),data,REML = T)
REMLmodel2 = lmer(Response ~ Fixed2 + (1|Random1),data,REML = T)
AIC(REMLmodel1,REMLmodel2)
summary(REMLmodel1)
summary(REMLmodel2)
# Model1 "wins" if REML=F:
MLmodel1 = lmer(Response ~ Fixed1 + (1|Random1),data,REML = F)
MLmodel2 = lmer(Response ~ Fixed2 + (1|Random1),data,REML = F)
AIC(MLmodel1,MLmodel2)
summary(MLmodel1)
summary(MLmodel2)
Conjunto de dados:
Response Fixed1 Fixed2 Random1
5.20 A A 1
32.50 A A 1
6.57 A A 2
24.77 A B 3
41.69 A B 3
34.29 A B 4
1.80 A B 4
10.00 A B 5
15.56 A B 5
4.44 A C 6
21.65 A C 6
9.20 A C 7
4.11 A C 7
12.52 B D 8
0.25 B D 8
27.34 B D 9
11.54 B E 10
0.86 B E 10
0.68 B E 11
4.00 B E 11