Un experimentador está estudiando el efecto de cinco posibles mezclas de material

explosivo para la fabricación de dinamita. El material proviene de cinco orígenes distintos y

es manipulado por 5 operarios para realizar las mezclas. Los datos que se recogen a la
fuerza al explotar.

Operarios
Origen 1 2 3 4 5
1 A (24) B (20) C (19) D (24) E (24)
2 B (17) C (24) D (30) E (27) A (36)
3 C (18) D (38) E (26) A (27) B (21)
4 D (26) E (31) A (26) B (23) C (22)
5 E (24) A (30) B (20) C (29) D (31)
a) Identifique los elementos del experimento. Escribir detalladamente el modelo.
La variable respuesta es la fuerza al explotar la mezcla.
Los factores a estudiar son: El origen (1, 2, 3, 4 y 5), el operario (1, 2, 3, 4 y 5) y la mezcla
(A, B, C, D y E).
Las unidades experimentales son el montaje de cada mezcla.
El modelo a plantear es un modelo estadístico lineal de la forma

{
i=1 , 2 , . .. , p
y ijk =μ+α i +τ j +β k +ε ijk j=1 , 2 , . .. , p
k=1 , 2 , .. . , p

donde
Y ijk es una variable aleatoria que denota la observación de la fuerza al explotar en la

fila i-ésima y la columna k-ésima para el tratamiento j-ésimo, μ es la media global de la

fuerza al explotar la mezcla,

α i es el efecto del origen i-ésimo, τ j es el efecto del tipo de

mezcla j-ésimo,
β k es el efecto del operador k-ésimo, y ε ijk es un componente del error
aleatorio.
b) Determine si influyen todas las variables consideradas. ¿Considera que debe modificar el
modelo?

Se deben conocer las sumas de cuadrados de los tratamientos (Mezcla), del origen, del

operador, del error, y el total, los cuales se determinan con las siguientes ecuaciones:
Total.
p
y 2.. .
SCT =∑ y 2ijk −
j =1 p2
Tratamientos.
p 2
1 y
SCTrat = ∑ y.2j.− .2..
p j=1 p (Letras)
Filas.
p 2
1 y
SC Filas= ∑ y i2..− .2. .
p i=1 p
Columnas.
p 2
1 y
SCColumnas = ∑
p k =1
y 2. . k − .2. .
p
Error.
SC E=SCT −SCTrat −SC Filas−SC Columnas

Los resultados previos a la suma de cuadrados se obtuvieron con una hoja de cálculo de

Excel y se muestran a continuación:

Totales ( y 2i . Promedios (
Origen Observaciones ri y i. ȳ i.
)
p )
24 20 19 24 24
1 5 111 2464.20 22.200
17 24 30 27 36
2 5 134 3591.20 26.800
18 38 26 27 21
3 5 130 3380.00 26.000
26 31 26 23 22
4 5 128 3276.80 25.600
24 30 20 29 31
5 5 134 3591.20 26.80

Totales (
yk . 109 143 121 130 134
)
y 2k . 2376. 2928. 3591.
4089.8 3380.0
P 2 2 2
Promedios 21.8 28.6 24.2 26.0 26.8
 ȳ k
( )

y . . .=637
ȳ . . .=25. 480

Suma de cuadrados.

Total.
t r
y . .2
SCT =∑ ∑ y 2i j k −
i=1 j=1 N
( 637 )2
SCT =242 +202 +19 2 +. ..+202 +292 +31 2−
25
405769
SCT =16897−
25
SCT =16897−16230 . 76
SCT =664 . 24

Origen (Filas).
t r
1 y . .2
SC Filas= ∑∑ y 2
−
p i=1 j=1 i. . N
(637 )2
SCOrigen =2464 . 20+3591 .20+3380. 00+3276 . 80+3591 . 20−
25
405769
SCOrigen =16303 . 40−
25
SCOrigen =16303 . 40−16230 . 76
SCOrigen =72. 64

Operador (Columnas).
p 2
1 y
SCColumnas = ∑ y 2. . k − .2. .
p k =1 p
(637)2
SCOperador =2376 . 20+4089. 80+2928 . 20+3380 . 0+3591. 2−
25
 405769
SCOperador =16365 . 40−
25
SCOperador =16365 . 40−16230 .76
SCOperador =134 . 64

Tratamientos (Letra)
p 2
1 y
SCTrat = ∑ y.2j.− .2..
p j=1 p

Total del
y 2. j.
Letra latina tratamiento,
p
y. j .
A y . A .=143 4089.8
B y . B .=101 2040.2
C y . C .=112 2508.8
D y . D .=149 4440.2
E y . E.=132 3484.8
1
p
y 2. ..
SCTrat = ∑
p j=1
y.2j.−
p 2

(637 )2
SCTratamiento =4089 . 8+2040 . 2+ 2508 .8+4440 .2+3484 .8−
25
405769
SCTratamiento =16563 . 80−
25
SCTratamiento =16563 . 80−16230 .76
SCTratamiento =333 . 04

Error.

SC E=SCT −SCTrat −SC Filas−SC Columnas

SC E=666 . 24−333 .04−72 . 64−134 . 64
SC E=125 . 92

Los grados de libertad.

Tratamiento.

g .l Tratmientos =p−1
g .l Tratmientos =5−1
g .l Tratmientos =4

Filas.

g .l Origen = p−1
g .l Origen =5−1
g .l Orien =4

Columnas.

g .l Operador = p−1
g .l Operador =5−1
g .l Operador =4

Error.

g .l Error =( p−2) ( p−1)

g .l Error =(5−2) (5−1)
g .l Error =3×4
g .l Error =12

Total.

g .l Total = p 2−1
g .l Total =(5)2 −1
g .l Total =25−1
g .l Total =24
Cuadrado medio.

El cuadrado medio se determina dividiendo la suma de cuadrados entre los grados de

libertad.

Tratamiento.

SC Tratamiento
CM Tratamiento=
p−1
333 . 04
CM Tratamiento=
4
CM Tratamiento=83 . 26

Filas.

SC Fila
CM Fila=
p−1
72 .64
CM Origen=
4
CM Origen=18 . 16

Columnas.
SC Columna
CM Columna=
p−1
134 .64
CM Operador =
4
CM Operador =33. 66

Error.
SC Error
CM E=
( p−2) ( p−1)
125 . 92
CM E=
12
CM E=10 . 49
Estadístico de prueba.

CM Tratamientos CM Filas CM Bloques

F= F= F=
CM E CM E CM E
83 . 26 18 .16 33 .66
F= F= F=
10 . 49 10 . 49 10 . 49
F=7 . 937 F=1 .731 F=3 .209

Los resultados se resumen en la tabla de análisis de varianza.

Fuente de Suma de Grados de Cuadrado

F
variación cuadrados libertad medio
Mezcla 333.04 4 83.26 7.937
Origen 72.64 4 18.16 1.731
Operador 134.64 4 33.66 3.209
Error 125.92 12 10.49
Total 24
F c=F α , p−1 , ( p−2) ( p−1 )=F 0 .05, 4 ,14 =3. 259
Siendo el valor crítico , Se rechaza H0 tanto para la
mezcla como para el operador, mientras que no se rechaza para el origen.

No todas las variables consideradas influyen en la fuerza al explotar. Es necesario

modificar el modelo considerando sólo la mezcla y el operador.
Los datos deben disponerse de la siguiente manera y realizar un análisis de bloques.
Operarios
Mezcla 1 2 3 4 5
A 24 30 26 27 36
B 17 20 20 23 21
C 18 24 19 29 22
D 26 38 30 24 31
E 24 31 26 27 24
Los resultados se resumen en la tabla siguiente:

Fuente de Suma de Grados de

Cuadrado medio F
variación cuadrados libertad
Tratamientos 333.04 4 83.2600 6.7091
 Bloques 134.64 4 33.6600 2.7123
Error 198.56 16 12.4100
Total 666.24 24
F c=F α , t−1 , ( t−1) ( b−1 )=F 0 .05 , 4 ,20=3. 112
Siendo el valor crítico , Se rechaza H0 tanto para la
mezcla más no para el operador.

Finalmente, se debe hacer un análisis de los tratamientos (Mezcla) sin considerar los
bloques.

Fuente de Suma de Grados de

Cuadrado medio F
variación cuadrados libertad
Tratamientos 333.04 4 83.2600 4.9976
Error 333.20 20 16.6600
Total 666.24 24
F c=F α , t−1 , t ( r−1)=F0 . 05 , 4 , 20=2 .866
Siendo el valor crítico , Se rechaza H0 tanto para la
mezcla, comprobándose que la fuerza al explotar, en media, no es la misma para las cinco
mezclas.
Un ingeniero diseña una batería para su uso en el motor de cilcomotores. Para ello dispone
de tres tipos diferentes de material. Como considera que la temperatura es un factor
influyente en la duración de la batería, decide diseñar el experimento combinando los 3
materiales con 3 temperaturas concretas: –10°C, 20°C y 50°C. Así, una vez sometidas las 4
baterías de cada material a cada temperatura, las duraciones efectivas de las baterías en
horas:
Material
1 2 3
-10° 20° 20°
130 122 150
-10° 20° 20°
155 126 174
50° -10° -10°
70 150 138
20° 50° -10°
34 70 110
-10° 50° 50°
80 25 60
20° 20° 50°
80 115 96
50° -10° 20°
82 126 120
20° 50° -10°
40 58 168
50° -10° -10°
20 159 160
-10° -10° 20°
74 188 139
20° 20° 50°
75 106 104
50° 50° 50°
58 45 82
a) Identifique los elementos del experimento. Escribir detalladamente el modelo
matemático.
La variable respuesta es la duración de la batería.
Los factores a estudiar son: El material (1, 2, 3) y la temperatura (-10°, 20°, 50°).
Las unidades experimentales son el montaje de cada batería.
El modelo a plantear es un modelo estadístico lineal de la forma

{
i=1 , 2 , .. . , a
y ijk =μ+τ i + β j +(τ β )i j +ε ijk j=1, 2 , . .. , b
k=1 , 2 , . . ., n
Yijk
donde es una variable aleatoria que denota la duración de la batería en la fila (tipo de
material) i-ésima y la columna (temperatura) j-ésima,  es el efecto promedio global de

duración de la batería,
τ i es el efecto del tipo de material i-ésimo, β j es el efecto de la

temperatura j-ésima, (τ β )i j es el efecto de interacción entre  i y  j , y  ijk es un

componente del error aleatorio.
b) Estudiar si el tipo de material y la temperatura son factores determinantes en la duración
de la batería. ¿Hay posibilidad de que un material sea más recomendado a una temperatura
en concreto y no lo sea a otra distinta?
Los datos se resumen en la siguiente tabla:
Tipo de Temperatura
material -10° 20° 50°
130 155 34 40 20 70
1
74 180 80 75 82 58
150 188 126 122 25 70
2
159 126 106 115 58 45
138 110 174 120 96 104
3
168 160 150 139 82 60
Se determinan los totales por tipo de material y temperatura.
Temperatura
-10° 20° 50°
130 155 34 40 20 70
1 539 229 230 998
74 180 80 75 82 58
150 188 126 122 25 70
2 623 469 198 1290
159 126 106 115 58 45
138 110 174 120 96 104
576 583 342 1501
3 168 160 150 139 82 60
1738 1281 770 3789
Las sumas de cuadrados se calculan de la siguiente manera:
Total.
 a b n
y .. .2
SCT =∑ ∑ ∑ y ijk
2
−
i=1 j=1 k =1 abn

2 2 2 2 2( 3789 )2
2
SC T =( 130 ) +(155 ) +( 34 ) +. ..+( 139) +( 82) +( 60 ) −
3×3×4
14356521
SCT =475927−
36
SCT =475927−398792 .25
SCT =77134 . 75
Material.
a
1 y . ..2
SC A = ∑ y 2
−
b n i=1 i. . a b n
1 2 2 2 (3789 )2
SC Material = [(998 ) +(1290 ) +(1501 ) ]−
3×4 3×3×4
1 14356521
SC Material = ( 4913105)−
12 36
SC Material =409425 . 42−398792. 25
SC Material =10633 .17
Temperatura.
a
1 y . ..2
SC B= ∑ y 2. j.−
a n i=1 abn
1 (3789)2
SCTemperatura= [(1738 )2 +(1281 )2 +(770)2 ]−
3×4 3×3×4
1 14356521
SCTemperatura= (5254505)−
12 36
SCTemperatura=437875 . 42−398792 .25
SCTemperatura=39083 . 17
Interacción.
a b
1 y . ..2
SC AB = ∑∑ y 2
−
n i=1 j=1 i j. a b n
−SC A −SC B
 1 (3789 )2
SC Material×Temperatura= [(539)2 +(229 )2 +. .. ..+(583)2 +(342 )2 ]− −10633 .17−39083 .17
4 3×3×4
1 14356521
SC Material×Temperatura= (1831785 )− −10633 . 17−39083 . 17
4 36
SC Material×Temperatura=457946 . 25−398792 . 25−10633 . 17−39083 . 17
SC Material×Temperatura=9437 . 66
Error.
SC E=SCT −SC A −SC B −SC AB
SC E=77134 . 75−10633. 17−39083. 17−9437 . 66
SC E=17980 . 75

Los resultados se resumen en la tabla de análisis de varianza.

Fuente de Suma de Grados de Cuadrado

F
variación cuadrados libertad medio
Material 10633.17 2 5316.585 7.983
Temperatura 39083.17 2 19541.585 29.344
Interacción 9437.66 4 2359.415 3.543
Error 17980.75 27 665.954
Total 77134.75 35
Decisiones.
Material.
F c=F α , t−1 , ( t−1) ( b−1 )=F 0 .05 , 4 ,20=3. 112
Siendo el valor crítico , Se rechaza H0
Temperatura.
F c=F α , t−1 , ( t−1) ( b−1 )=F 0 .05 , 4 ,20=3. 112
Siendo el valor crítico , Se rechaza H0
Interacción.
F c=F α , t−1 , ( t−1) ( b−1 )=F 0 .05 , 4 ,20=3. 112
Siendo el valor crítico , Se rechaza H0