Análisis de varianza para DCA (usando R)

Análisis de varianza para Diseño Completamente al Azar (DCA). Ejemplo de Análisis de Varianza tomado de Procedimientos Estadísticos para Investigaciones Agrícolas. 1984. Gómez, K.A.; Gómez, A.A. Segunda Edición; IRRI—John Wiley & Sons. pp. 14-17.

Tabla 2.1 Rendimiento de arroz en granos usando diferentes foliares e insecticidas granulados para el control de Saltamonte Marrón, (Nilaparvata lugens [Stal]) y Taladrador del Arroz Amarillo (Scirpophaga incertulas [Walker]); experimento CRD con 4 (r) repeticiones y 7 (t) tratamientos.

Tratamiento	Rendimiento de grano, kg/ha				Total Tratamiento	Media Tratamientos
Foliar A (1 kg)	2,537	2,069	2,104	1,797	8,507	2126.75
Foliar B (2 kg)	3,366	2,591	2,211	2,544	10,712	2678.00
Insecticida A	2,536	2,459	2,827	2,385	10,207	2551.75
Insecticida B	2,387	2,453	1,556	2,116	8,512	2128.00
Insecticida C	1,997	1,679	1,649	1,859	7,184	1796.00
Insecticida D	1,796	1,704	1,904	1,320	6,724	1681.00
Control	1,401	1,516	1,270	1,077	5,264	1316.00
Grand total (G)					57,110
Grand media						2039.64

# Analisis de varianza para Diseno Completamente al Azar (DCA)

# Cargando datos de table2.1.txt
arroz <- span=""> read.table("C:/Users/Administrator/Desktop/table2.1.txt", header = T)



#Medias por tratamiento
aggregate(arroz$RENDIMIENTO, list(TRATAMIENTO = arroz$TRATAMIENTO), mean)

    TRATAMIENTO    x
1       Control 1316
2 Foliar.A(1kg) 2127
3 Foliar.B(2kg) 2678
4 Insecticida.A 2552
5 Insecticida.B 2128
6 Insecticida.C 1796
7 Insecticida.D 1681

#Boxplot 
require(ggplot2)

Loading required package: ggplot2

ggplot(arroz, aes(x = TRATAMIENTO, y = RENDIMIENTO)) +
  geom_boxplot(fill = "grey80", colour = "blue") +
  scale_x_discrete() + xlab("Tratamiento") +
  ylab("Rendimiento, kg/ha")

# Analisis de varianza 
anova.arroz = aov(RENDIMIENTO ~ TRATAMIENTO, data=arroz)
summary(anova.arroz)

            Df  Sum Sq Mean Sq F value  Pr(>F)    
TRATAMIENTO  6 5587175  931196    9.83 3.3e-05 ***
Residuals   21 1990238   94773                    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1