Fábio,<div>Não tem bloco dentro do local, foi em DIC, apenas uma repetição dentro de cada local.</div><div>Att<br><br><div class="gmail_quote">Em 22 de outubro de 2012 07:05, Fabio Mathias Corrêa <span dir="ltr"><<a href="mailto:fabio.ufla@yahoo.com.br" target="_blank">fabio.ufla@yahoo.com.br</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-size:12pt;font-family:times new roman,new york,times,serif"><div><span>Não tem bloco dentro de local?</span></div>
<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif"><span><br></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif">
<span><br></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif"><span>Foi inteiramente casualizado?</span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif">
<span><br></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif">Valeu!</div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'times new roman','new york',times,serif">
<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.: <a href="tel:73-3680-5076" value="+557336805076" target="_blank">73-3680-5076</a><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 face="Arial"> <hr size="1"> <b><span style="font-weight:bold">De:</span></b> Walmes Zeviani <<a href="mailto:walmeszeviani@gmail.com" target="_blank">walmeszeviani@gmail.com</a>><br>
<b><span style="font-weight:bold">Para:</span></b> <a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a> <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><div><div class="h5"> <br><div><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><br></div><div>À disposição.</div><div>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: <a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573" target="_blank">(+55) 41 3361 3573</a></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" href="mailto:walmes@ufpr.br" target="_blank">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" href="http://www.leg.ufpr.br/~walmes" target="_blank">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></div></div><div class="im">_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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></div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a 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></blockquote></div><br></div>