Cap Viii
Cap Viii
Cap Viii
Captulo VIII
Experimentos Factoriales
Las caractersticas del diseo siempre triunfan sobre las caractersticas del anlisis.
G. E. Dallal
1. Introduccin
2. Ventajas y Desventajas
Ventajas:
- Permite obtener ms informacin que en un experimento de un solo factor, ya que se
estudian los efectos principales, los efectos simples, los efectos cruzados y de interaccin
entre los factores.
- Todas las unidades intervienen en la estimacin de los efectos principales y de
interaccin, por lo que el nmero de repeticiones es elevado para estos casos.
Desventajas:
- Se requiere un mayor nmero de unidades experimentales que en los experimentos con un
solo factor.
Ing. Ral Eyzaguirre Prez 112
reyzaguirre@lamolina.edu.pe
Departamento de Estadstica e Informtica
Curso: Mtodos Estadsticos para la Investigacin I Experimentos Factoriales
- Dado que todos los niveles de un factor se combinan con todos los niveles de los otros,
por requerimientos del anlisis estadstico, se tendr que algunas combinaciones que no
son de inters para el investigador, sern tambin incluidas en el experimento.
- El anlisis estadstico y la interpretacin de los resultados son ms complicados que en los
experimentos con un solo factor, y la dificultad aumenta considerablemente conforme ms
factores son incluidos.
3. Notacin y Definiciones
3.1. Factor
Los factores son designados por letras maysculas. Por ejemplo en un experimento en el que
se evalan 3 cantidades de semilla con 4 dosis de nitrgeno por parcela y 2 variedades de
maz destinado a chala, el factor cantidad de semilla se puede denotar por A, el factor dosis de
nitrgeno por B y el factor variedad de maz por C.
Los niveles de un factor son denotados por letras minsculas con subndices. Por ejemplo, las
3 cantidades de semilla podran ser denotadas por a1, a2 y a3, las 4 dosis de nitrgeno por b1,
b2, b3 y b4, y las 2 variedades de maz por c1 y c2.
Una combinacin de letras minsculas con sus respectivos subndices es utilizada para
denotar una combinacin de los niveles de los factores. Por ejemplo la combinacin a2b2c1
denotar el tratamiento conformado por la aplicacin de la cantidad a2 de semilla con la dosis
b2 de nitrgeno y la variedad c1 de maz.
Dependiendo de la naturaleza de los niveles de los factores, estos pueden ser cualitativos o
cuantitativos. En el ejemplo, los factores A y B son cuantitativos y el factor C cualitativo. En
el caso de factores cuantitativos estos pueden ser igualmente espaciados o no. As por
ejemplo, para el factor B, niveles de 0, 10, 20 y 30 kg/parcela y de 10, 20, 40 y 80 kg/parcela
constituiran niveles igualmente espaciados y no igualmente espaciados respectivamente.
Adicionalmente, los factores pueden ser fijos o al azar, dependiendo de la forma en que son
seleccionados sus niveles (ver el captulo II). Un experimento factorial con todos sus factores
fijos corresponder a un modelo I o de efectos fijos, un experimento factorial con todos sus
factores aleatorios corresponder a un modelo II o de efectos aleatorios y un experimento
factorial con algunos factores fijos y otros aleatorios corresponder a un modelo III o de
efectos mixtos. En el desarrollo de este captulo se considerar que todos los factores son
fijos.
Un experimento factorial queda definido por el nmero de factores y niveles de cada factor.
Un experimento factorial puede ser denotado utilizando las letras correspondientes a los
factores antecedidas por el nmero de niveles correspondiente a cada uno. Por ejemplo, el
experimento con 3 niveles del factor A, 4 del factor B y 2 del factor C puede ser denotado por
3A4B2C o simplemente 3x4x2.
Efecto Principal: Es el efecto de un factor en promedio sobre los niveles de los otros
factores.
Efecto Simple: Es el efecto de un factor, en un nivel de los dems factores.
Efecto de Interaccin: Est dado por la variacin que tiene un efecto simple de un factor al
pasar de un nivel a otro de otro factor.
Efecto Cruzado: Est dado por las combinaciones cruzadas de dos factores.
Efectos Simples:
- de A en b1: ES(A(b1)) = a1b1 a2b1 = 54 45 = 9
- de A en b2: ES(A(b2)) = a1b2 a2b2 = 38 56 = -18
- de B en a1: ES(B(a1)) = a1b1 a1b2 = 54 38 = 16
- de B en a2: ES(B(a2)) = a2b1 a2b2 = 45 56 = -11
Efectos Principales:
1 9 18
- de A: EP(A) = [ ES( A(b1 ) + ES( A(b2 )] = = 4.5
2 2
1 16 11
- de B: EP(B) = [ ES( B(a1 ) + ES( B(a2 )] = = 2.5
2 2
Efecto de interaccin:
1 9 + 18
- de AB: EI(AB) = [ ES( A(b1 ) ES( A(b2 )] = = 6.75
2 2 4
1 16 + 11
EI(AB) = [ ES( B(a1 ) ES( B(a2 )] = = 6.75
2 2 4
Efectos cruzados:
- entre a1b1 y a2b2: EC(a1b1 - a2b2) = a1b1 a2b2 = 54 56 = -2
- entre a1b2 y a2b1: EC(a1b2 - a2b1) = a1b2 a2b1 = 38 45 = -7
Cada lnea corresponde a un efecto simple, y la interaccin puede notarse cuando las lneas
tienen pendientes diferentes, esto es, cuando el efecto simple de un factor no es el mismo en
todos los niveles del otro. En el primer grfico se presentan los efectos simples de A en b1
(lnea punteada) y de A en b2 (lnea continua) y es claro que el efecto simple de A depende del
nivel de B (esto es, que existe interaccin). A continuacin se presenta otro ejemplo de
interaccin:
Note que no es necesario que las lneas se crucen para evidenciar una interaccin entre ambos
factores. Un caso de factores sin interaccin aparente sera el siguiente:
Ahora, dado que los tratamientos son generados por las combinaciones entre los niveles de
dos factores, el efecto de los tratamientos se descompone en el efecto del factor A, el efecto
del factor B y el efecto de la interaccin entre los dos factores. As, el modelo aditivo lineal
para un factorial pxq en DCA ser:
Yijk = + i + j + ( )ij + ijk i = 1,..., p j = 1,..., q k = 1,..., rij
donde:
Yijk es el valor o rendimiento observado con el i-simo nivel del factor A, j-simo nivel
del factor B, k-sima repeticin.
es el efecto de la media general.
i es el efecto del i-simo nivel del factor A.
j es el efecto del j-simo nivel del factor B.
()ij es el efecto de la interaccin en el i-simo nivel del factor A, j-simo nivel del
factor B.
ijk es el efecto del error experimental en el i-simo nivel del factor A, j-simo nivel del
factor B, k-sima repeticin.
p es el nmero de niveles del factor A.
q es el nmero de niveles del factor B.
rij es el nmero de repeticiones en el i-simo nivel del factor A, j-simo nivel del factor B.
donde:
k es el efecto del k-simo bloque.
b es el nmero de bloques.
Los supuestos del modelo sern los mismos que para el DCA o DBCA con un solo factor
vistos en captulos anteriores. Los clculos y procedimientos presentados de aqu en adelante
correspondern al caso del experimento factorial pxq en DBCA. El caso del experimento
factorial en DCA es similar al del DBCA y ser ilustrado mediante un ejemplo.
a1 a2
Campos b1 b2 b3 b1 b2 b3
1 1.9 1.8 2.7 1.8 2.9 3.0
2 2.3 2.1 2.4 2.2 2.7 3.2
3 2.0 2.4 2.9 2.0 3.2 2.9
4 2.1 2.9 2.8 2.4 3.5 3.4
Total 8.3 9.2 10.8 8.4 12.3 12.5
donde:
Yijk es el rendimiento de maz en Tm/Ha obtenido con el i-simo tipo de abono, j-sima
dosis, k-simo campo de cultivo.
es el efecto de la media general.
i es el efecto del i-simo tipo de abono.
j es el efecto de la j-sima dosis de abono.
()ij es el efecto de la interaccin en el i-simo tipo de abono, j-sima dosis.
k es el efecto del k-simo campo de cultivo.
ijk es el efecto del error experimental en el i-simo tipo de abono, j-sima dosis, k-
simo campo de cultivo.
p = 2 es el nmero de niveles del factor A.
q = 3 es el nmero de niveles del factor B.
b = 4 es el nmero de bloques.
Los efectos del modelo, , i, j, ()ij y k son estimados de modo que se minimice la
siguiente expresin (Mtodo de Mnimos Cuadrados):
p q b p q b
Q = = (Yijk i j ( )ij k )
2 2
ijk
i =1 j = i k =1 i =1 j = i k =1
i = 0
i =1
j = 0
j =1
( )ij = 0
i =1
( )ij = 0
j =1
k =1
k =0
= Y i = Yi Y j = Y j Y
( )ij = Yij Yi Y j + Y k = Yk Y ijk = Yijk Yij Yk + Y
Ejemplo 2 (Cont.): Con los datos del ejemplo anterior, la media estimada es:
= Y = 2.5625
El efecto estimado de la interaccin entre el nivel 1 del factor A y el nivel 2 del factor B:
( )12 = Y12 Y1 Y2 + Y = 2.3 2.3583 2.6875 + 2.5625 = 0.1833
Las sumas de cuadrados para los factores A y B, para la interaccin, bloques y error se
calculan de la siguiente manera:
p
Yi2
SC(A) = TC
i =1 qb
q
Y2j
SC(B) = TC
j =1 pb
Cuadro ANVA
Hiptesis:
Para el Modelo I (Efectos fijos) las hiptesis son, en trminos de los efectos de los niveles de
los factores las siguientes:
Para el efecto principal de A: H0: i = 0 i
H1: i 0 para al menos algn i
Para el efecto principal de B: H0: j = 0 j
H1: j 0 para al menos algn j
Para el efecto de la interaccin AB: H0: ()ij = 0 i, j
H1: ()ij 0 para al menos algn i, j
Para el Modelo II (Efectos aleatorios) las hiptesis sern planteadas en trminos de la
variancia de los niveles de los factores:
Estadstico de Prueba:
CM(A)
Para el efecto principal de A: Fc = ~ F( gl( A),gl(Error) )
CM(Error)
CM(B)
Para el efecto principal de B: Fc = ~ F( gl( B ),gl(Error) )
CM(Error)
CM(AB)
Para el efecto de la interaccin AB: Fc = ~ F( gl( AB ),gl(Error) )
CM(Error)
Regla de Decisin:
Las hiptesis nulas se rechazan con un nivel de significacin si el Fc resulta mayor que el
valor de tabla F(1 ) con los grados de libertad correspondientes a cada caso.
q
Y2j
SC( B) = TC
j =1 pb
16.7 2 21.52 23.32 61.52
= + + = 2.91
(2)(4) (2)(4) (2)(4) (2)(3)(4)
Cuadro ANVA
Fuentes de Variacin gl SC CM Fc
Bloques 3 0.8046 0.2682
A 1 1.0004 1.0004 18.81
B 2 2.9100 1.4550 27.35
AB 2 0.5633 0.2817 5.30
Error Experimental 15 0.7979 0.0532
Total 23 6.0763
CME 0.0532
cv = = = 9%
Y 2.5625
Este anlisis debe ser efectuado en el caso que la interaccin resulte significativa y consiste en
evaluar a cada factor en cada uno de los niveles del otro. Las hiptesis a contrastar en este
caso, asumiendo un Modelo I (Efectos fijos) son las siguientes:
1. Para el efecto simple de A en el j-simo nivel de B:
H0: 1 j = 2 j = " pj
H1: Al menos un ij es diferente.
Los grados de libertad para cada efecto simple sern iguales a los grados de libertad del
correspondiente efecto principal y las sumas de cuadrados son calculadas de acuerdo con las
siguientes frmulas:
1. Para el factor A en el j-simo nivel del factor B:
p
Yij2 Y2j
SC(Abj) =
i =1 b pb
Para cada efecto simple el estadstico de prueba Fc se calcula dividiendo el cuadrado medio
del efecto simple entre el cuadrado medio del error. El efecto ser significativo con un nivel
de significacin si es que el Fc es mayor que el valor F de tabla con los grados de libertad
del efecto y del error.
p
Yi 32 Y23 10.82 12.52 23.32
SC(Ab3) = = + = 0.36125
i =1 b pb 4 4 (2)(4)
q
Y12j 2
Y1 8.32 9.22 10.82 28.32
SC(Ba1) = = + + = 0.80167
j =1 b qb 4 4 4 (3)(4)
q
Y22j Y22 8.42 12.32 12.52 33.22
SC(Ba2) = = + + = 2.67167
j =1 b qb 4 4 4 (3)(4)
Fuentes de Variacin gl SC CM Fc
A b1 1 0.00125 0.00125 0.02 n.s.
A b2 1 1.20125 1.20125 22.58 *
A b3 1 0.36125 0.36125 6.79 *
B a1 2 0.80167 0.40083 7.54 *
B a2 2 2.67167 1.33583 25.11 *
Error Experimental 15 0.79792 0.05319
Total 23 6.07625
Los efectos simples del factor A son comparados con el valor de tabla F( 0.95, 1,15) = 4.54 y
los efectos simples del factor B con F( 0.95, 2,15) = 3.68. Note que solo el efecto simple de A en
b1 resulta no significativo. Las conclusiones en este experimento seran las siguientes:
- No existe suficiente evidencia estadstica para aceptar que con los dos tipos de abono
se obtengan resultados diferentes en el rendimiento de maz cuando se aplican en la
dosis b1 (20 kg/ha).
Ing. Ral Eyzaguirre Prez 124
reyzaguirre@lamolina.edu.pe
Departamento de Estadstica e Informtica
Curso: Mtodos Estadsticos para la Investigacin I Experimentos Factoriales
- Existe suficiente evidencia estadstica para aceptar que con los dos tipos de abono se
obtienen resultados diferentes en el rendimiento de maz cuando se aplican en las dosis
b2 (30 kg/ha) y b3 (40 kg/ha).
- Existe suficiente evidencia estadstica para aceptar que con al menos una de las dosis
se obtienen resultados diferentes en el rendimiento de maz tanto con el abono a1 como
con el a2.
Estos resultados pueden apreciarse en el siguiente grfico:
Para comparar las medias de los niveles i y j de un factor sobre todos los niveles del otro
utilice las siguientes frmulas para las desviaciones estndar:
Para comparar las medias de los niveles k y l de un factor en un nivel del otro utilice las
siguientes frmulas para las desviaciones estndar:
CME 0.05319
ALS(T) = AES(T) = 3.67 = 0.4232
b 4
A continuacin se presentan los resultados para las 3 comparaciones:
El valor de tabla con = 5% y 15 grados de libertad para el error es t(0.05, 15) = -1.753.
Como el estadstico de prueba resulta menor que el valor de tabla, se rechaza H0 y se
concluye que existe suficiente evidencia estadstica para aceptar que con la dosis 3 se
obtienen mejores resultados que con la dosis 1 cuando se aplica el abono tipo 2.
Ejemplo 3 (Experimento factorial en DCA): Los siguientes datos muestran los aumentos
de peso en gramos de ratas machos bajo 6 tratamientos de alimentacin en un experimento
conducido en un DCA. Los factores en estudio y los resultados del experimento son:
A: Fuente de protena (a1 = Res, a2 = cereal, a3 = cerdo)
B: Dosis de protena (b1 = Baja, b2 = alta).
a1 a2 a3
b1 b2 b1 b2 b1 b2
88 83 91 107 64 93
74 102 81 95 83 78
90 108 75 97 74 95
83 104 101 80 73 97
335 397 348 379 294 363
Modelo Aditivo Lineal: En este caso el diseo utilizado es un DCA por lo que el modelo
aditivo lineal ser el siguiente:
Yijk = + i + j + ( )ij + ijk i = 1,..., p j = 1,..., q k = 1,..., r
donde:
Yijk es la ganancia de peso obtenida con la i-sima fuente de protena en la j-sima dosis,
k-sima repeticin.
es el efecto de la media general.
i es el efecto de la i-sima fuente de protena.
j es el efecto de la j-sima dosis de protena.
()ij es el efecto de la interaccin en la i-sima fuente de protena, j-sima dosis.
ijk es el efecto del error experimental en la i-sima fuente de protena, j-sima dosis, k-
sima repeticin.
p = 3 es el nmero de niveles del factor A.
q = 2 es el nmero de niveles del factor B.
r = 4 es el nmero de repeticiones.
= Y = 88.17
Anlisis de Variancia:
p 2
q r
Y
SC(Total) = Y 2
ijk
i =1 j =1 k =1 pqr
21162
= (882 + 742 + ... + 97 2 ) = 3329.33
(3)(2)(4)
p q
Yij2
SC(Comb. AB) = TC
i =1 j =1 r
3352 397 2 3632 21162
= + +"+ = 1635.33
4 4 4 (3)(2)(4)
p
Yi2
SC( A) = TC
i =1 qr
Cuadro ANVA
Fuentes de Variacin gl SC CM Fc
A 2 439.58 219.79 2.34
B 1 1093.50 1093.50 11.62
AB 2 102.25 51.13 0.54
Error Experimental 18 1694.00 94.11
Total 23 3329.33
CME 94.11
cv = = = 11%
Y 88.17
Ejemplo 3
General Linear Model: Y versus A, B
Factor Type Levels Values
A fixed 3 1 2 3
B fixed 2 1 2
Ejercicios
v1 v2 v3
Rep. a1 a2 a3 a1 a2 a3 a1 a2 a3
1 4.5 6.3 3.5 4.1 4.0 3.6 4.2 4.8 3.9
2 4.4 6.9 3.5 3.9 3.5 3.1 4.3 4.3 3.6
3 4.2 6.4 4.0 3.6 4.0 3.5 3.8 3.9 4.0
4 3.9 6.5 3.2 4.1 4.1 3.9 4.7 4.2 4.1
a) Presente el modelo aditivo lineal e interprete cada uno de sus componentes en trminos
del enunciado.
b) Efecte el anlisis grfico de la interaccin.
c) Efecte el anlisis de variancia. Analice los efectos principales o simples segn
corresponda.
d) Efecte la prueba de Tukey para evaluar si existen diferencias entre los venenos cuando se
aplica el antdoto a2.
e) Se cree que el antdoto a2 es ms efectivo que el a1 para contrarrestar al veneno v1.
Efecte la prueba correspondiente.
s1 (Machos) s2 (Hembras)
Bloque d1 d2 d3 d1 d2 d3
1 45 57 57 40 49 52
2 50 53 55 45 52 56
3 42 63 65 43 53 55
4 46 60 61 39 50 55
a) Presente el modelo aditivo lineal e interprete cada uno de sus componentes en trminos
del enunciado
b) Efecte el anlisis grfico de la interaccin.
c) Efecte el anlisis de variancia. Analice los efectos principales o simples segn
corresponda.
d) Efecte la prueba de Tukey para evaluar si existen diferencias entre las dosis.
envoltura, sobre las ventas de uno de sus productos. Se utilizaron dos tipos de material
(a1=Papel encerado y a2=Plstico), en tres colores (b1=Amarillo, b2=Rojo y b3=Verde). Se
seleccionaron 4 supermercados para el experimento y despus de estar el producto en el
mercado por una semana se registr la venta total (en miles de soles) para cada una de las
6 combinaciones. Los resultados se muestran a continuacin:
a) Presente el modelo aditivo lineal e interprete cada uno de sus componentes en trminos
del enunciado
b) Efecte el anlisis grfico de la interaccin.
c) Efecte el anlisis de variancia. Analice los efectos principales o simples segn
corresponda.
d) Efecte la prueba de Tukey donde sea necesario.
4. Se realiz un experimento para evaluar el efecto de tres densidades de siembra (d1, d2, d3)
con tres variedades de frijol (v1, v2 y v3), en el rendimiento de frijol en kg/parcela. El
diseo utilizado para el experimento fue un DBCA.
v1 v2 v3
Bloque d1 d2 d3 d1 d2 d3 d1 d2 d3
I 10.05 9.66 9.14 10.71 10.35 11.42 9.03 10.46 13
II 8.71 8.45 9.02 9.45 10.24 12.91 8.54 10.5 10.1
III 9.9 8.05 8.01 9.25 11.1 11.5 7.24 8.85 11.57
a) Presente el modelo aditivo lineal e interprete cada uno de sus componentes en trminos
del enunciado
b) Efecte el anlisis grfico de la interaccin.
c) Efecte el anlisis de variancia. Analice los efectos principales o simples segn
corresponda.
d) Efecte la prueba de Tukey en donde sea necesario.
e) A partir de los resultados obtenidos en las preguntas anteriores, presente sus
recomendaciones.
a1 a2 a3
Campo b1 b2 b1 b2 b1 b2
1 9.56 8.26 9.18 8.90 8.26 9.82
2 9.32 8.50 8.86 8.50 8.64 9.84
3 8.96 8.42 8.22 9.82 8.10 9.7
4 8.78 8.26 8.70 9.78 8.72 10.04
a) Presente el modelo aditivo lineal e interprete cada uno de sus componentes en trminos
del enunciado.
b) Efecte el anlisis grfico de la interaccin.
c) Efecte el anlisis de variancia. Analice los efectos principales o simples segn
corresponda.
d) Mediante la prueba DLS compare a1b2 con a3b2.
6. Con la finalidad de estudiar el efecto de tres niveles de Nitrgeno (a1, a2, a3) y dos niveles
de fsforo (b1, b2), en el cultivo de una variedad de papa se realiz un experimento con un
arreglo factorial conducido en el DCA con 4 repeticiones. Los resultados obtenidos en
kg/parcela son los siguientes:
a1 a2 a3
Repeticin b1 b2 b1 b2 b1 B2
1 31 43 42 45 48 51
2 32 41 38 46 50 47
3 34 43 36 44 48 50
4 35 39 41 43 51 52
a) Presente el modelo aditivo lineal e interprete cada uno de sus componentes en trminos
del enunciado.
b) Efecte el anlisis grfico de la interaccin.
c) Efecte el anlisis de variancia. Analice los efectos principales o simples segn
corresponda.
d) Efecte la prueba de Tukey en donde sea necesario.