Tema: Ejercicios de diseño factorial A*B
18. A continuación se muestra parte del ANOVA para un diseño factorial 3X5 con dos
réplicas, el factor A con tres niveles y B con cinco.
Factores de Suma de Grados de Cuadrado Fca Ftabl
variación cuadrados libertad medio l a
A 800
B 900
AB 300
Error 400
Total
a=3, b=5, n=2
Grados de libertad para A = a-1
Grados de libertad para B = b-1
Grados de libertad para A*B = (a-1)(b-1)
Grados de libertad para Error = ab(n-1)
Grados de libertad para Total = abn-1
Cuadrado medio A = (suma de cuadrados de A)/(grados de libertad de A)
Cuadrado medio B = (suma de cuadrados de B)/(grados de libertad de B)
Cuadrado medio AB = (suma de cuadrados de AB)/(grados de libertad de AB)
Cuadrado medio Error = (suma de cuadrados de Error)/(grados de libertad de Error)
FcalA = (cuadrado medio de A)/(cuadrado medio del error)
FcalB = (cuadrado medio de B)/(cuadrado medio del error)
FcalAB = (cuadrado medio de AB)/(cuadrado medio del error)
A) Suponga efectos fijos, anote el modelo estadístico apropiado y formule las hipótesis a
probar para este experimento.
Modelo estadístico
❑❑ ❑❑❑❑( )❑❑❑
Si: i=3
j=5
k =
❑❑ ❑❑❑❑( )❑❑❑ 2
Hipótesis H O =no hay efecto factor A=0
❑❑
❑❑
❑❑
❑❑
❑❑
B) Agregue en esta tabla los grados de libertad, el cuadro medio y la razón F para cada una
de las fuentes de variación.
Factores de Suma de Grados de Cuadrado Fca Ftabl
variación cuadrados libertad medio l a
A 800 2 400 15 3,68
B 900 4 225 8,4 3,06
AB 300 8 37,5 1,4 2,64
Error 400 15 26,7
Total 2400 29
C) Con la ayuda de un software calcule el valor-p para los efectos.
D) Fcal > Ftabla (existe efecto del factor)
Fcal < Ftabla (no hay efecto del factor)
FcalA > Ftabla (Rechazo Ho. Existe efecto del factor A)
FcalB > Ftabla (Rechazo Ho. Existe efecto del factor B)
FcalAB < Ftabla (Se cumple Ho. No existe efecto de la interacción de los factores A y B
Como conclusión no existe diferencia significativa de la interacción entre los factores A y
B, ya que el Fcal de la interacción es menor al Ftabla. Mientras que el Fcal del factor A y
del factor B independientemente tienen valores mayores al Ftabla, lo cual nos conlleva a
rechazar la hipótesis nula y posteriormente a comprobar la hipótesis alternativa.
19. Se corre un diseño factorial de 3*2 con 10 réplicas para investigar el
hinchamiento del catalizador después de la extrusión den la fábrica de botellas
de polietileno de alta densidad. El catalizador se utiliza en la obtención de dicho
polietileno. Los factores investigados son: molde (con dos niveles) y B:
catalizador (con tres niveles). Los datos obtenidos se muestran en la siguiente
tabla.
Catalizador
❑❑ ❑❑ ❑❑
93 92 93 92 90 92 95 94 94
❑❑ 92 91 90 94 91 91 94 97 96
90 90 90 92 94 95
91 91 91 92 94 96
Molde 88 88 87 90 89 88 91 97 91
❑❑ 88 87 88 88 90 89 90 89 91
87 87 88 89 92 90
87 87 88 88 90 91
a) Plantee la hipótesis de interés en este problema y el modelo estadístico correspondiente.
b) Construya una tabla de análisis de varianza y determine cuales efectos están activos.
c) Dibuje las tablas de medias para los dos efectos principales con los métodos LSD y de
Tukey. Compare los resultados de ambos métodos.
d) Haga una gráfica de interacción con intervalos de confianza sobrepuestos.
e) Determine cuál es el mejor tratamiento. ¿Cuál es el hinchamiento predicho en el mejor
tratamiento?
f) Verifique los supuestos de normalidad y varianza constante.
g) Utilice la gráfica de residuos contra los factores para detectar los posibles efectos sobre la
dispersión del hinchamiento. ¿En cuál molde parece que es menor la dispersión?
Suma total
Y= 5450 d=2 b=3 n=10
( )❑ ()❑ ( )❑
❑❑
❑ ❑
( 1787 )2+ ( 1802 )2+ ( 1861 )2 (5450 )2
SC B= − =153.0333
2∗10 60
(913)2+(915)2 +(949)2 +(874)2+(887)2 +( 912)2 ( 5450 )2
SC AB= − −180.26−153.0333=3.4343
10 60
()❑
❑❑
❑
SC E=424.333333−180.2666−153.0333−3.4343=87.6
Fv SC GL CM ❑❑ valor – p
Efecto A 180.2666 1 180.267 111.12 0.000
Efecto B 153.0333 2 76.517 47.17 0.000
Efecto AB 3.4343 2 1.717 1.06 0.0354
Error 87.6 54 1.622
Total 424.333333 59
Ejercicio 21.
Se desea investigar de qué manera afecta el tiempo de curado y el tipo del acelerante a la
resistencia de caucho vulcanizado. Se realiza un experimento y se obtienen los siguientes
datos:
a) Señale el nombre del diseño de experimento utilizado y su modelo estadístico.
• Diseño factorial mixto: AxB, 3x3, n=2
• Y=µ+αi+βi+(α*β)ij+Eijk
b) Formule claramente todas las hipótesis que se pueden probar.
Hipótesis
• Tiempo de curado
Ho: El tiempo de curado no interviene en la resistencia del caucho vulcanizado.
H1: El tiempo de curado interviene en la resistencia del caucho vulcanizado.
• Acelerante
Ho: El acelerante no afecta en la resistencia del caucho vulcanizado.
H1: El acelerante afecta en la resistencia del caucho vulcanizado.
• Interacción
Ho: El efecto del tiempo de curado no depende del acelerante.
H1: El efecto del tiempo de curdo depende del acelerante.
c) Realice el análisis estadístico apropiado para probar las hipótesis que formuló.
ACELERANTE
A B C Total Yi.. Medias Yi..
3900 4300 3700
40 7500 8000 7800 23300 3883,33
3600 3700 4100
4100 4200 3900
60 7600 8100 7900 23600 3933,33
3500 3900 4000
4000 4300 3600
80 7800 7900 7400 23100 3850,00
3800 3600 3800
Total Y.j.= 22900 24000 23100 Y..=70000
Media Y.j.= 3816,67 4000 3850
❑ ❑
❑❑ ❑ ❑
❑ ❑
❑ ❑ ❑
❑ ❑ () () ( ) ❑❑
❑ ❑ ❑
❑ ❑ ❑
❑ ❑ () () ( ) ❑❑
❑ ❑ ❑
❑ ❑ ❑ ❑❑
❑ ❑ ❑
d) ¿Hay algún tiempo de cura que es mejor para aumentar la resistencia? Argumente
su respuesta
LSD A =t CM E
. 820000
(α ¿¿2 ,ab (n−1))×
√ nA √
=t (0,025 ;9) ×
6
=2,26 × 123,23=278,76 ¿
y´40=3883,33 ± 278,76 LS 40+¿ 4162,09 LI 40−¿ 3604,57
y´60=3933,33 ± 278,76 LS60 +¿ 4212,09 LI 60−¿ 3654,57
y´80=3850 ± 278,76 LS80 +¿ 4128,76 LI 80−¿ 3571,24
Grfico de A
4400
4162.09 4212.09
4200 4128.76
Resistencia al caucho
4000 3883.33 3933.33
3850
3800 3654.57
3604.57 3571.24
3600
3400
3200
0.5 1 1.5 2 2.5 3 3.5
Niveles A
Como podemos observar en la gráfica los tres tiempos de cura serian iguales ya que ay un excelente
traslape entre los tres tiempos.
e) ¿Algún acelerante es mejor? Explique.
LSD B =t CM E 820000
(α ¿¿2 , ab(n−1))×
√ nB √
=t (0,025 ;9 )×
6
=2,26 ×123,23=278,76 ¿
y´A =3816,67 ±278,76 LS A +¿ 4095,43 LI A−¿=3537,91¿
y´B =4000 ± 278,76 LS B +¿ 4278,76 LI B−¿3721,24
y´C =3850 ± 278,76 LSC +¿ 4128,76 LI C −¿ 3571,24
Grafica de B
4500 4278.76 4128.76
4095.43 4000
3816.67 3721.24 3850
4000 3537.91 3571.24
3500
Resistencia al caucho
3000
2500
2000
1500
1000
500
0
0.5 1 1.5 2 2.5 3 3.5
Acelerantes
Los tres tipos de acelerantés son iguales al existir un excelente traslape entre ellos como se puede
observar en la gráfica.
f) ¿Hay alguna combinación de tiempo y acelerante que sea mejor?
Grafico de AxB
4100 4050
4000
4000 3950 3950
Resistencia al caucho
3900 3900
3900
3800
3800 3750
3700
3700
3600
3500
0.5 1 1.5 2 2.5 3 3.5
Tiempo de curado
La mejor resistencia al caucho seria a un tiempo de 60 min con el acelerarte del tipo B.
22. En una fábrica de aceites vegetales comestibles, la calidad resulta afectada por la
cantidad de impurezas dentro del aceite, ya que estas causan oxidación, y ello repercute
a su vez en las características de sabor y color del producto final. El proceso de blanqueo
es el responsable de eliminar tales impurezas y una forma de medir la eficacia es
midiendo el color del aceite. Para generar una primera aproximación a la solución del
problema se decide estudiar el efecto de la temperatura y el porcentaje de arcilla en el
color del aceite inicialmente a nivel laboratorio. El diseño y los datos de las pruebas
experimentales se muestran a continuación.
Factor A = Temperatura (A1: 90, A2: 100, A3: 110).
Factor B = % de arcilla (B1: 0,8; B2: 0,9; B3: 1,0; B4: 1,1).
Porcentaje de arcilla
Temperatura 0,8 0,9 1,0 1,1
90 5, 5, 5,4 5, 4,9 5, 4, 4,4
8 0 5 1 5
100 5, 4, 4,8 4, 4,6 4, 4, 4,3
0 9 7 4 1
110 4, 4, 4,4 4, 4,1 4, 3, 3,6
7 6 4 0 7
A) Construya el modelo estadístico y formule las hipótesis pertinentes.
Modelo estadístico
Y ijk =μ+α i+ β j+ ( αβ )ij + ϵ ijk
i=3
Si:
j=4
k =
2
Y 3 4 2=μ+ α 3 + β 4 + ( αβ )3 4 +ϵ 3 4 2
Hipótesis
H O =no hay efecto temperatura=0
H A =efecto temperatura ≠ 0
H O =no hay efecto % de arcilla=0
H A =efecto % de arcilla ≠ 0
H O =no hay efecto interacción=0
H A =efecto interacción ≠0
B) ¿Cuál es el nombre del diseño utilizado?
Diseño factorial A*B
C) Por lo general, en condiciones reales se utiliza 1.1% de arcilla y 100 grados de
temperatura. ¿Por qué cree que se eligieron precisamente esos niveles de prueba para el
experimento?
Porque este tipo de tratamiento ha funcionado hasta el momento.
D) Realice un análisis de varianza para probar las hipótesis y obtenga conclusiones.
E)
Y= 111.8 a=3 b=4 n=2
❑ ❑
❑ ❑ ❑
❑
❑ ❑
❑ ❑ ❑❑ ❑❑❑❑ ❑❑
❑ ❑ ❑
❑ ❑ ❑ ❑
❑ ❑ () () ( ) ( ) ❑❑
❑ ❑ ❑
❑ ❑ ❑ ❑❑
❑ ❑ ❑
Pruebas de efectos inter-sujetos
Variable dependiente: Calidad de aceites vegetales
Tipo III de suma Media
Origen de cuadrados gl cuadrática F Sig.
Modelo 528,780a 12 44,065 5287,800 ,000
Temperatura 4,041 2 2,020 242,450 ,000
PorArci 3,702 3 1,234 148,067 ,000
Temperatura * PorArci ,236 6 ,039 4,717 ,011
Error ,100 12 ,008
Total 528,880 24
a. R al cuadrado = 1.000 (R al cuadrado ajustada = 1.000)
Se puede concluir que en cuanto a la interacción A*B no hay efecto debido a que
el valor de p es mayor a 0,05 por tanto se acepta la hipótesis nula, mientras que para cada factor
es decir A (temperatura) y B (% de arcilla), existe efecto de los factores, debido a que el
valor de p es menor a 0,05 de la misma forma se rechaza la hipótesis nula y compruebo la
hipótesis alternativa.
F) Apoyándose en las gráficas de efectos. ¿Cuál es la relación general entre el color y los
factores controlados en su rango de experimentación?
La nitidez del color es menor cuando los niveles en ambos factores son altos y es menor en
su nivel más bajo.
G) A partir de la gráfica de interacciones, ¿cree que haya un efecto no lineal?
No, los efectos se consideran lineales.
H) Considerando que el nivel mínimo aceptable de blancura es de 4,8. ¿Qué tratamiento
utilizaría?
Se utilizaría un tratamiento de 0,8% de arcilla a una temperatura de 90, 100 y 110 ºC, para
obtener un valor de blancura mínima de 4,8; de la misma forma el tratamiento con 0,9%
de arcilla a una temperatura de 90 y 100ºC proporciona un nivel de blancura mínimo de
4,8.
I) ¿Valdría la pena planear el estudio a nivel planta?
Sí, para observar resultados más notorios.
J) ¿Qué cambio le haría al experimento si lo corriera a nivel planta?
Se añadiría más réplicas.
23. En un laboratorio de microbiología se realiza un experimento para investigar si influye el
tipo de verdura (lechuga-L, cilantro-C, zanahoria-Z) y la temperatura (8 y 20ºC) de
almacenamiento en la sobrevivencia del Vibrio cholerae. Se hicieron varias réplicas. A
continuación se muestra el porcentaje de sobrevivencia obtenido después de 24 horas
de inoculado el alimento:
Aliment Temperaturas
o 8ºC 20ºC
15, 33, 42, 12, 41, 35, 23,
L 13,1 35,5 11,1 6,2 28,5 25,0 79,0 41,6
0 6 0 8 0 9 8
19, 66, 11, 49, 84, 68, 30, 11,
C 19,0 66,6 11,0 49,0 68,7 30,5 11,0 20,0
0 6 0 0 3 7 5 0
0, 0, 25, 16, 16, 15,
Z 1,2 1,2 2,0 0,1 0,3 0,2 1,0 0,4 21,8 20,1 13,3 25,2
2 3 8 0 0 4
A) Señale el nombre del diseño empleado y formule las hipótesis que pueden ser
probadas.
El nombre es Diseño factorial A*B
H O =no hay efecto alimento=0
H A =efecto alimento ≠ 0
H O =no hay efecto temperatura=0
H A =efecto temperatura ≠ 0
H O =no hay efecto interacción=0
H A =efecto interacción ≠0
Factor A = tipo de verdura (A1: lechuga, A2: cilantro, A3: zanahoria)
Factor B = temperatura (B1: 8ºC, B2: 20ºC)
B) Haga un análisis de varianza e interprételo con detalle.
Pruebas de efectos inter-sujetos
Variable dependiente: Datos
Tipo III de suma Media
Origen de cuadrados gl cuadrática F Sig.
Temperatura 2001,415 1 2001,415 6,361 ,016
Alimento 7990,363 2 3995,182 12,697 ,000
Temperatura * Alimento 230,675 2 115,338 ,367 ,695
Error 13215,040 42 314,644
Total 54494,200 48
a. R al cuadrado = .757 (R al cuadrado ajustada = .723)
Se puede determinar que no hay efecto de la interacción ya que el valor de p es
mayor a 0,05, mientas que el valor de p de los factores A (tipo de verdura) y B
(temperatura) es menor a 0,05 por lo que existe efecto de los factores, por tanto se rechaza
la hipótesis nula y se comprueba la hipótesis alternativa.
C) Verifique el supuesto de igual varianza entre los tratamientos (varianza constante).
¿Se cumple satisfactoriamente?
No se cumple el supuesto de igual varianza, ya que al observar la gráfica de alimento
vs. Residuo, solo la verdura Z (A3) tiene una distribución de sus residuos estable,
siendo L y C (A1 y A2 respectivamente) las verduras con mayor variabilidad del factor
alimento.
Al analizar el gráfico de temperatura vs. Residuo, se observa que ninguno de los datos
tiene una variabilidad estable, siendo 20°C (B2) el que mayor variabilidad tiene.
D) En caso de que no se cumpla el supuesto anterior, ¿cómo afecta esto a sus
conclusiones?
La variabilidad de los residuos es un indicador de la presencia de error experimental.
Este indicador muestra que por diferentes factores no controlados, se pudo alterar los
datos obtenidos y por ende los resultados analizados.
