Teste de Tukey para comparações múltiplas

No R, o teste de Tukey é apresentado através de intervalos de confiança. A interpretação é: se o intervalo de confiança para a diferença entre duas médias não incluir o valor zero, siginifica que rejeita-se a hipótese nula, caso contrário, não rejeita-se. O resultado pode ser visto através de uma tabela e/ou graficamente:

ex01.tu <- TukeyHSD(ex01.av)
ex01.tu
plot(ex01.tu)

Uma outra maneira é utilizar a seguinte função para estes dados.

## Diferença entre médias para um fator e igual número de repetições
#r=número de repetições
#t=número de tratamentos

dif.medias<-function(dados=ex01,r=6, t=9,alpha=0.95 )
{ attach(dados)
  modelo<-aov(resp~trat,data=dados)
  trat.m<-tapply(resp,trat,mean) 
  trat.m1<-trat.m 
  m1d<-outer(trat.m1,trat.m1,"-") 
  m1d<-m1d[lower.tri(m1d)] 
  m1n<-outer(names(trat.m1),names(trat.m1),paste,sep="-") 
  names(m1d)<-m1n[lower.tri(m1n)] 
  s2<-sum(resid(modelo)^2)/modelo$df.res 
  n<-r 
  dif.t<-qtukey(alpha,t,modelo$df.res)*sqrt(s2/n) 
  data.frame(dif=m1d,sig=ifelse(abs(m1d)>dif.t,"*","ns")) 
}