Comparación de medias

Una vez que se ha determinado que existen diferencias entre las medias, las pruebas de rango post hoc o de clasificación de medias permiten determinar como medias difieren. La prueba de rango post hoc identifica subconjuntos homogéneos de medias que no se diferencian entre sí. Las pruebas más comunes en investigación agropecuaria son:

  • Diferencia mínima significativa (Fisher)
  • Prueba de Fisher ajustada (Bonferroni)
  • Diferencia honesta significativa (Tukey)
  • Prueba de rangos múltiples (Duncan)

Para hacer la comparación de medias volvemos a utilizar el set de datos forraje que se utilizó para el ANOVA.

  • Europe1<-read.csv("C:\\Users\\Rafael\\Documents\\Visiting Scientists\\forraje.csv")
  • fix (Europe1)   #Cambiamos de nombre al sitio por bloques
  • ANOVA3 = aov(Rendimiento ~ Bloque+Tratamiento, data=Europe1) #Efecto del tratamiento y los bloques
  • summary (ANOVA3)
  • library (agricolae)
  • #Comparación de medias de Fisher
  • Fisher <- LSD.test (ANOVA3, "Tratamiento") 
  • Fisher
  • #Comparación de medias de Fisher con el ajuste de Bonferroni
  • Bonferroni <- LSD.test (ANOVA3, "Tratamiento", p.adj = "bonferroni")  
  • Bonferroni
  • #Comparación de medias de Tukey
  • Tukey <- HSD.test(ANOVA3,"Tratamiento")  
  • Tukey
  • #Comparación de medias de Duncan
  • Duncan<- duncan.test (ANOVA3, "Tratamiento")   
  • Duncan

Ahora hacemos los gráficos

  • par(mfrow=c(2,2))
  • bar.group (Fisher$groups, density=8, col="blue", border="red" , main="Fisher", ylim=c(0,4000))
  • bar.group (Bonferroni$groups,density=8,col="blue",border="red",main="Bonferroni", ylim=c(0,4000))
  • bar.group(Tukey$groups,density=8,col="blue",border="red",main="Tukey", ylim=c(0,4000))
  • bar.group(Duncan$groups,density=8,col="blue",border="red",main="Duncan", ylim=c(0,4000))