[R-br] metodologia para validação de tratamentos

[Felipe]

Talvez o pacote multcomp possa te ajudar:

http://cran.r-project.org/web/packages/multcomp/multcomp.pdf

Exemplo:
library(multcomp)
Contr=c(123,179,121,185,183,191,182,187)
A=c(96,79,91,85,83,91,82,87)
B=c(77,76,74,73,78,71,80)
C=c(66,73,69,66,77,73,71,70,74)

  resp=c(Contr,A,B,C)
  Metodo=c(rep(1,length(Contr)), rep(2,length(A)), rep(3,length(B)), 
rep(4,length(C)))
  Metodo=factor(Metodo)
  anova=aov(resp~Metodo)
  summary(anova)

test.dunnett=glht(anova,linfct=mcp(Metodo="Dunnett"))
confint(test.dunnett)
plot(test.dunnett)