19 Mar
2015
19 Mar
'15
20:49
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)