[R-br] Resíduos e MCMCglmm
Elias Teixeira Krainski
eliaskrainski em yahoo.com.br
Segunda Novembro 23 13:22:35 BRST 2015
Estou confuso agora se você quer fazer análise Bayesiana ou não...
Primeiro: estou pensando em modelos hierárquicos enquanto deduzo que o
Wagner está pensando em modelos marginais.
Segundo: Resíduo (se houver esse tipo de coisa no modelo) segue
distribuição Normal se você assumiu isso para comecar.
Esse termo (resíduo) é comum em modelos do tipo
Y = X beta + Z b + e
X: covariáveis/desenho/features
beta: coeficientes de regressão
Z: matriz dos efeitos aleatórios
b: efeitos aleatórios
r: os tais resíduos,
Matematicamente, coeficientes de regressão e efeitos aleatórios podem
ser tratados juntos (basta juntar X e Z). Geralmente, temos que b ~ N(0,
alguma.matriz), r ~ N(0, sigma^2 * I). Ou seja, aqui assume-se
Normalidade para r. Neste modelo, ambos, r|y e b|y, também têm
distribuição Normal.
Se invés do modelo acima você tiver
E(Y) = eta = X beta + Z b
Y|eta ~ Distribuição.não.Gaussiana(...)
Ou seja, não existe resíduos Normais. A distribuição Normal é
(geralmente) ainda assumida para os efeitos aleatórios, mas,
condicionado aos dados isso não é mais o caso, isto é, b|y não segue
distribuição Normal.
Este último modelo recai no anterior se Y|eta ~ N(0, sigma^2 * I).
Assim, apenas quando a verossimilhança é Normal temos b|y e r|y com
distribuição Normal.
Exite a possibilidade de se adicionar um termo de "resíduo" considerando
E(Y) = eta = X beta + Z b + r
em situações em que y|eta pode ter uma variação adicional que não torne
o modelo não-identificável. Nestes casos, mesmo que r ~ Normal(0,
sigma^2 * I), r|y não é Normal.
Elias
Mais detalhes sobre a lista de discussão R-br