DISEÑO DE EXPERIMENTOS

Diseños con bloques aleatorizados

DISEÑO DE BLOQUES COMPLETOS ALEATORIZADOS

En muchos problemas de experimentos, es necesario

hacer un diseño de tal manera que la variabilidad
proveniente de fuentes conocidas pueda ser
sistemáticamente controlada. El diseño con bloques
completos aleatorizados pretende reducir el efecto de la
variabilidad proveniente de causas propias del
experimento pero independiente del efecto que se desea
estudiar.

Para los fines del análisis de varianza el bloqueo introduce

un efecto adicional ficticio, cuyo objetivo es separar del
error experimental, alguna fuente de variabilidad
conocida.
DISEÑO DE BLOQUES COMPLETOS ALEATORIZADOS

Ejemplo:

Se desea saber si el tipo de la herramienta de medida

tiene efecto en las lecturas de dureza de cierto material.
La prueba de dureza consiste en someter a cierta presión
la herramienta sobre la muestra de metal y medir la
profundidad del orificio producido en la muestra. No se
sabe como la muestra de material podría afectar la
medida registrada, pero se sabe que las muestras no
necesariamente provienen de un material homogéneo. Para
aislar el efecto del material sobre la medida, se somete la
misma muestra a la prueba de dureza con todos los
instrumentos y se analizan los resultados. Los resultados
podemos verlos en la siguiente tabla:
DISEÑO DE BLOQUES COMPLETOS ALEATORIZADOS

RESULTADOS DE PRUEBAS DE DUREZA

PARA CUATRO MUESTRAS
(Escala Rockwell)

Material de prueba
Herramienta a b c d
1 9.3 9.4 9.6 10.0
2 9.4 9.3 9.8 9.9
3 9.2 9.4 9.5 9.7
4 9.7 9.6 10.0 10.2
DISEÑO DE BLOQUES COMPLETOS ALEATORIZADOS

TABLA ANOVA

Fuente de Suma de Grados de Cuadrados F0

Variación Cuadrados Libertad medios
Tratamiento a y
2
y
2
a-1 SS Tratamiento MS
Tratamiento

i 1 b
i.

N
..
a 1 MS E

2
Bloque b y y
2
b-1 SS Bloque

 
.j ..

a N b 1
j 1

Error SS (Por E
(a - 1)(b - SS E

1) (a  1)(b  1)
diferencia)
2
Total a b y N-1
 y
2
 ..

i 1 j 1
ij N
DISEÑO DE BLOQUES COMPLETOS ALEATORIZADOS

Para simplificar restemos 9.5 de cada una de las observaciones y

multipliquemos por 10 para formar la tabla siguiente:

RESULTADOS DE PRUEBAS DE DUREZA

PARA CUATRO MUESTRAS
(Escala Rockwell)

Material de Prueba
Herramienta a b c d y i.

1 -2 -1 1 5 3
2 -1 -2 3 4 4
3 -3 -1 0 2 -2
4 2 1 5 7 15
y  20
y .j -4 -3 9 18 ..
DISEÑO DE BLOQUES COMPLETOS ALEATORIZADOS

Cálculo de resultados

 4  2  15 20
2 2 2 2 2

SS  3   38.50
Tratamiento
4 16

4  3  9  18 20
2 2 2 2 2

SS Bloque  4

16
 82.50

SS T
 154  20  129.00
16

SS  SS  SS
E T tratamiento
 SSbloque  129.00  3850
.  82.50  8
DISEÑO DE BLOQUES COMPLETOS ALEATORIZADOS

TABLA ANOVA

Fuente de Suma de Grados de Cuadrados

Variación Cuadrados Libertad Medios F0

Tratamiento 38.50 3 12.83 14.44

Bloque 82.50 3 27.50
Error 8.00 9 0.89
Total 129.00 15

Buscando en la tabla F para un nivel de significación de 5% y 3 y 9

grados de libertad, tenemos que F 0.05,3,9 = 3.86, dado que 14.44 >
3.86, concluimos en que el tipo de herramienta si afecta la lectura de
la dureza del material.
Ejercicio 4.1 con MINITAB (1 de 7)
Ejercicio 4.1 con MINITAB (2 de 7)
Ejercicio 4.1 con MINITAB (3 de 7)
Ejercicio 4.1 con MINITAB (4 de 7)
Ejercicio 4.1 con MINITAB (5 de 7)
Ejercicio 4.1 con MINITAB (6 de 7)
Ejercicio 4.1 con MINITAB (7 de 7)
DISEÑO DE BLOQUES COMPLETOS ALEATORIZADOS

EJERCICIO
Utilizar un diseño aleatorio con bloques completos
para determinar si el compuesto químico afecta la
resistencia a la tensión de piezas de tela.

Piezas de Ropa
Compuesto A B C D E
1 73 68 74 71 67
2 73 67 75 72 70
3 75 68 78 73 68
4 73 71 75 75 69
 DISEÑOS CUADRADOS LATINOS

Estos tipos de diseño son una ampliación de los diseños

con bloques completos aleatorizados, en los que se trata
de aislar los efectos de dos fuentes de variabilidad de
los efectos atribuibles a los tratamientos. Estos diseños
pueden modelarse usando tablas de ANOVA en los
cuales hay tres fuentes de variación diferenciadas.
 DISEÑOS CUADRADOS LATINOS

EJEMPLO

Se está preparando un explosivo con cinco diferentes

formulaciones y se desea saber si el poder del mismo está
afectado por la constitución de la formulación. Se forma el
explosivo a partir de cinco lotes de materias primas, y la
mezcla es preparada por cinco diferentes operarios. Se diseñó
un experimento para medir el efecto de la formulación en el
poder explosivo de la mezcla, y se quiere separar las
variaciones atribuibles a los operadores y al lote de materia
prima. En la siguiente tabla vemos las diferentes corridas y los
resultados de las pruebas del poder explosivo de las mezclas,
según el diseño experimental seleccionado.
 DISEÑOS CUADRADOS LATINOS

Diseño cuadrado latino para la formulación explosiva

Operador
Lote 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=22 A=30 B=20 C=29 D=31

Las letras A,B,C,D, y E, representan las cinco formulaciones a ser

evaluadas. En estos diseños se debe probar cada formulación en
cada bloque solo una vez. Debe haber el mismo número de
ocurrencias en cada bloque, de ahí viene el nombre de diseños
cuadrados. El adjetivo “latino” proviene de las letras usadas para
los niveles del factor a estudiar.
 DISEÑOS CUADRADOS LATINOS

TABLA ANOVA PARA DISEÑOS CUADRADOS LATINOS

Fuente de Suma de Grados de Cuadrado

Variación Cuadrados Libertad Medio F0
Tratamientos SS tratamiento p-1 SS tratamiento MS tratamiento
p 1 MS e

Filas SS fila
p-1 SS fila

p 1

Columnas SS columna
p-1 SS columna

p 1
SS e
Error (p -2)(p -1) SS e

( p  1)( p  2)
SS t 2
Total p -1
 DISEÑOS CUADRADOS LATINOS

Donde:
2 2
p y y
SS tratamiento    p = Número de filas o columnas
. j. ...

j 1 p N

2 2 N = Número total de observaciones

p
y y
SS fila   i 1 p
i ..

N
...
i = índice para la fila (lotes)

2 2
j = índice para el tratamiento
p
y y
SS  ..k
 ...
columna
k 1 p N (formulación)
2
p p p
y k = índice para la columna (operador)
SS     y
2
t
 ...

i 1 j 1 k 1
ijk N
 DISEÑOS CUADRADOS LATINOS

Codificando los datos para simplificar las operaciones, restando 25

de cada observación, tenemos la siguiente tabla:

Diseño cuadrado latino para la formulación explosiva

Tabla codificada Los totales para
Operador los tratamientos
Lote 1 2 3 4 5 yi.. son:
1 A=-1 B=-5 C=-6 D=-1 E=-1 -14
2 B=-8 C=-1 D=5 E=2 A=11 9 A y.1. = 18
3 C=-7 D=13 E=1 A=2 B=-4 5 B y.2. = -24
4 D=1 E=6 A=1 B=-2 C=-3 3 C y.3. = -13
5 E=-3 A=5 B=-5 C=4 D=6 7
D y.4. = 24
y..k -18 18 -4 5 9 y...=10
E y.5. = 5
 DISEÑOS CUADRADOS LATINOS

Las ecuaciones de estos resultados son:

2

 680 
(10)  676.0
SS t
25

(14)  9  5  3  7  10
2 2 2 2 2 2

SS  lote
 68.0
5 25

(18)  18  (4)  5  9
2 2 2 2 2 2

SS operador
  10  150.0
5 25

(18)  (24)  (13)  24  5  10

2 2 2 2 2 2

SS tratamiento
  330.0
5 25

SS  SS  SS
e t lote
 SS operario  SS tratamiento  676.0  68.0  150.0  330.0  128.0
Ejercicio 4.2 con MINITAB (1 de 3)
Ejercicio 4.2 con MINITAB (2 de 3)
Ejercicio 4.2 con MINITAB (3 de 3)
DISEÑO CON BLOQUES INCOMPLETOS BALANCEADOS

En ciertos experimentos con frecuencia no es posible

probar todas las combinaciones de los tratamientos para
cada bloque. Situaciones como esta suelen ocurrir
debido a falta de material para pruebas o el tamaño
físico del bloque. En casos como estos recurrimos al uso
de diseños con bloques incompletos.

Se dice que el diseño es balanceado, cuando la

combinación de tratamientos en cada bloque es igual, es
decir, cada par de tratamientos ocurre junto el mismo
número de veces en cada bloque.
DISEÑO CON BLOQUES INCOMPLETOS BALANCEADOS

RESULTADOS DE EXPERIMENTO DE CATALÍTICO

LOTE DE MATERIA PRIMA

Catalítico 1 2 3 4 yi.
1 73 74 - 71 218
2 - 75 67 72 214
3 73 75 68 - 216
4 75 - 72 75 222
y.j 221 224 207 218 870= y...
DISEÑO CON BLOQUES INCOMPLETOS BALANCEADOS

TABLA ANOVA PARA DISEÑOS

DE BLOQUES INCOMPLETOS BALANCEADOS

Fuente de Suma de Grados de Cuadrados F0

Variación Cuadrados Libertad medios
kQ
2
Tratamientos
i
a-1 SS tratamiento ( ajustado MS tratamiento ( ajustado
(ajustados)  a a 1 MS e

2
Bloque y y
2
b-1 SS Bloque

 
.j ..
b 1
k N
Error SS e
N -a - b + 1 SS e

N  a b 1
2
Total y N-1
 y
2
 ..
ij N
DISEÑO CON BLOQUES INCOMPLETOS BALANCEADOS

Donde:
2
y
 y
2
SS T
  ..

i j
ij N

SS T
 SS tratamiento ( ajustado )
 SS bloque  SS E

2 2
b y y
  
.j ..
SS bloque
j 1 k N

a
kQ
2

SS  i 1
tratamiento ( ajustado )
 a

1 b
Q  y   nij y i=1,2,...,a
i i. k j 1 .j
 DISEÑO CON BLOQUES INCOMPLETOS BALANCEADOS

Sustituyendo en las fórmulas tenemos:

2

 63,156 
(870)  81.00
SS T
12

  224  218
2 2 2 2 2

SS bloque
 221 207  870  55.00
3 12

Q  218  13 (221  224  218)  9 / 3 Q  214  13 (207  224  218)  7 / 3

1 2

 216  13 (221  207  224)  4 / 3 Q  222  13 (221  207  218)  20 / 3

Q 3
4

SS tratamiento ( ajustado )


3 (9 / 3)  (7 / 3)  (4 / 3)  (20/ 3)
2 2 2 2
  22.75
(2)(4)

SSE = 81.00 - 22.75 - 55.00 = 3.25

DISEÑO CON BLOQUES INCOMPLETOS BALANCEADOS

Sustituyendo en la tabla ANOVA tenemos:

Fuente de Suma de Grados de Cuadrado F0

Variación Cuadrados Libertad medio
Tratamient 22.75 3 7.58 11.66
o (ajustado)
Bloque 55.00 3 -
Error 3.25 5 0.65
Total 81.00 11

Del análisis de varianza mostrado en la tabla anterior, vemos

que para un nivel de significación de 0.05, F 0.05,3.5 = 5.41, lo
que indica que hay diferencia notable entre los catalíticos.
Ejercicio 4.3 con MINITAB (1 de 2)
Ejercicio 4.3 con MINITAB (2 de 2)