PROBLEMA
Una empresa embotelladora de refrescos está interesada en obtener alturas de llenado más
uniformes en las botellas que se fabrican en su proceso de manufactura. Teóricamente la máquina
de llenado llena cada botella a la altura objetivo correcto, pero en la práctica, existe variación en
torno a este objetivo, y a la embotelladora le gustaría entender mejor las fuentes de esta
variabilidad y, en última instancia, reducirla.
El ingeniero del proceso puede controlar tres variables durante el proceso de llenado: el porcentaje
de carbonatación (A), la presión de operación en el llenador (B) y las botellas producidas por
minuto o rapidez de línea (C). Es sencillo controlar la presión y la rapidez, pero el porcentaje de
carbonatación es más difícil de controlar durante la manufactura real debido a que varía con la
temperatura. Sin embargo, para los fines de un experimento, el ingeniero puede controlar la
carbonatación en tres niveles: 10,12 y 14 por ciento. Elige dos niveles para la presión (25 y 30
psi) y dos niveles para la rapidez de línea (200 y 250 bpm). El ingeniero decide correr dos replicas
con estos tres factores, haciendo las 24 corridas de manera aleatoria. En la tabla siguiente se
muestran los datos que resultaron de este experimento. Las desviaciones + son alturas de llenado
arriba del objetivo, mientras que las desviaciones – son alturas de llenado abajo del objetivo.
% de 𝒃𝟏 = 𝟐𝟓 𝒑𝒔𝒊 𝒃𝟐 = 𝟑𝟎 𝒑𝒔𝒊
carbonatación 𝒄𝟏 = 𝟐𝟎𝟎 𝒄𝟐 = 𝟐𝟓𝟎 𝒄𝟏 = 𝟐𝟎𝟎 𝒄𝟐 = 𝟐𝟓𝟎
-3 -1 -1 1
𝒂𝟏 = 𝟏𝟎
-1 0 0 1
0 2 2 6
𝒂𝟐 = 12
1 1 3 5
5 7 7 10
𝒂𝟑 = 14
4 6 9 11
a. Plantee el modelo estadístico lineal en términos del enunciado del problema.
b. Calcular el ANVA respectivo, utilice ∝ = 𝟏%. De su conclusión.
c. Determine el coeficiente de variabilidad.
d. Graficar los tres efectos principales y la interacción AB (carbonatación – presión) e
interprete el experimento.
e. Realice la prueba de Tukey para los factores y/o interacciones que tengan
significancia. Use ∝ = 𝟏% y concluya.
Solución:
a. Plantee el modelo estadístico lineal en términos del enunciado del problema.
Modelo aditivo lineal de los experimentos factoriales de tres factores
yijkl Ai B j Ck AB ij AC ik BC jk ABC ijk ijkl
Donde:
i = 1,2,…,a
j = 1,2,…,b
k = 1,2,…,c
yijkl Respuesta correspondiente al i–ésimo nivel de A, j–ésimo nivel de B, k-
ésimo nivel de C y l-ésima repetición.
Efecto de la media general.
Ai Efecto del i-ésimo nivel de A.
Bj Efecto del j-ésimo nivel de B.
Ck Efecto del k-ésimo nivel de C.
AB ij Efecto de la interacción correspondiente al i-ésimo nivel de A y j-ésimo
nivel de B.
AC ik Efecto de la interacción correspondiente al i-ésimo nivel de A y k-ésimo
nivel de C.
BC jk Efecto de la interacción correspondiente al j-ésimo nivel de B y k-ésimo
nivel de C.
ABC ijk Efecto de la interacción correspondiente al i-ésimo nivel de A, j-ésimo
nivel de B y k-ésimo nivel de C.
ijkl Error experimental correspondiente al i-ésimo nivel de A, j-ésimo nivel de
B, k-ésimo nivel de C y l-ésima repetición.
𝒂𝟏 𝒂𝟐 𝒂𝟑
𝒃𝟏 𝒃𝟐 𝒃𝟏 𝒃𝟐 𝒃𝟏 𝒃𝟐
𝒄𝟏 𝒄𝟐 𝒄𝟏 𝒄𝟐 𝒄𝟏 𝒄𝟐 𝒄𝟏 𝒄𝟐 𝒄𝟏 𝒄𝟐 𝒄𝟏 𝒄𝟐
-3 -1 -1 1 0 2 2 6 5 7 7 10
-1 0 0 1 1 1 3 5 4 6 9 11
Σ -4 -1 -1 2 1 3 5 11 9 13 16 21
Interacción AB
B
𝒃𝟏 𝒃𝟐 𝜮𝒂
A
𝒂𝟏 -5 1 -4
𝒂𝟐 4 16 20
𝒂𝟑 22 37 59
𝜮𝒃 21 54 75
Interacción AC
C
𝒄𝟏 𝒄𝟐 𝜮𝒂
A
𝒂𝟏 -5 1 -4
𝒂𝟐 6 14 20
𝒂𝟑 25 34 59
𝜮𝒄 26 49 75
Interacción BC
C
𝒄𝟏 𝒄𝟐 𝜮𝒃
B
𝒃𝟏 6 15 21
𝒃𝟐 20 34 54
𝜮𝒄 26 49 75
b. Calcular el ANVA respectivo, utilice ∝ = 𝟏%. De su conclusión.
Cálculos:
𝒚𝟐 … 752
1) TC = = = 234,375
𝒂.𝒃.𝒄.𝒓 24
(−4)2 +(20)2 +(59)2
2) 𝑺𝑪𝑨 = − 𝑇𝐶 = 252,750
𝑏.𝑐.𝑟=8
212 + 542
3) 𝑺𝑪𝑩 = − 𝑇𝐶 = 45,375
𝑎.𝑐.𝑟=12
262 +492
4) 𝑺𝑪𝑪 = − 𝑇𝐶 = 22,042
𝑎.𝑏.𝑟=12
5) 𝑺𝑪𝑨 = 𝑺𝑪𝑪𝒐𝒎𝒃 𝑨𝑩 − 𝑺𝑪𝑨 − 𝑺𝑪𝑩
(−5)2 + ⋯ + (37)2
𝑆𝐶𝐶𝑜𝑚𝑏 𝐴𝐵 = − 𝑇𝐶
𝑐. 𝑟 = 4
𝑆𝐶𝐼𝑛𝑡.𝐴𝐵 = 5,250
6) 𝑺𝑪𝑰𝒏𝒕 𝑨𝑪 = 𝑺𝑪𝑪𝒐𝒎𝒃 𝑨𝑪 − 𝑺𝑪𝑨 − 𝑺𝑪𝑪
(−5)2 + ⋯ + (34)2
𝑆𝐶𝐶𝑜𝑚𝑏 𝐴𝐶 = − 𝑇𝐶
4
𝑆𝐶𝐼𝑛𝑡.𝐴𝐶 = 0,583
7) 𝑺𝑪𝑰𝒏𝒕 𝑩𝑪 = 𝑺𝑪𝑪𝒐𝒎𝒃 𝑩𝑪 − 𝑺𝑪𝑩 − 𝑺𝑪𝑪
(6)2 + ⋯ + (34)2
𝑆𝐶𝐶𝑜𝑚𝑏 𝐵𝐶 = − 𝑇𝐶
6
𝑆𝐶𝐼𝑛𝑡.𝐵𝐶 = 1,042
𝑺𝑪𝑰𝒏𝒕 𝑨𝑩𝑪 = 𝑺𝑪𝑪𝒐𝒎𝒃 𝑨𝑩𝑪 − 𝑺𝑪𝑨 − 𝑺𝑪𝑩 − 𝑺𝑪𝑪 − 𝑺𝑪𝑰𝒏𝒕 𝑨𝑩 − 𝑺𝑪𝑰𝒏𝒕 𝑨𝑪 − 𝑺𝑪𝑰𝒏𝒕 𝑩𝑪
(−4)2 + ⋯ + (21)2
𝑺𝑪𝑪𝒐𝒎𝒃 𝑨𝑩𝑪 = − 𝑇𝐶 = 328,125
𝑟=2
𝑆𝐶𝐼𝑛𝑡.𝐴𝐵𝐶 = 1,083
8) 𝑺𝑪𝑻𝒐𝒕𝒂𝒍 = (−3)2 + (−1)2 + ⋯ + (11)2 – TC
𝑺𝑪𝑻𝒐𝒕𝒂𝒍 = 571 - 234,375
𝑺𝑪𝑻𝒐𝒕𝒂𝒍 = 336,625
9) 𝑺𝑪𝒆𝒓𝒓𝒐𝒓 = 𝑺𝑪𝒕𝒐𝒕𝒂𝒍 − 𝑺𝑪𝑨 − 𝑺𝑪𝑩 − 𝑺𝑪𝑨𝑩 − 𝑺𝑪𝑨𝑪 − 𝑺𝑪𝑩𝑪 − 𝑺𝑪𝑨𝑩𝑪
𝑆𝐶𝐶𝑜𝑚𝑏 𝐴𝐵𝐶
𝑺𝑪𝒆𝒓𝒓𝒐𝒓 = 336,625 − 328,125
𝑺𝑪𝒆𝒓𝒓𝒐𝒓 = 8,5
ANVA
F.V G.L. S.C. C.M. Fc Sig.
%Carb.(A) (a-1) = 2 252,750 126,375 178,420 **
Presión (B) (b-1) = 1 45,375 45,375 64,062 **
Rapidez (C) (c-1) = 1 22,042 22,042 31,119 **
Int AB (a-1)( b-1)=2 5,250 2,625 3,706 N.S.
Int AC (a-1)( c-1) =2 0,583 0,2915 0,412 N.S.
Int BC (b-1)(c-1) = 1 1,042 1,042 1,471 N.S.
Int ABC (a-1)(b-1)(c-1) = 2 1,083 0,5415 0,765 N.S.
Error exp. abc (r-1) =12 8,500 0,7083
Total abcr-1 =23 336,625
𝑭𝒕 𝐹𝑡 (2, 12, 0,01) = 6,93
𝐹𝑡 (1, 12, 0,01) = 9,33
Conclusión: el porcentaje de carbonatación, la presión de operación y la rapidez de línea afecta
significativamente el volumen de llenado.
c. Determine el coeficiente de variabilidad.
√𝑪𝑴𝒆𝒆
C.V. = ∗ 𝟏𝟎𝟎%
𝐘
75
Ῡ= = 3,125
24
0,7083
C.V. = √ ∗ 100
3,125
C.V. = 26,93%
d. Graficar los tres efectos principales y la interacción AB (carbonatación – presión) e
interprete el experimento.
10
DPLL
0
0 5 10 15
-2
% de carbonatación A
El número se da por el número de repeticiones
-4/8 = -0,5
20/8 = 2,5
59/8 = 7,375
10
DPLL
8
0
0 10 20 30 40
-2
Presión (B)
N° de repeticiones 12
21/12 =1,75
54/12 = 4,5
6
DPLL
5
4
3
2
1
0
0 50 100 150 200 250 300
RAPIDEZ DE LINEA
N° de repeticiones 12
26/12 = 2,17
De la tabla: 25 PSI y 30 PSI
−5 1
a1 = = −1,25 a1 = = 0,25
4 4
4 16
a2 = =1 a2 = =4
4 4
21 16+21
a3 = = 5,25 a3 = = 9,25
4 4
- Las tres variables tienen efectos principales positivos, es decir; el incremento de la
variable mueve hacia arriba la desviación promedio de llenado objetivo.
- La interacción entre la carbonatación y la presión es bastante pequeña como lo indica
las curvas de la figura.
- Como la empresa quiere que la desviación promedio de llenado objetivo esté cerca
de cero, el ingeniero decide recomendar el nivel bajo de presión de 30 psi y el nivel
alto de rapidez de línea (250 bpm) que maximizará la rapidez de producción.
e. Realice la prueba de Tukey para los factores y/o interacciones que tengan significancia.
Use ∝ = 𝟏% y concluya.
Tukey para porcentaje de carbonatación:
ALST = 𝐴𝐸𝑆𝑇 ∗ 𝑆𝑑
𝑪𝑴𝒆𝒆
(Glee = 12, ∝= 5%, 𝑝 = 3) √ = 0,2976
𝐛∗𝐜∗𝐫
AEST = 5,04
0,7083
ALS = 5,04*√ = 1,4997
8
Sacando promedios y ordenando:
ā3 = 59/8 = 7,375 a
ā2 = 20/8 = 2,500 b
ā1 = -4/8 = -0,500 c
/ā3 -ā2 / = 4,875 > 1,4997 *
/ā3 -ā1 / = 7,875 > 1,4997 *
/ā2 -ā1 / = 2,5-(-0,5) = 3 > 1,4997 *
Conclusión:
A un nivel de significancia de 0,01 existe suficiente evidencia estadística para rechazar
HO.
Por lo tanto, podemos afirmar que el nivel a3 supera a los demás tratamientos, según el
estudio realizado, la presión más baja es la más precisa.
