Capitulo 6
Capitulo 6
Capitulo 6
FACTORIAL 2 k
CAPITULO 6
Diseo y Anlisis de Experimentos
Montgomery
DISEO FACTORIAL 2k
DISEO FACTORIAL 2k (cont.)
i 1 j 1 k 1
SSE=SST-SSA-SSB-SSAB=323-208.33-75-8.33=31.34
ANALISIS DE VARIANZA
Promedio
Fuente de Suma de Grados de
de los F0 Valor P
variacin cuadrados libertad
cuadrados
A 208.33 1 208.33 53.18 0.0001
B 75.00 1 75.00 19.14 0.0024
AB 8.33 1 8.33 2.13 0.1829
Error 31.34 8 3.92
Total 323.00 11
EL MODELO DE REGRESION
Para el experimento del proceso qumico
que se viene tratando, se tiene:
y = b0 + b1x1 + b2x2 + e
Donde:
x1 es la variable codificada que
representa el reactivo.
x2 es la variable codificada que
representa la cantidad del catalizador.
Las b son los coeficientes de regresin.
EL MODELO DE REGRESION (cont.)
La relacin entre las variables naturales la
concentracin del reactivo y la cantidad de
catalizador y las variables codificadas es:
Catalizador + Catalizador / 2
alta baja
Catalizado r
x Catalizador Catalizador / 2
bajo alto
2
alto bajo
8
RESIDUALES Y ADECUACION DEL
MODELO
El modelo de regresin puede usarse para obtener el valor
predicho o ajustado de y en los cuatro puntos del modelo.
Los residuales son la diferencia entre el valor observado y el
valor ajustado de y.
Cuando la concentracin del reactivo est en el nivel bajo
(x1 = -1) y el catalizador est en el nivel bajo (x2 = -1), el
rendimiento predicho es:
^
8.33 1 + 5.00 1 25.835
y 27.5+
2 2
Los residuales para las tres observaciones en esta
combinacin, son:
e1 = 28 25.835 = 2.165
e2= 25 25.835 = -0.835
e3 = 27 25.835 = 1.165
RESIDUALES Y ADECUACION DEL
MODELO (continuacin)
Para el nivel alto de la concentracin del reactivo y el
nivel bajo del catalizador se tiene:
^
8.33 + 1 + 5.00 1 34.165
y 27.5+
2
2
Los residuales para las tres observaciones en esta
combinacin, son:
e1 = 36 34.165 = 1.835
e2 = 32 34.165 = -2.165
e3 = 32 34.165 = -2.165
RESIDUALES Y ADECUACION
DEL MODELO (continuacin)
Para el nivel bajo de la concentracin del reactivo y el
nivel alto del catalizador se tiene:
^
8.33 1 + 5.00 + 1 20.835
y 27.5+
2
2
Los residuales para las tres observaciones en esta
combinacin, son:
e1 = 18 20.835 = -2.835
e2 = 19 20.835 = -1.835
e3 = 23 20.835 = 2.165
RESIDUALES Y ADECUACION
DEL MODELO (continuacin)
Para el nivel alto de la concentracin del reactivo y el
nivel alto del catalizador se tiene:
^
8.33 + 1 + 5.00 + 1 29.165
y 27.5+
2
2
Los residuales para las tres observaciones en esta
combinacin, son:
e1 = 31 29.165 = 1.835
e2 = 30 29.165 = 0.835
e3 = 29 29.165 = -0.165
GRAFICA DE PROBABILIDAD
NORMAL DE LOS RESIDUALES
Serie Residual % Prob Norm GRAFICA DE PROBABILIDAD NORMAL
2 -2.165 12.50
3 -2.165 20.83 100.00
4 -1.835 29.17
% PROBABILIDAD NORMAL
80.00
5 -1.835 37.50
6 -0.835 45.83 60.00
7 -0.165 54.17
40.00
8 0.835 62.50
9 1.165 70.83 20.00
10 1.835 79.17
0.00
11 2.165 87.50 -4 -3 -2 -1 0 1 2 3
RESIDUAL
12 2.165 95.83
SUPERFICIE DE RESPUESTA
SUPERFICIE DE RESPUESTA (cont.)
SUPERFICIE DE RESPUESTA (cont.)
SUPERFICIE DE RESPUESTA (cont.)
EL DISEO 23
Suponga tres factores A, B y C, cada uno con dos niveles.
EL DISEO 23 - NOTACIONES
Existen tres notaciones diferentes para las
corridas de diseo 23 que son de uso
general.
La primera es la notacin + y -, llamada con
frecuencia notacin geomtrica.
La segunda es el uso de las etiquetas en
letras minsculas para identificar las
combinaciones de los tratamientos.
La tercera notacin utiliza 1 y 0 para denotar
los niveles alto y bajo, respectivamente, de
los factores, en lugar de + y -.
EL DISEO 23 (continuacin)
Las combinaciones de los tratamientos en el orden estndar
se escriben como: (1), a, b, ab, c, ac, bc, abc.
Las diferentes notaciones para el diseo 23 se ilustran as:
Corrida A B C Etiquetas A B C
1 - - - (1) 0 0 0
2 + - - a 1 0 0
3 - + - b 0 1 0
4 + + - ab 1 1 0
5 - - + c 0 0 1
6 + - + ac 1 0 1
7 - + + bc 0 1 1
8 + + + abc 1 1 1
EL DISEO 23 - GRADOS DE LIBERTAD
Hay siete grados de libertad entre las
ocho combinaciones de tratamientos del
diseo 23, as:
Tres grados de libertad se asocian con los
efectos principales de A, B, y C.
Cuatro grados de libertad se asocian con las
interacciones; uno con cada una de la
interacciones AB, AC, BC y uno con la
interaccin ABC.
EL DISEO 23 EFECTO de A
El efecto de A en el nivel bajo de B y C es: [a-(1)]/n
El efecto de A en el nivel alto de B y en nivel bajo de
C es: [ab-b]/n
El efecto de A en el nivel alto de C y en nivel bajo de
B es: [ac-c]/n
El efecto de A en el nivel alto de B y C es: [abc-bc]/n
El efecto principal promedio de A es:
A = [a + ab + ac + abc - (1) b c - bc] / 4n
EL DISEO 23 EFECTO de B
El efecto de B en el nivel bajo de A y C es: [b-(1)]/n
El efecto de B en el nivel alto de A y en nivel bajo de
C es: [ab-a]/n
El efecto de B en el nivel alto de C y en nivel bajo de
A es: [bc-c]/n
El efecto de B en el nivel alto de B y C es: [abc-ac]/n
El efecto principal promedio de B es:
B = [b + ab + bc + abc - (1 ) a c - ac] / 4n
EL DISEO 23 EFECTO de C
El efecto de C en el nivel bajo de A y C es: [c-(1)]/n
El efecto de C en el nivel alto de A y en nivel bajo de
B es: [ac-a]/n
El efecto de C en el nivel alto de B y en nivel bajo de
A es: [bc-b]/n
El efecto de C en el nivel alto de A y B es: [abc-ab]/n
El efecto principal promedio de C es:
C = [c + ac + bc + abc - (1) a b - ab] / 4n
EL DISEO 23 EFECTO INTERACCION AB
Una medida de la interaccin AB es la
diferencia entre los efectos promedio de A
con los dos niveles de B:
B Efecto promedio de A
Alto (+) [(abc bc) + (ab b)] / 2n
Bajo (-) {(ac c) + [a (1)]} / 2n
Diferencia [abc+ab+c+(1)acbcab] / 4n
EL DISEO 23 EFECTO INTERACCIONES
AB = [abc+ab+c+(1)acbcab] / 4n
AC = [abc+ac+b+(1)abbcac] / 4n
BC = [abc+bc+a+(1)abacbc] / 4n
ABC = [abc+a+b+cabacbc(1)] / 4n
EL DISEO 23 EFECTOS (Otro mtodo de obtenerlos)
bc abc
c ac
b ab
(1) a
Signos algebraicos para calcular
los efectos del diseo 23
Combinacin de Efecto Factorial
Tratamientos I A B AB C AC BC ABC
(1) + - - + - + + -
a + + - - - - + +
b + - + - - + - +
ab + + + + - - - -
c + - - + + - - +
ac + + - - + + - -
bc + - + - + - + -
abc + + + + + + + +
Ejemplo 6-1 Pag. 232
Una empresa embotelladora de refrescos est
interesada en obtener alturas de llenado ms uniforme
en las botellas que se fabrican en su proceso de
manufactura. Tericamente, la mquina de llenado
llena cada botella a la altura objetivo correcta, pero en
la prctica existe variacin en torno a este objetivo y a
la embotelladora le gustara entender mejor las
fuentes de esta variabilidad y, en ltima instancia
reducirla.
El ingeniero de proceso puede controlar tres variables
durante el proceso de llenado: el porcentaje de
carbonatacin (A), la presin de operacin en el
llenado (B) y las botellas producidas por minuto
rapidez de lnea (C).
Para cada factor slo se usan dos niveles.
Los datos del experimento se muestran a
continuacin:
Ejemplo 6-1 (continuacin)
Factores Desviacin de la
Nivel del factor
Corrida Codificados altura de llenado
A B C Rplica 1 Rplica 2 Bajo (-1) Alto (+1)
1 -1 -1 -1 -3 -1 A (%) 10 12
2 1 -1 -1 0 1 B (psi) 25 30
3 -1 1 -1 -1 0 C (bpm) 200 250
4 1 1 -1 2 3
5 -1 -1 1 -1 0
6 1 -1 1 2 1
7 -1 1 1 1 1
8 1 1 1 6 5
El diseo 23 para el experimento
de la desviacin de la altura de
llenado del ejemplo 6-1
=2 =11
= -1
=3
Velocidad (C)
=- 1
=5
Presin (B)
=-4 =1
Carbonatacin (A)
Efectos de los factores
A = [a - (1) + ab b + ac c + abc - bc] / 4n
A = [ 1 (- 4) + 5 (- 1) + 3 (- 1) + 11 2] / 8 = 24 / 8 = 3.00
B = [b + ab + bc + abc - (1 ) a c - ac] / 4n
B = [- 1 + 5 + 2 + 11 (- 4) 1 (- 1) -3] / 8 = 18 / 8 = 2.25
C = [c + ac + bc + abc - (1) a b - ab] / 4n
C = [- 1 + 3 + 2 + 11 (- 4) 1 (- 1) 5] / 8 = 14 / 8 = 1.75
AB = [ab a b + (1) + abc bc ac + c] / 4n
AB = [5 1 (- 1) + (- 4) + 11 2 3 + (- 1) ] / 8 = 6 / 8 = 0.75
AC = [(1) - a + b - ab c ac bc + abc] / 4n
AC = [- 4 1 + (- 1) - 5 - (- 1) + 3 2 + 11] / 8 = 2 / 8 = 0.25
BC = [ (1) + a b ab c ac + bc + abc] / 4n
BC = [- 4 + 1 - (- 1) 5 - (- 1) 3 + 2 + 11] / 8 = 4 / 8 = 0.50
ABC = [abc bc ac + c ab + b + a (1)] / 4n
ABC = [11 2 3 + (- 1) 5 + (- 1) + 1 (- 4)] / 8 = 4 / 8 = 0.50
SUMAS DE CUADRADOS - SS
SS = (Contraste)2 / 8n
SSA = (24)2 / 16 = 36.00
SSB = (18)2 / 16 = 20.25
SSC = (14)2 / 16 = 12.25
SSAB = (6)2 / 16 = 2.25
SSAC = (2)2 / 16 = 0.25
SSBC = (4)2 / 16 = 1.00
SSABC = (4)2 / 16 = 1.00
Resumen de la estimacin de
los efectos
Estimacin Sumas de Contribucin
Factor
del efecto cuadrados porcentual
A 3.00 36.00 46.153846
B 2.25 20.25 25.961538
C 1.75 12.25 15.705128
AB 0.75 2.25 2.884615
AC 0.25 0.25 0.320513
BC 0.50 1.00 1.282051
ABC 0.50 1.00 1.282051
Error puro 5.00 6.410256
Total 78.00
Anlisis de Varianza de los datos
de la altura de llenado
120.000
PROBABILIDAD NORMAL
100.000
=((I-0.5)/N)*100
80.000
60.000
40.000
20.000
0.000
-1.5 -1 -0.5 0 0.5 1 1.5
RESIDUAL
EL DISEO GENERAL 2k
k
3. interacciones de tres factores.
3
4. ...............
k
5. k = una interaccin de k factores.
ANALISIS DE VARIANZA DE UN DISEO
GENERAL 2k (continuacin)
Suma de Grados de
Fuente de variacin
cuadrados libertad
k efectos principales
A SSA 1
B SSB 1
: : :
K SSK 1
k
interacciones de dos factores
2
AB SSAB 1
AC SSAC 1
: : :
JK SSJK 1
kk
interacciones de tres factores
3
3
ABC SSABC 1
ABD SSABD 1
: : :
IJK SSIJK 1
: : :
k
= 1 interaccin de k factores
k
ABC...K SSABC...K 1
k
Error SSE 2 (n - 1)
Total SST n2k - 1
En general, el contraste del
efecto ABK se determina as:
ContrasteABK = (a 1) (b 1) (k 1)
Para expandir esta ecuacin se usa el lgebra
ordinaria reemplazando 1 con (1) en la
expresin final.
El signo de cada grupo de parntesis es
negativo si el factor est incluido en el efecto y
es positivo si el factor no est incluido.
Ejemplo, considere un diseo factorial 23, el
contraste AB sera:
ContrasteAB = (a - 1)(b - 1)(c + 1)
= abc + ab + c + (1) ac bc a - b
Ejemplo 6-2 Pag. 246
Un producto qumico se fabrica en un envase presurizado. Se lleva
a cabo un experimento factorial en la planta piloto para estudiar los
factores que se piensan influyen en el ndice de filtracin de este
producto.
Los cuatro factores son: la temperatura (A), la presin (B), la
concentracin del formaldehdo (C) y la velocidad de agitacin (D).
Cada factor est presente con dos niveles. La matriz de diseo y los
datos de la respuesta obtenidos de una sola rplica de un diseo 24
se muestran a continuacin. Las 16 corridas se hacen de manera
aleatoria.
El ingeniero del proceso est interesado en maximizar el ndice de
filtracin. Las condiciones actuales del proceso producen ndice de
filtracin de alrededor de 75 gal/h. As mismo, en el proceso actual
la concentracin de formaldehdo, factor C, se usa en el nivel alto.
Al ingeniero le gustara reducir la concentracin de formaldehdo lo
ms posible, pero no ha podido hacerlo porque siempre produce
ndices de filtracin ms bajos.
Tabla 6-10 Experimento del ndice de filtracin en la planta piloto
Nmero Factor Etiqueta de Indice de
de A B C D Corrida Filtracin
1 - - - - (1) 45
2 + - - - a 71
3 - + - - b 48
4 + + - - ab 65
5 - - + - c 68
6 + - + - ac 60
7 - + + - bc 80
8 + + + - abc 65
9 - - - + d 43
10 + - - + ad 100
11 - + - + bd 45
12 + + - + abd 104
13 - - + + cd 75
14 + - + + acd 86
15 - + + + bcd 70
16 + + + + abcd 96
El diseo 24 para el experimento del ndice de
filtracin del ejemplo 6-2
c ac cd acd
b ab bd abd
(1) a d ad
Estimacin de los efectos de los factores y
sumas de cuadrados del diseo factorial 24
Trmino del Estimacin Con- Suma de Contribucin
Modelo del Efecto traste Cuadrados Porcentual
A 21.625 173 1870.5625 32.6397
B 3.125 25 39.0625 0.6816
AB 0.125 1 0.0625 0.0011
C 9.875 79 390.0625 6.8063
AC -18.125 -145 1314.0625 22.9293
BC 2.375 19 22.5625 0.3937
ABC 1.875 15 14.0625 0.2454
D 14.625 117 855.5625 14.9288
AD 16.625 133 1105.5625 19.2911
BD -0.375 -3 0.5625 0.0098
ABD 4.125 33 68.0625 1.1876
CD -1.125 -9 5.0625 0.0883
ACD -1.625 -13 10.5625 0.1843
BCD -2.625 -21 27.5625 0.4809
ABCD 1.375 11 7.5625 0.1320
Suma 5730.9375 100.00000
GRAFICA DE PROBABILIDAD NORMAL DE LOS
EFECTOS PARA EL DISEO FACTORIAL 24
Efecto % Prob.
-18.1 3.33 Grfica de Probabilidad Normal de los Efectos para el Diseo
-2.63 10.00 Factorial 2^4 Ejemplo 6-2
-1.63 16.67
120.00
-1.13 23.33
-0.38 30.00 100.00 A
% de Probabilidad Normal
0.125 36.67 AD
80.00 D
1.375 43.33 C
1.875 50.00 60.00
2.375 56.67
3.125 63.33 40.00