<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Não tem bloco dentro de local?</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><span>Foi inteiramente casualizado?</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><span><br></span></div><div style="color: rgb(0,
0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; ">Valeu!</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal; "><span><br></span></div><div> </div><div> Fábio Mathias Corrêa<br><br></div><div> Universidade Estadual de Santa Cruz<br>Departamento de Ciências Exatas e da Terra - DCET</div><br><br><div>Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna<br>CEP 45662-900. Ilhéus-Bahia</div><div><br><br></div><div>Tel.: 73-3680-5076<br> <div style="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; "> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span
style="font-weight:bold;">De:</span></b> Walmes Zeviani <walmeszeviani@gmail.com><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 Outubro de 2012 16:30<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] Blocos com algumas repetições do material<br> </font> </div> <br><div id="yiv1479791909"><div><font face="trebuchet ms,sans-serif">Éder,</font></div><div><font face="trebuchet ms,sans-serif"><br></font></div><div><font face="trebuchet ms,sans-serif">Agora entendi suas colocações. Bem, dessa maneira acredito que o modelo deva ser declarado dessa maneira.</font></div>
<div><font face="trebuchet ms,sans-serif"><br></font></div><div><div><font face="courier new, monospace">#------------------------------------------------------------------------------------------</font></div><div><font face="courier new, monospace"># modelo com efeito aleatório e local, cultivar, local:cultivar e resídual</font></div>
<div><font face="courier new, monospace"># só é possível estimar resídual se haver repetições dentro de locais:cultivar</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"># número de níveis</font></div>
<div><font face="courier new, monospace">nl <- 10; nc <- 8; nr <- 3</font></div><div><font face="courier new, monospace"># efeitos aleatórios</font></div><div><font face="courier new, monospace">l <- rnorm(nl,0,2); # de local</font></div>
<div><font face="courier new, monospace">c <- rnorm(nc,0,1.5); # de cultivar</font></div><div><font face="courier new, monospace">lc <- rnorm(nl*nc,0,1); # de local:cultivar (interação genótipo:ambiente)</font></div>
<div><font face="courier new, monospace">e <- rnorm(nl*(nc+nr),0,0.5) # residual, var entre rep no mesmo local:cult</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"># dados artificiais</font></div>
<div><font face="courier new, monospace">da1 <- expand.grid(local=gl(nl,1), cultivar=gl(nc,1), r=1)</font></div><div><font face="courier new, monospace">da2 <- expand.grid(local=gl(nl,1), cultivar=gl(nr,1), r=2)</font></div>
<div><font face="courier new, monospace">da <- rbind(da1, da2)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"># matriz de delineamento</font></div><div><font face="courier new, monospace">Z <- with(da, cbind(model.matrix(~-1+local),</font></div>
<div><font face="courier new, monospace"> model.matrix(~-1+cultivar),</font></div><div><font face="courier new, monospace"> model.matrix(~-1+local:cultivar)))</font></div><div><font face="courier new, monospace">dim(Z)</font></div>
<div><font face="courier new, monospace">ranef <- c(l, c, lc)</font></div><div><font face="courier new, monospace">length(ranef)</font></div><div><font face="courier new, monospace">y <- Z%*%ranef+e</font></div><div>
<font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">require(lme4)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">m0 <- lmer(y~(1|local)+(1|cultivar)+(1|local:cultivar), data=da)</font></div>
<div><font face="courier new, monospace">summary(m0)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">#------------------------------------------------------------------------------------------</font></div>
<div style=""><br></div><div style="">À disposição.</div><div style="">Walmes.</div></div>
<font face="trebuchet ms,sans-serif"><br clear="all"></font><span style="font-family: 'trebuchet ms', sans-serif; ">==========================================================================</span><br style="font-family: 'trebuchet ms', sans-serif; ">
<span style="font-family: 'trebuchet ms', sans-serif; ">Walmes Marques Zeviani</span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family: 'trebuchet ms', sans-serif; ">
<span style="font-family: 'trebuchet ms', sans-serif; ">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">fone: (+55) 41 3361 3573</span><br style="font-family: 'trebuchet ms', sans-serif; ">
<span style="font-family: 'trebuchet ms', sans-serif; ">VoIP: (3361 3600) 1053 1173</span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">e-mail: <a rel="nofollow" ymailto="mailto:walmes@ufpr.br" target="_blank" href="mailto:walmes@ufpr.br">walmes@ufpr.br</a><br>
skype: walmeszeviani<br style="font-family: 'trebuchet ms', sans-serif; "></span><span style="font-family: 'trebuchet ms', sans-serif; ">twitter: @walmeszeviani</span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">homepage: <a rel="nofollow" target="_blank" href="http://www.leg.ufpr.br/~walmes">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family: 'trebuchet ms', sans-serif; ">
<span style="font-family: 'trebuchet ms', sans-serif; ">linux user number: 531218</span><br style="font-family: 'trebuchet ms', sans-serif; "><span style="font-family: 'trebuchet ms', sans-serif; ">==========================================================================</span><br>
</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>