<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),<wbr>rep("Trat2",20), rep("Trat3",20),rep("Trat4",<wbr>20),rep("Trat5",20), rep("Trat6",20),rep("Trat7",<wbr>20),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="gmail-yj6qo"></div><div class="gmail-adL">tapply(p_larvas_mortas, list(trat, tempo), mean, na.rm=T)<br></div></div>

<br></div>