Diseño Experimental: Análisis Factorial
Diseño Experimental: Análisis Factorial
CHIMBORAZO.
FACULTAD DE MECÁNICA
ESCUELA DE INGENIERÍA INDUSTRIAL
DISEÑO EXPERIMENTAL
TEMA:
DISEÑO FACTORIAL
NOMBRE:
• Cristhian Ilbay
• Diego Loor
• Alex Ilbis
• Alex López
• Segundo Macas
DOCENTE:
Ing. Iván Acosta
FECHA DE ENTREGA:
2 021 - 02 - 02
PERÍODO LECTIVO:
PERIODO ACADÉMICO ESPECIAL OCTUBRE 2020 – MARZO
2021
Ejercicio #1:
HIPÓTESIS:
MODELO ESTADÍSTICO:
CATALIZADOR
B1 B2 B3 𝑌𝑖..
93 92 95
92 90 94
93 92 94
92 94 94
91 91 97
A1 2777
90 91 96
90 90 94
90 92 95
91 91 94
91 92 96
MOLDE
88 90 91
88 89 97
87 88 91
88 88 90
87 90 89
A2 2673
88 89 91
87 88 92
87 89 90
87 88 90
87 88 91
𝑌.𝑗. 1787 1802 1861 5450
𝑌𝑖𝑗2 495466 495041,667
FV SM GL CM F Valor crítico CD
MOLDE 180,27 1,00 180,27 111,12 4,02 S
CATALIZADOR 153,03 2,00 76,52 47,17 3,17 S
INTERACCIÓN 3,43 2,00 1,72 1,06 3,17 NS
ERROR 87,60 54,00 1,62
TOTAL 424,3 59
CONCLUSIÓN: Al determinar la tabla ANOVA, observamos que existen dos factores
activos o significativos, los cuales son el factor del tipo del molde y el factor del tipo de
catalizador, dado que 112,12>4,02 y 47,17>3,17 respectivamente. Entonces no podemos
rechazar la hipótesis alternativa la cual nos dice que, el molde influye en la hinchazón del
catalizador y a su vez también aceptamos la otra hipótesis alternativa la cual nos
establecía que el catalizador influye en el hinchamiento del catalizador.
c. Dibuje las gráficas de medias para los dos efectos principales con los métodos LSD
y de Tukey. Compare los resultados de ambos métodos.
𝐶𝑀𝐸 1,62
𝑇𝐵 = 𝑞(𝛼,𝑘,𝑔𝑙) *√ = 4,646 ∗ √ = 1,192
𝑛 20
CONCLUSIÓN: Mediante las gráficas LSD tanto para el molde como para el
catalizador, se observó lo que se concluyó en el ANOVA, es decir, que existe diferencia
significativa entre ellos. Por parte de la gráfica LSD del Molde observamos que no existe
traslape, y en la gráfica del LSD del catalizador, observamos que tan solo existe traslape
en el catalizador B1 Y B2.
PRUEBA DE NORMALIDAD:
CONCLUSIÓN: Dado que el valor p es menor que nuestro nivel de significancia de 0,05
podemos decir que los datos no se distribuyen de forma normal, por lo que se tendría que
realizar un ajuste en los datos para que tenga mayor veracidad en las conclusiones
obtenidas.
PRUEBA DE HOMOCEDASTICIDAD:
CONCLUSIÓN: Dado que mediante la prueba de Levene nos resultó un valor p de 0,55
lo cual es mayor a nuestro nivel de significancia de 0,05 nos indica que si cumplen con el
supuesto de homogeneidad de varianza.
f. Utilice la gráfica de residuos contra factores para detectar posibles efectos sobre
la dispersión del hinchamiento. ¿En cuál molde parece que es menor la dispersión?
CONCLUSIÓN: Una vez realizada la gráfica de residuos vs los efectos, se observa que
existe una mayor dispersión en molde de tipo A1 y en el caso de los catalizadores se
observa una semejanza entre el B1 y B2 mientras que el que posee mayor dispersión
parece ser el B3.
EJERCICIO #2.
MODELO ESTADÍSTICO:
CURADO
PEGAMENTO
60 80 100 𝑌𝑖..
2,5 3,8 4
PEGAMENTO A1 20,7
2,8 3,4 4,2
1,6 3,2 4,3
PEGAMENTO A2 17,82
1,22 2,8 4,7
𝑌.𝑗. 8,12 13,2 17,2 38,52
𝑌𝑖𝑗2 136,4384 123,6492
CURADO
PEGAMENTO 60 80 100
PEGAMENTO A1 5,3 7,2 8,2
PEGAMENTO A2 2,82 6 9
Cálculo de la suma de cuadrados:
SCT 12,79
SCA 0,69
SCB 10,35
SCAxB 1,37
SCE 0,38
ANOVA:
FV SM GL CM F Valor crítico CD
PEGAMENTO 0,69 1,00 0,69 10,99 5,99 S
CURADO 10,35 2,00 5,18 82,35 5,14 S
INTERACCIÓN 1,37 2,00 0,68 10,87 5,14 S
ERROR 0,38 6,00 0,06
TOTAL 12,8 11
En este caso como queremos elevar la resistencia a la torsión, para ello el mejor
tratamiento se observa que es mediante el tipo de PEGAMENTO 2 y una temperatura de
CURADO de 100, por lo que debemos de tomar en cuenta estos parámetros para aumentar
la resistencia a la torsión.
e. Verifique residuos.
EJERCICIO #3.
Tiempo
Producto
4 horas 8 horas 12 horas
Control 3,84 3,72 27,63 27,58 39,95 39
A 4 3,91 22 21,83 46,2 45,6
B 3,61 3,61 21,94 21,85 46,58 42,98
C 3,57 3,5 20,5 20,32 45,14 44,89
D 3,64 3,61 20,3 20,19 44,36 44,02
a) Formule el modelo estadístico más apropiado al problema y las hipótesis
estadísticas que se pueden probar.
HIPÓTESIS
Para el factor (ANTIOXIDANTES)
Ho: La relación que existe entre los antioxidantes y el tiempo es nula.; α𝛽𝑖𝑗 =
0 , 𝑝𝑎𝑟𝑎 𝑡𝑜𝑑𝑜 𝑖𝑗
HA: La relación que existe entre La relación que existe entre los antioxidantes y el tiempo
es significativa.; α𝛽𝑖𝑗 ≠ 0, 𝑎𝑙 𝑚𝑒𝑛𝑜𝑠 𝑢𝑛𝑜
MODELO MATEMÁTICO
Tiempo
Producto Suma Yi..^2
4 horas 8 horas 12 horas
Control 3,84 3,72 27,63 27,58 39,95 39
A 4 3,91 22 21,83 46,2 45,6 143,54 20603,73
B 3,61 3,61 21,94 21,85 46,58 42,98 140,57 19759,92
C 3,57 3,5 20,5 20,32 45,14 44,89 137,92 19021,93
D 3,64 3,61 20,3 20,19 44,36 44,02 136,12 18528,65
Suma 29,45 168,93 359,77 558,15 77914,24
Y.j.^2 867,30 28537,34 129434,45 158839,10
∑Y.j. ^2 19869,92
Tiempo
Producto ∑Yij^2
4 horas 8 horas 12 horas
Control 7,56 55,21 78,95
A 7,91 43,83 91,8
B 7,22 43,79 89,56 39726,25
C 7,07 40,82 90,03
D 7,25 40,49 88,38
ANOVA
F-
SC GL CM FO CD
FV CRÍTICO
SCantioxidante (A) 5,23 3 1,74 3,08 3,490 NS
SCtiempo(B) 6874,41 2 3437,21 6068,56 3,885 S
SCantioxidante X
3,01 6 0,50 0,88 2,996 NS
tiempo
SCErrror 6,80 12 0,57
TOTAL 6889,44 23 299,54
CONCLUSIÓN
En estadístico de prueba FO del tiempo (3.08) > F- critico (3,490); esto quiere decir que
el tiempo de exposición si influye en el retraso de la oxidación del aceite vegetal, por lo
tanto, se acepta HA y se rechaza Ho. Mientras que para el tipo de producto antioxidante
utilizado y la relación entre producto vs tiempo ambas resultaron ser no significativas,
por lo que se acepta Ho y se rechaza HA de manera que estos factores no influyen en el
índice de los peróxidos en otras palabras no tienen una influencia considerable reduciendo
el tiempo de oxidación.
b) ¿Los supuestos del modelo se cumplen?
En este caso la gráfica de efectos con respecto al producto se puede comprobar que no
hay significancia ya que la gráfica se asemeja a una constante, mientras que para el tiempo
se puede observar que mientras va aumentando el tiempo, el retraso de la oxidación es
mayor. Por lo que se puede concluir que la mejor opción para tener un mayor retardo de
oxidación seria exponer el producto a un determinado antioxidante a 12 horas.
40
Desemp o
30
eñ
21,915 21,895
20 20,41 20,245
10
CONCLUSIÓN
En este caso para cumplir con el propósito que es seleccionar el producto que retrase más
la oxidación, el producto ideal será el A con un tiempo de 12 horas, ya que se puede
observar en las gráficas que el valor promedio de A a 12 horas de exposición es de 45, 9,
siendo el mejor resultado de todo el tratamiento y la mejor opción para cumplir con los
requerimientos de oxidación.
EJERCICIO #4.
En la siguiente tabla se presenta la vida útil (en horas) observada sobre un tipo de batería.
Se considera que los factores que influyen son el tipo de material y la temperatura a la
que está expuesta:
TIPO DE TEMPERATURA
MATERIAL 15 70 125
1 130 155 174 180 34 40 80 75 50 70 82 58
2 150 188 159 126 136 122 106 115 45 72 60 48
3 138 110 168 160 174 120 150 139 96 104 82 78
Ho: El tipo de material con el que es construida la batería no influye en la vida útil de la
misma.; 𝑯𝒐: 𝛼1 = 𝛼2 = 𝛼3
HA: El tipo de material con el que es construida la batería si influye en la vida útil de la
misma.; 𝑯𝑨 : 𝛼1 ≠ 𝛼2 ≠ 𝛼3
Ho: La relación que existe entre el tipo de material y la temperatura a la que está
expuesta la batería es nula.; 𝑯𝒐: α𝛽𝑖𝑗 = 0 , 𝑝𝑎𝑟𝑎 𝑡𝑜𝑑𝑜 𝑖𝑗
HA: La relación que existe entre el tipo de material y la temperatura a la que está
expuesta la batería es significativa.; 𝑯𝑨 : α𝛽𝑖𝑗 ≠ 0, 𝑎𝑙 𝑚𝑒𝑛𝑜𝑠 𝑢𝑛𝑜
Modelo matemático
TIPO DE TEMPERATURA
Suma ∑Yi..^2
MATERIAL 15 70 125
1 130 155 174 180 34 40 80 75 50 70 82 58 1128 1272384
2 150 188 159 126 136 122 106 115 45 72 60 48 1327 1760929
3 138 110 168 160 174 120 150 139 96 104 82 78 1519 2307361
Suma 1838 1291 845 3974 5340674
Y.j.^2 3378244 1666681 714025 5758950
∑Y.j.^2 510130
TIPO DE TEMPERATURA
MATERIAL 15 70 125
1 639 229 260
2 623 479 225
3 576 583 360 1997822
N=36
A=3
B=3
AxB=4
Factor corrección 438685,44
SCtotal 71444,56
SCmaterial (A) 6370,72
SCtemperatura(B) 41227,06
SCmaterial X temperatura 13172,28
SCErrror 10674,50
ANOVA
F-
FV SC GL CM FO CD
CRÍTICO
SCmaterial (A) 6370,72 2 3185,36 8,057 3,354 SIGNIFICATIVO
SCtemperatura(B) 41227,06 2 20613,53 52,14 3,354 SIGNIFICATIVO
SCmaterial X
13172,28 4 3293,07 8,329 2,728 SIGNIFICATIVO
temperatura
SCErrror 10674,5 27 395,35
SCtotal 71444,56 35 2041,27
CONCLUSIÓN
CONCLUSIÓN
En este caso se puede observar en la gráfica de efectos del tipo de material que cuando la
batería es construida con el tipo de material 1 el nivel de la vida útil es bajo, mientras que
si se le construye con el tipo de material 3 la vida útil aumenta. En la gráfica de efectos
de la temperatura a la que está expuesto la batería se puede concluir que a una temperatura
de 15 °C la vida útil va a ser la máxima, ya que si la batería está expuesto a una
temperatura de 125 °C su vida útil va a ser muy corta.
CONCLUSIÓN
En este caso se puede observar a través de la gráfica de interacción del tipo de material
vs la temperatura, para que la batería para que tenga una vida útil al máximo se tendrá
que construir con el tipo de material 1 y la temperatura más adecuada a la que tendrá que
estar expuesta es 15 °C.
175
159,75
155,75
150 145,75
144
Desemp o
eñ
125
119,75
100
90
75
65
57,25 56,25
50
En este caso se puede observar que existe traslape en algunos puntos, Tipo 1 con una
temperatura de 15 °C, Tipo 2 con una temperatura de 15 °C y Tipo 3 a una temperatura
de 15 y 70 °C; mediante los valores promedio se puede concluir que la opción más
adecuada para establecer la vida útil de la batería al máximo es construir la batería con
tipo de material 1 y exponer la misma a una temperatura máxima de 15 °C ( valor
promedio = 159,75), ya que si se pasa de ese valor la vida útil va a ser menor.
𝐶𝑀𝑇 𝐶𝑀𝐸
𝑹𝟐 = ∗ 100%
𝐶𝑀𝑇
2041,27 395,35
𝑹𝟐 = ∗ 100% = 𝟖𝟎, 𝟔𝟑%
2041,27
EJERCICIO #6.
Réplica
A B I II III
- - 82 80 84
+ - 78 82 79
- + 71 70 66
+ + 89 88 93
De acuerdo a los datos de la tabla suministrados por el ejercicio, podemos afirmar que se
trata de un diseño factorial 22 ya que tienen dos factores A Y B y dos niveles en cada uno
(+; -)
Las corridas debieron ser tomadas de manera aleatoria para equilibrar el efecto de las
condiciones externas o no controlables que siempre existen en los diferentes procesos, de
manera que los datos puedan cumplir con el supuesto de independencia de manera que
los resultados obtenidos van ser mucho más fiables.
Réplica
TOTAL YATES
A B A*B I II III
- - + 82 80 84 246 (1)
+ - - 78 82 79 239 a
- + - 71 70 66 207 b
+ + + 89 88 93 270 ab
MODELO ESTADÍSTICO
Hipótesis para A
9622
𝑺𝑪𝑻 = 77840 = 719,667
3∗4
562
𝑺𝑪𝑨 = = 261
3∗4
82
𝑺𝑪𝑩 = = 5,33
3∗4
702
𝑺𝑪𝑨𝑩 = = 408,33
3∗4
𝑺𝑪𝑬 = 𝑆𝐶𝑇 𝑆𝐶𝐴 𝑆𝐶𝐵 𝑆𝐶𝐴𝐵
𝑺𝑪𝑬 = 719,667 261 5,33 408,33 = 45,007
ANOVA
Criterio de
ANOVA
decisión
FV SC GL CM Fo Fα (95%) Fα (99%) 95% 99%
Se identifica que hay menos variabilidad en el factor B y que se obtiene el mayor valor
en A+ con B+ y el menor valor con A- Y B- Y que nuestra mejor opción en la opción A
Resultado = 80,167 - 4,667 A_- + 4,667 A_+ + 0,667 B_- - 0,667 B_+ + 5,833 A*B_- -
- 5,833 A*B_- + - 5,833 A*B_+ - + 5,833 A*B_+ +
EJERCICIO 7
En una empresa lechera se han tenido problemas con la viscosidad de cierta bebida de
chocolate. Se cree que con dos ingredientes que se agregan en pequeñas cantidades se
puede resolver este problema, por lo que es necesario explorar la situación; para ello se
corre un experimento con cuatro réplicas. A continuación, se aprecian los resultados
obtenidos:
Sumatoria de cuadrados.
Contrastes
Contraste A 10,1
Contraste B 24,5
Contraste
A*B 2,5
Efectos
Efecto A 1,2625
Efecto B 17,7125
Efecto A*B 0,3125
ANOVA
ANOVA
FV SC GL CM Fo Foo CD
A 6,375625 1 6,375625 15,22 4,747 S
B 37,515625 1 37,515625 89,55 S
INTERACCION A*B 0,390625 1 0,390625 0,93 NS
ERROR 5,03 12 0,42
TOTAL 49,31 15 3,29 3,29
El l 87% de los efectos se explica por los tres factores que estamos estudiando, y el 13%
son errores naturales.
Ejercicio 8
Y2 = tiempo real (en segundos) para tomar y "colocar" los 500 componentes.
Total 8900 7
• En cambio podemos observar que el tiempo para tomar y colocar se ven afectados
por la velocidad de la mesa y la velocidad del cam, pero sus valores tambien se
encuentran dentro de un promedio alto dandonos a denotar que no varian mucho
entre ellos.
EJERCICIO #9
En una fábrica de dientes se tiene problemas con la calidad: porosidad (burbujas de aire
dentro de los dientes), manchas blancas, dientes sucios, dientes quebrados. En los intentos
por resolver los problemas han hecho cambios en algunos factores o variables del proceso.
Con base en la metodología del DDE se decide correr un diseño de experimentos. Los
factores y niveles son: temperatura de prensado (90, 130 °C), tiempo de prensado (8 y 15
minutos) y tamaño de partícula (sin tamizar y con tamizado), la variable de respuesta fue
porcentaje de diente bueno en cada corrida (un lote de producción). Los datos son los
siguientes:
Temperatura (A) Tiempo (B) T. de partícula (C) Porcentaje de diente bueno
90 8 Sin Tamizar 76.4 76.9
130 8 Sin Tamizar 76.3 76.9
90 15 Sin Tamizar 80.4 81
130 15 Sin Tamizar 77.9 79.6
90 8 Con Tamizar 84.4 84.6
130 8 Con Tamizar 84.7 84.5
90 15 Con Tamizar 82.7 83.2
130 15 Con Tamizar 85 84.7
Como se observa los efectos mas significativos resultaron ser por parte del tiempo de
prensado, del tamaño de la partícula, de la interacción doble entre AxC y la interacción
triple entre AxBxC. De manera que se observa la influencia que tiene el facto C (tamaño
de particula), en la variable respuesta (% de diente bueno).
b. Realice un análisis de varianza de estos datos y obtenga conclusiones generales.
HIPÓTESIS
Para el factor A (Temperatura de prensado)
Ho: La relación que existe entre la temperatura y tiempo de prensado es nula por lo tanto
no tiene impacto en el porcentaje de diente bueno por corrida.; α𝛽𝑖𝑗 = 0 , 𝑝𝑎𝑟𝑎 𝑡𝑜𝑑𝑜 𝑖𝑗
HA: La relación que existe entre la temperatura y tiempo de prensado tiene impacto en
el porcentaje de diente bueno por corrida.; α𝛽𝑖𝑗 ≠ 0, 𝑎𝑙 𝑚𝑒𝑛𝑜𝑠 𝑢𝑛𝑜
Ho: La relación que existe entre el tiempo de prensado y el tamaño de partícula es nula
por lo tanto no tiene impacto en el porcentaje de diente bueno por corrida.; 𝛽𝛾𝑗𝑘 =
0 , 𝑝𝑎𝑟𝑎 𝑡𝑜𝑑𝑜 𝑖𝑗
HA: La relación que existe entre el tiempo de prensado y el tamaño de partícula tiene
impacto en el porcentaje de diente bueno por corrida.; 𝛽𝛾𝑗𝑘 ≠ 0, 𝑎𝑙 𝑚𝑒𝑛𝑜𝑠 𝑢𝑛𝑜
MODELO MATEMÁTICO
FV SC GL CM Fo Fcritico CD
A 0 1 0 0 NS
B 6.0025 1 6.0025 22.4392523 S
C 146.41 1 146.41 547.327103 S
AxB 0.0025 1 0.0025 0.00934579 5.31765507 NS
BxC 14.0625 1 14.0625 52.5700935 S
AxC 4 1 4 14.953271 S
AxBxC 3.4225 1 3.4225 12.7943925 S
Error 2.14 8 0.2675
TOTAL 176.04 15 11.736
R^2 98%
CONCLUSIÓN
Una vez analizado el ANOVA se observa que los efectos producidos por los factores
principales B (tiempo) y C (tamaño de particula), las interacciones dobles entre AxC y
BxC así como el efecto de la interacción triple AxBxC, son significativos con un 98% de
manera que para estos se acepta la hipótesis alternativa y se rechaza la nula; pues si
producen un efecto sobre la variable Y (porcentaje de diente bueno).