Pessoal,<div>dado um modelo y = Xb + Zu + e , modelo misto com Xb parte fixa e Zu parte aleatoria, estou tentando entender melhor a analise feita neste documento, <a href="http://www.infoteca.cnptia.embrapa.br/bitstream/doc/883425/1/Doc210.pdf">http://www.infoteca.cnptia.embrapa.br/bitstream/doc/883425/1/Doc210.pdf</a> , pagina 54, ao meu entendimento é um modelo misto, onde a mudança é a forma que é construída a matriz Z, que como o autor faz é produzido um efeito aleatório para cada marcador. desta forma gostaria de saber se alguém ja declarou um matriz Z em uma analise de modelo misto no R em especifico na lmer, essa analise é GWS (Seleção genômica ampla), caso alguém tenha familiaridade ficaria grato de umas dicas.</div>
<div><br></div><div><div>dados <- data.frame(ind=c(1:5),</div><div>                    diametro=c(9.87,14.48,8.91,14.64,9.55),</div><div>                    M1=c(2,1,0,1,1),</div><div>                    M2=c(0,1,2,0,0),</div>
<div>                    M3=c(0,0,0,1,0),</div><div>                    M4=c(0,0,0,0,1),</div><div>                    M5=c(2,1,0,1,1),</div><div>                    M6=c(0,1,0,0,1),</div><div>                    M7=c(0,0,2,0,0))</div>
<div>dados</div><div>require(lme4)</div><div><br></div><div>Z <- model.matrix(~-1+M1+M2+M3+M4+M5+M6+M7,dados)</div><div>Z ### Esta é a matriz que querro por no modelo</div><div>m0 <- lmer(diametro~1+(1|?????),dados)</div>
<div>summary(m0)</div><div>ranef(m0)</div><div>fixef(m0)</div></div><div># Meu interesse e ter os componentes de variância e o efeito fixo(Media geral), e os efeitos aleatórios (7, uma para cada marcador)</div><div><br></div>
<div>OBS: Creio que não chegaremos ao valor exato do autor, visto que ele fixou (sigma.erro/(sigma.g/n)) =1.</div><div>Obrigado</div>