Se llevó a cabo un experimento bajo un Diseño Completamente al Azar (DCA) con número desigual de repeticiones, con el objetivo de comparar los rendimientos (tn/ha) de cuatro variedades de aceituna de mesa: Sevillana, Ascolana, Empeltre e Ilo.
El presente documento presenta el análisis estadístico completo, ordenado de la siguiente manera:
Antes de interpretar los resultados del ANOVA es fundamental verificar que se cumplan los supuestos de:
##
## Shapiro-Wilk normality test
##
## data: residuos
## W = 0.94672, p-value = 0.2299
Interpretación:
El valor-p de la prueba de Shapiro-Wilk es 0.2299.
Como p > 0.05, no se rechaza la hipótesis nula de normalidad. Los
residuales se distribuyen de forma aproximadamente normal.
Diagnóstico de normalidad de residuales
Interpretación:
El valor-p de la prueba de Levene es 0.1623.
Como p > 0.05, no se rechaza la homogeneidad de varianzas. Se cumple
el supuesto de homocedasticidad.
##
## Bartlett test of homogeneity of variances
##
## data: rendimiento by variedad
## Bartlett's K-squared = 8.4254, df = 3, p-value = 0.03799
Residuales vs valores ajustados
Conclusión de supuestos:
Los residuales cumplen razonablemente los supuestos de normalidad y
homogeneidad de varianzas, por lo que es válido proceder con el ANOVA
paramétrico.
| Df | Sum Sq | Mean Sq | F value | Pr(>F) | |
|---|---|---|---|---|---|
| variedad | 3 | 82.8898 | 27.6299 | 16.4218 | 0 |
| Residuals | 20 | 33.6502 | 1.6825 | NA | NA |
Interpretación del ANOVA:
Distribución del rendimiento por variedad de aceituna
Medias de rendimiento con barras de error estándar
| Variedad | n | Media | DE | EE | Mín | Máx | CV (%) |
|---|---|---|---|---|---|---|---|
| Sevillana | 6 | 8.417 | 2.031 | 0.829 | 4.9 | 10.5 | 24.133 |
| Ascolana | 7 | 5.886 | 1.220 | 0.461 | 4.7 | 8.2 | 20.726 |
| Empeltre | 5 | 4.900 | 0.474 | 0.212 | 4.3 | 5.6 | 9.680 |
| Ilo | 6 | 9.633 | 0.799 | 0.326 | 8.7 | 10.6 | 8.296 |
Observaciones descriptivas:
Dado que el ANOVA resultó significativo, se realizan pruebas de comparación de medias.
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = rendimiento ~ variedad, data = datos)
##
## $variedad
## diff lwr upr p adj
## Ascolana-Sevillana -2.5309524 -4.5508023 -0.5111025 0.0109566
## Empeltre-Sevillana -3.5166667 -5.7150721 -1.3182613 0.0012178
## Ilo-Sevillana 1.2166667 -0.8794306 3.3127640 0.3880384
## Empeltre-Ascolana -0.9857143 -3.1115449 1.1401164 0.5748123
## Ilo-Ascolana 3.7476190 1.7277692 5.7674689 0.0002396
## Ilo-Empeltre 4.7333333 2.5349279 6.9317387 0.0000378
Intervalos de confianza de las diferencias (Tukey HSD)
| Media | Grupo | |
|---|---|---|
| Ilo | 9.633333 | a |
| Sevillana | 8.416667 | a |
| Ascolana | 5.885714 | b |
| Empeltre | 4.900000 | b |
## === Prueba LSD (Fisher) ===
##
## Study: modelo ~ "variedad"
##
## LSD t Test for rendimiento
##
## Mean Square Error: 1.682512
##
## variedad, means and individual ( 95 %) CI
##
## rendimiento std r se LCL UCL Min Max Q25
## Ascolana 5.885714 1.2198751 7 0.4902641 4.863041 6.908387 4.7 8.2 5.000
## Empeltre 4.900000 0.4743416 5 0.5800883 3.689957 6.110043 4.3 5.6 4.700
## Ilo 9.633333 0.7991662 6 0.5295457 8.528720 10.737946 8.7 10.6 8.925
## Sevillana 8.416667 2.0311737 6 0.5295457 7.312054 9.521280 4.9 10.5 7.675
## Q50 Q75
## Ascolana 5.60 6.35
## Empeltre 4.90 5.00
## Ilo 9.75 10.20
## Sevillana 8.95 9.70
##
## Alpha: 0.05 ; DF Error: 20
## Critical Value of t: 2.085963
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Treatments with the same letter are not significantly different.
##
## rendimiento groups
## Ilo 9.633333 a
## Sevillana 8.416667 a
## Ascolana 5.885714 b
## Empeltre 4.900000 b
##
## === Prueba de Duncan ===
##
## Study: modelo ~ "variedad"
##
## Duncan's new multiple range test
## for rendimiento
##
## Mean Square Error: 1.682512
##
## variedad, means
##
## rendimiento std r se Min Max Q25 Q50 Q75
## Ascolana 5.885714 1.2198751 7 0.4902641 4.7 8.2 5.000 5.60 6.35
## Empeltre 4.900000 0.4743416 5 0.5800883 4.3 5.6 4.700 4.90 5.00
## Ilo 9.633333 0.7991662 6 0.5295457 8.7 10.6 8.925 9.75 10.20
## Sevillana 8.416667 2.0311737 6 0.5295457 4.9 10.5 7.675 8.95 9.70
##
## Groups according to probability of means differences and alpha level( 0.05 )
##
## Means with the same letter are not significantly different.
##
## rendimiento groups
## Ilo 9.633333 a
## Sevillana 8.416667 a
## Ascolana 5.885714 b
## Empeltre 4.900000 b
Con base en el análisis realizado se concluye lo siguiente:
Los supuestos del ANOVA se cumplen de manera adecuada (normalidad y homogeneidad de varianzas), por lo que los resultados del análisis paramétrico son válidos.
Existen diferencias altamente significativas entre los rendimientos de las cuatro variedades de aceituna de mesa (ANOVA, p < 0.05).
Según la prueba de Tukey HSD, las variedades se agrupan de la siguiente forma (letras distintas indican diferencias significativas a α = 0.05):
## rendimiento groups
## Ilo 9.633333 a
## Sevillana 8.416667 a
## Ascolana 5.885714 b
## Empeltre 4.900000 b
La variedad Ilo muestra el mayor rendimiento promedio, mientras que Empeltre presenta el menor desempeño bajo las condiciones del experimento.
Se recomienda considerar la variedad Ilo cuando el objetivo principal sea maximizar el rendimiento en tn/ha.
Documento generado en R Markdown – Análisis estadístico de DCA con repeticiones desiguales. ```