febrero 09, 2013

Análisis de varianza para diseño cuadrado latino en R

Datos tomados de Tabla 2.7: Rendimiento en grano de híbridos de maíz prometedores (A, B, y D) y una variedad control (C). Resultados de experimento con diseño Cuadrado Latino. Procedimientos Estadísticos para Investigaciones Agrícolas. 1984. Gómez, K.A.; Gómez, A.A. Segunda Edición; IRRI-John Wiley & Sons.  pp. 33–39.

# Análisis de varianza para diseño cuadrado latino
# Cargando datos
maiz = read.csv("C:/Users/Administrator/Desktop/tabla2.8.csv", header = T)

# Mostrando contenido de archivo "maiz"
maiz
   hilera columna hibrido yield
1       1       1       B 1.640
2       2       1       C 1.475
3       3       1       A 1.670
4       4       1       D 1.565
5       1       2       D 1.210
6       2       2       A 1.185
7       3       2       C 0.710
8       4       2       B 1.290
9       1       3       C 1.425
10      2       3       D 1.400
11      3       3       B 1.665
12      4       3       A 1.655
13      1       4       A 1.345
14      2       4       B 1.290
15      3       4       D 1.180
16      4       4       C 0.660
str(maiz)
'data.frame':   16 obs. of  4 variables:
 $ hilera : int  1 2 3 4 1 2 3 4 1 2 ...
 $ columna: int  1 1 1 1 2 2 2 2 3 3 ...
 $ hibrido: Factor w/ 4 levels "A","B","C","D": 2 3 1 4 4 1 3 2 3 4 ...
 $ yield  : num  1.64 1.48 1.67 1.56 1.21 ...

# Convirtiendo "hileras" y "columnas" en variables categóricas
maiz$hilera = factor(maiz$hilera)
maiz$columna = factor(maiz$columna)

# Visualizando tratamientos
windows (height=4, width=4)
with(maiz, stripchart(yield ~ columna, xlab = "columna", vert = T))
plot of chunk unnamed-chunk-1
with(maiz, stripchart(yield ~ hilera, xlab = "hilera", vert = T))
plot of chunk unnamed-chunk-1
with(maiz, stripchart(yield ~ hilera, xlab = "hilera", vert = T))
with(maiz, stripchart(yield ~ hibrido, xlab = "hibrido", vert = T))
plot of chunk unnamed-chunk-1


# Análisis de variance
attach(maiz)
maiz.lm = lm(yield ~ hilera + columna + hibrido)
anova(maiz.lm)
Analysis of Variance Table

Response: yield
          Df Sum Sq Mean Sq F value Pr(>F)   
hilera     3  0.030  0.0101    0.47 0.7170   
columna    3  0.827  0.2758   12.77 0.0051 **
hibrido    3  0.427  0.1423    6.59 0.0251 * 
Residuals  6  0.130  0.0216                  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1




No hay comentarios.: