Bom dia pessoal,
 
Estou trabalhando com a estimação de componentes de variância via modelos mistos usando o procedimento de máxima verossimilhança restrita (REML), no entanto, em minhas pesquisas pela internet encontrei o algoritmo de Fisher's Scoring
para solução do seguinte modelo:
 
y = Xa + Zg + e 
 
em que:
 
y é o vetor de dados
 
a e g são vetores de efeito fixo e aleatório respectivamente.
 
X e Z são as matrizes de incidência dos efeitos fixo e aleatório respectivamente.
 
 
Gostaria de poder adicionar mais efeitos aleatórios no modelo como por exemplo y = Xa + Zg + Wp + e, e encontrar as soluções dos componentes de variância pelo algoritmo de Fisher's Scoring.
 
Alguém me indica algum material onde possa encontrar algo a respeito?
 
Att.
 
Tiago.