<div dir="ltr">Veja se na Task View de farmacocinética algum dos pacotes indicados não lhe serve.<div><br></div><div>HTH</div><div>--</div><div>Cesar Rabak</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-12-20 22:26 GMT-02:00 Maurício Lordêlo via R-br <span dir="ltr"><<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Caros,</div><div>Foi realizado um experimento no qual deseja-se saber qual composto extraído de uma determinada espécie vegetal é mais eficiente para matar uma determinado tipo de inseto.</div><div>Foram utilizadas 32 placas cada uma contendo 10 insetos vivos.</div><div>Em cada placa, foram realizadas 4 repetições das 8 diferentes doses do composto de cada espécie.</div><div>O número de insetos mortos foi registrado após um determinado período de tempo.</div><div>O objetivo é encontrar qual a dose, para cada um das três espécies, que mata 50% dos insetos.</div><div>Alguma sugestão de como encontrar os modelos para as três espécies? </div><div><br></div><div><br></div><div><br></div><div>dose = rep(c(0,0.15625,0.31250,0.<wbr>62500,1.2500, 2.5000, 5.0000, 10.000), each=4)</div><div>n_insetos = rep(10,32)</div><div>m_especie1 = c(1,3,4,0,5,2,5,5,4,4,2,2,0,3,<wbr>3,5,10,10,7,0,10,10,10,10,10,<wbr>10,10,10,10,10,10,10)</div><div>m_especie2 = c(0,3,4,4,1,0,4,0,0,1,0,0,3,2,<wbr>0,3,3,2,0,3,3,3,3,6,9,4,7,2,<wbr>10, 10, 10, 10)</div><div>m_especie3 = c(4,2,2,3,4,6,8,6,4,6,6,5,5,9,<wbr>5,8,10,8,5,5,10,10,10,10,10,<wbr>10,10,10,10,10,10,10)</div><div>pm_especie1 = m_especie1/n_insetos</div><div>pm_especie2 = m_especie2/n_insetos</div><div>pm_especie3 = m_especie3/n_insetos</div><div><br></div><div>dose.fator = as.factor(dose)</div><div><br></div><div>boxplot(pm_especie1~dose.<wbr>fator, xlab= "Dose", ylab="Proporção de mortos espécie 1",cex.axis=.85)</div><div>boxplot(pm_especie2~dose.<wbr>fator, xlab= "Dose", ylab="Proporção de mortos espécie 2",cex.axis=.85)</div><div>boxplot(pm_especie3~dose.<wbr>fator, xlab= "Dose", ylab="Proporção de mortos espécie 3",cex.axis=.85)</div><div><br></div><div>bartlett.test(pm_especie1~<wbr>dose.fator)</div><div>bartlett.test(pm_especie2~<wbr>dose.fator)</div><div>bartlett.test(pm_especie3~<wbr>dose.fator)</div><div><br></div><div>modelo1 = glm(cbind(good=m_especie1, bad=n_insetos-m_especie1)~<wbr>dose,family="binomial", data=dados)</div><div>par(mfrow=c(2,2))<br></div><div>plot(modelo1)</div><div>shapiro.test(modelo1$<wbr>residuals)</div><div><br></div><div>modelo2 = glm(cbind(good=m_especie2, bad=n_insetos-m_especie2)~<wbr>dose,family="binomial", data=dados)</div><div><div>par(mfrow=c(2,2))<br></div><div>plot(modelo2)</div></div><div>shapiro.test(modelo2$<wbr>residuals)</div><div><br></div><div>modelo3 = glm(cbind(good=m_especie3, bad=n_insetos-m_especie3)~<wbr>dose,family="binomial", data=dados)</div><div><div>par(mfrow=c(2,2))<br></div><div>plot(modelo2)</div></div><div>shapiro.test(modelo3$<wbr>residuals)</div><div><br></div></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/<wbr>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-<wbr>guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>