Ola, pessoal eu estou apanhando aqui para conseguir realizar uma analise de contrastes pra poder dizer que os resultados do Dia 1 são iguais do Dia 2, ai so no Dia 3 começa a aumentar.<div>Eu vi um exemplo usando a função gls do pacote nlme, mas não consegui reproduzir pra conseguir verificar o que falo acima.</div>
<div>Se alguém puder indicar um tutorial ou algum texto não muito complicado sobre como avaliar contrastes para anova de medidas repetidas, muito me ajudaria. </div><div><br></div><div><div>Segue um exemplo dos dados:</div>
<div><br></div><div><div>#Criandos dados</div><div>set.seed(1)</div><div>Resposta&lt;-c(rnorm(20,2),rnorm(10,5),rnorm(20,5),rnorm(10,8))</div><div>Perigo&lt;-rep(c(&quot;Baixo&quot;,&quot;Alto&quot;),each=30)</div><div>Id&lt;-c(rep(1:10,3),rep(11:20,3))</div>
<div>Dia&lt;-rep(rep(1:3,each=10),2)</div><div>dados&lt;-data.frame(cbind(Perigo,Dia,Id,Resposta))</div><div>dados$Resposta&lt;-as.numeric(levels(dados$Resposta))[dados$Resposta]</div><div><br></div><div>#Olhando os dados com um grafico</div>
<div>str(dados)</div><div>library(lattice)</div><div>xyplot(Resposta~Dia|Perigo,data=dados, pch=16, col=&quot;black&quot;, cex=1.3,</div><div>       strip=strip.custom(bg=&quot;white&quot;),type=c(&quot;p&quot;,&quot;a&quot;),scale=list(tck=c(1,0),alternating=1))</div>
<div><br></div><div>#Anova de medidas repetidas</div><div>modelo&lt;-aov(Resposta~Perigo*Dia+Error(Id),data=dados)</div><div>summary(modelo)</div></div><div><br>-- <br><div>Grato<br>Augusto C. A. Ribas</div>
<div> </div>
<div>Site Pessoal: <a href="http://augustoribas.heliohost.org" target="_blank">http://augustoribas.heliohost.org</a></div>
<div>Lattes: <a href="http://lattes.cnpq.br/7355685961127056" target="_blank">http://lattes.cnpq.br/7355685961127056</a><br></div><br>
</div></div>