<div dir="ltr">Só para testar meu entendimento, as pupas e as larvas são objetos <i>diferentes</i>, então formalmente poder-se-ia dizer que foram dois experimentos feitos "em paralelo", correto?<div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-02-10 20:06 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, <br>Foi realizado um experimento para avaliar o efeito de oito <br>diferentes tratamentos na mortalidade de insetos (foram duas variáveis respostas).<br>O pesquisador colocou em cada unidade experimental 10 larvas<br>e registrou o número de larvas mortas e de pupas inviáveis<br>em quatro diferentes tempos sendo cinco repetições.<br>Este registro no tempo não foi feito de forma independente e sim de forma acumulada.<br>Ou seja, se o número de mortos no tempo t1 é igual 2 e no tempo t2<br>é igual a 5, indica entre t1 e t2 "morreram" 3 insetos (na mesma unidade experimental) .<br>O objetivo é avaliar o efeito do tratamento e do tempo.<br>Como não há independência do tempo, alguém tem alguma sugestão de como analisar estes dados?<br><br>rm(list = ls(all=TRUE))<br></div>#Tratamentos<br><div>trat = as.factor(c(rep("Trat1",20),re<wbr>p("Trat2",20), rep("Trat3",20),rep("Trat4",20<wbr>),rep("Trat5",20), rep("Trat6",20),rep("Trat7",20<wbr>),rep("Trat8",20)))<br><br></div><div>#Tempo<br></div><div>tempo = as.factor(rep(c("t1","t2","t3"<wbr>,"t4"), each = 5, times =8))<br><br>#Número total de larvas em cada unidade experimental<br>n_total_larvas= rep(10,160) <br><br>#Número de larvas mortas<br>n_larvas_mortas = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 3, 0, 0, 1, 1, 3, 0, 3, 0, 0, 0, 0, 3, 1, 1, <br> 1, 1, 5, 3, 1, 2, 2, 5, 3, 3, 4, 3, 0, 1, 0, 0, 0, 0, 2, 1, 2, 1, 0, 4, 1, 2, 1, 4,<br> 4, 2, 2, 2, 0, 0, 1, 1, 0, 0, 0, 1, 2, 0, 2, 2, 1, 2, 0, 2, 3, 2, 2, 4, 0, 0, 0, 0,<br> 0, 0, 0, 1, 0, 0, 0, 2, 1, 1, 0, 0, 2, 1, 1, 0, 0, 1, 1, 1, 1, 0, 3, 1, 2, 1, 0, 3,<br> 1, 2, 2, 1, 3, 1, 3, 2, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 2, 0, 1, 0, 1, 3, 0, 1,<br> 8, 8, 7, 9, 8, 8, 8, 7, 9, 8, 8, 9, 7, 9, 8, 8, 9, 7, 9, 8)<br>#Número de pupas inviáveis<br>n_pupas_inv = c(3, 4, 3, 6, 1, 3, 8, 6, 6, 1, 10, 9, 9, 4, 6, 10, 9, 9, 7, 10, 0, 0, <br> 0, 0, 0, 0, 0, 0, 1, 0, 5, 7, 4, 6, 5, 5, 7, 7, 6, 7, 2, 1, 0, 2,<br> 0, 2, 4, 2, 5, 1, 4, 4, 2, 6, 7, 6, 6, 8, 6, 8, 0, 3, 0, 2, 0, 4,<br> 5, 0, 4, 3, 7, 6, 5, 5, 5, 8, 7, 7, 8, 5, 1, 0, 2, 6, 0, 3, 3, 2,<br> 9, 9, 6, 7, 7, 9, 9, 10, 8, 9, 9, 10, 0, 1, 0, 0, 0, 5, 4, 3, 2, 3,<br> 9, 7, 8, 4, 7, 9, 7, 9, 6, 7, 1, 2, 0, 0, 1, 6, 8, 2, 2, 1, 10, 8,<br> 6, 9, 6, 10, 9, 7, 10, 8, 2, 1, 1, 0, 0, 2, 1, 3, 1, 1, 2, 1, 3, 1,<br> 2, 2, 1, 3, 1, 2)<br><br>p_larvas_mortas = n_larvas_mortas/n_total_larvas<br>table(p_larvas_mortas)<br>p_pupas_inv = n_pupas_inv/n_total_larvas<br>table(p_pupas_inv)<br>dados_larvas = data.frame(trat, tempo, n_total_larvas, n_larvas_mortas,n_pupas_inv,p_<wbr>larvas_mortas,p_pupas_inv)<br>str(dados_larvas)<div class="m_5061034050186456644gmail-yj6qo"></div><div class="m_5061034050186456644gmail-adL">tapply(p_larvas_mortas, list(trat, tempo), mean, na.rm=T)<br></div></div>
<br></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>