<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:10pt"><div>glmer1 <- glmer(cbind(morto, total - morto) ~ tratamento + tempo + tratamento:tempo + (1 | placapetri), family = binomial, data = dados))</div><div><br></div><div style="color: rgb(0, 0, 0); font-size: 13px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; background-color: transparent; font-style: normal;">Use glmer em lme4 + step para anova em lmeTest</div><div style="color: rgb(0, 0, 0); font-size: 13px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; background-color: transparent; font-style: normal;"><br></div><div style="color: rgb(0, 0, 0); font-size: 13px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; background-color: transparent; font-style:
normal;">att.</div><div> </div><div>Alisson Lucrécio da Costa</div><div class="yahoo_quoted" style="display: block;"> <br> <br> <div style="font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-size: 10pt;"> <div style="font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-size: 12pt;"> <div dir="ltr"> <font size="2" face="Arial"> On Monday, January 20, 2014 9:50 AM, walmes . <walmeszeviani@gmail.com> wrote:<br> </font> </div> <div class="y_msg_container"><div id="yiv8085261503"><div><div dir="ltr"><div class="yiv8085261503gmail_default" style="font-family: 'trebuchet ms', sans-serif;">Vejamos se entendi corretamente. A placa de petri é sua unidade experimental que recebe um nível do fator (contínuo) ração que tem 6 níveis. A placa é subdividida e 4 partes e cada uma delas é observada respeitando níveis do segundo fator (contínuo) tempo. Se
são 4 partes então são 4 níveis de tempo. A resposta avaliada é o número de mortos e número de vivos. Sob certas suposições, a probabilidade de mortos pode ser associada à distribuição binomial (x=número de mortos, n=número total=mortos+vivos). Ainda pode-se considerar o efeito aleatório da unidade experimental (ue) pois trata-se de um experimento de parcela subdividida, muito embora não seja de medida repetida por existe uma fatia na placa para cada tempo. Então o modelo candidato seria<br clear="none">
<br clear="none"></div><div class="yiv8085261503gmail_default" style="font-family: 'trebuchet ms', sans-serif;">y ~ binomial(p=f(...), n=n)<br clear="none"></div><div class="yiv8085261503gmail_default" style="font-family: 'trebuchet ms', sans-serif;">f(...) = ração+tempo+ração:tempo+ue<br clear="none">
</div><div class="yiv8085261503gmail_default" style="font-family: 'trebuchet ms', sans-serif;">ue ~ normal(0, sigma_u)<br clear="none"><br clear="none"></div><div class="yiv8085261503gmail_default" style="font-family: 'trebuchet ms', sans-serif;">Se o n for grande e x/n não estiver tão nas bordas (perto de zero ou um), então não é tão contraindicado usar a distribuição normal no lugar da binomial. Lembrando que é possível aplicar uma transformação, como logito, gompito, probito, para "ter mais normalidade". Mesmo diante dessas alternativas, eu considero a abordagem glmm (binomial com efeitos aleatórios) mais apropriada. Verifique a função lme4::glmer().<br clear="none">
<br clear="none">À disposição.<br clear="none">Walmes.<div class="yiv8085261503yqt4385259769" id="yiv8085261503yqtfd18836"><br clear="none"></div></div></div></div></div><br><div class="yqt4385259769" id="yqtfd35728">_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a shape="rect" 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 clear="none">Leia o guia de postagem (<a shape="rect" 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.</div><br><br></div> </div> </div> </div> </div></body></html>