<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><span>Éder,</span></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><span><br></span></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; ">Os efeitos aleatórios devem ser especificados como no exemplo abaixo.</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><br></div><div>m1 <- lmer(diametro~1+(1|M1) + (1|M2) + (1|M3) + (1|M4) + (1|M5) + (1|M6) + (1|M7),dados)<br></div><div>summary(m1)</div><div>ranef(m1)</div><div><br></div><div><div>dotplot(ranef(m1, postVar=TRUE))$M1 # Plota os efeitos aleatórios e seus IC's.</div><div><br></div><div>Valeu!!!</div><div><br></div></div><div style="font-family: 'times new
roman', 'new york', times, serif; font-size: 12pt; "> </div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "> Fábio Mathias Corrêa<br><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "> Universidade Estadual de Santa Cruz<br>Departamento de Ciências Exatas e da Terra - DCET</div><br><br><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; ">Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna<br>CEP 45662-900. Ilhéus-Bahia</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "><br><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; ">Tel.: 73-3680-5076<br><div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "><div style="font-size: 12pt; font-family: 'times
new roman', 'new york', times, serif; "><font size="2" face="Arial"><hr size="1"><b><span style="font-weight:bold;">De:</span></b> Eder David Borges da Silva <eder@leg.ufpr.br><br><b><span style="font-weight: bold;">Para:</span></b> R-br@listas.c3sl.ufpr.br<br><b><span style="font-weight: bold;">Enviadas:</span></b> Domingo, 21 de Agosto de 2011 15:52<br><b><span style="font-weight: bold;">Assunto:</span></b> [R-br] Declarar matriz Z modelos mistos<br></font><br><div id="yiv1552882754">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 rel="nofollow" target="_blank" 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>
</div><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br></div></div></div></div></body></html>