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<-c(rnorm(20,2),rnorm(10,5),rnorm(20,5),rnorm(10,8))</div><div>Perigo<-rep(c("Baixo","Alto"),each=30)</div><div>Id<-c(rep(1:10,3),rep(11:20,3))</div>
<div>Dia<-rep(rep(1:3,each=10),2)</div><div>dados<-data.frame(cbind(Perigo,Dia,Id,Resposta))</div><div>dados$Resposta<-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="black", cex=1.3,</div><div>       strip=strip.custom(bg="white"),type=c("p","a"),scale=list(tck=c(1,0),alternating=1))</div>
<div><br></div><div>#Anova de medidas repetidas</div><div>modelo<-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>