DISEO DE EXPERIMENTOS

Ing. Felipe Llaugel

Diseos con bloques aleatorizados
Ing. Felipe Llaugel
En muchos problemas de experimentos, es necesario
hacer un diseo de tal manera que la variabilidad
proveniente de fuentes conocidas pueda ser
sistemticamente controlada. El diseo 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.
DISEO DE BLOQUES COMPLETOS ALEATORIZADOS
Para los fines del anlisis de varianza el bloqueo introduce
un efecto adicional ficticio, cuyo objetivo es separar del
error experimental, alguna fuente de variabilidad
conocida.
Ing. Felipe Llaugel
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 presin
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 podra afectar la
medida registrada, pero se sabe que las muestras no
necesariamente provienen de un material homogneo. 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:
Ejemplo:
DISEO DE BLOQUES COMPLETOS ALEATORIZADOS
Ing. Felipe Llaugel
DISEO 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
Ing. Felipe Llaugel
DISEO DE BLOQUES COMPLETOS ALEATORIZADOS
TABLA ANOVA
Fuente de
Variacin
Suma de
Cuadrados
Grados de
Libertad
Cuadrados
medios
F
0
Tratamiento
i
i
a
y y
b N
. ..
2
1
2

a - 1
Tratamiento SS
a 1
Tratamiento
E
MS
MS
Bloque
.
..
j
j
b
y
y
a N
2
1
2

b - 1 Bloque SS
b 1
Error E SS (Por
diferencia)
(a - 1)(b -
1)
E SS
a b ( )( ) 1 1
Total
ij
j
b
i
a
y
y
N
2
1 1
2



..
N - 1
Ing. Felipe Llaugel
DISEO 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 i
y
.
1 -2 -1 1 5 3
2 -1 -2 3 4 4
3 -3 -1 0 2 -2
4 2 1 5 7 15
. j
y -4 -3 9 18
..
y
20
Ing. Felipe Llaugel
DISEO DE BLOQUES COMPLETOS ALEATORIZADOS
Clculo de resultados
Tratamiento SS

2 2 2 2 2
3 4 2 15 20
4 16
3850 .
Bloque SS




2 2 2 2 2
4 3 9 18 20
4 16
8250 .
T SS
154
16
129 00
2
20
.
E T tratamiento bloque SS SS SS SS
129 00 3850 82 50 8 . . .
Ing. Felipe Llaugel
DISEO DE BLOQUES COMPLETOS ALEATORIZADOS
TABLA ANOVA
Fuente de
Variacin
Suma de
Cuadrados
Grados de
Libertad
Cuadrados
Medios F
0
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 significacin 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)
Ing. Felipe Llaugel
DISEO DE BLOQUES COMPLETOS ALEATORIZADOS
Utilizar un diseo aleatorio con bloques completos
para determinar si el compuesto qumico afecta la
resistencia a la tensin de piezas de tela.
EJERCICIO
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
Ing. Felipe Llaugel
DISEOS CUADRADOS LATINOS
Estos tipos de diseo son una ampliacin de los diseos
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 diseos
pueden modelarse usando tablas de ANOVA en los
cuales hay tres fuentes de variacin diferenciadas.
Ing. Felipe Llaugel
Se esta preparando un explosivo con cinco diferentes
formulaciones y se desea saber si el poder del mismo esta
afectado por la constitucin de la formulacion. Se forma el
explosivo a partir de cinco lotes de materias primas, y la
mezcla es preparada por cinco diferentes operarios. Se diseo
un experimento para medir el efecto de la formulacion 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,
segn el diseo experimental seleccionado.
DISEOS CUADRADOS LATINOS
EJEMPLO
Ing. Felipe Llaugel
DISEOS CUADRADOS LATINOS
Diseo cuadrado latino para la formulacion 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 diseos se debe probar cada formulacion en
cada bloque solo una vez. Debe haber el mismo numero de
ocurrencias en cada bloque, de ah viene el nombre de diseos
cuadrados. El adjetivo latino proviene de las letras usadas para
los niveles del factor a estudiar.
Ing. Felipe Llaugel
DISEOS CUADRADOS LATINOS
TABLA ANOVA PARA DISEOS CUADRADOS LATINOS
Fuente de
Variacin
Suma de
Cuadrados
Grados de
Libertad
Cuadrado
Medio F
0
Tratamientos tratamiento SS p - 1
tratamiento SS
p 1
tratamiento
e
MS
MS
Filas
fila SS
p - 1 fila SS
p 1
Columnas
columna SS
p - 1 columna SS
p 1
Error
e SS
(p -2)(p -1) e SS
p p ( )( ) 1 2
Total
t SS
p
2
- 1
Ing. Felipe Llaugel
DISEOS CUADRADOS LATINOS
Donde:
tratamiento
j
j
p
SS
y
y
p N

. .
...
2
1
2
fila
i
i
p
SS
y y
p N

.. ...
2
1
2
columna
k
k
p
SS
y y
p N

.. ...
2
1
2
t
ijk
k
p
j
p
i
p
SS
y
y
N



2
1 1
2
1
...
p = Numero de filas o columnas
N = Numero total de observaciones
I = ndice para la fila (lotes)
j = ndice para el tratamiento (formulacion)
k = ndice para la columna (operador)
Ing. Felipe Llaugel
DISEOS CUADRADOS LATINOS
Codificando los datos para simplificar las operaciones, restando 25
de cada observacin, tenemos la siguiente tabla:
Diseo cuadrado latino para la formulacion explosiva
Tabla codificada
Operador
Lote 1 2 3 4 5 y
i..
1 A=-1 B=-5 C=-6 D=-1 E=-1 -14
2 B=-8 C=-1 D=5 E=2 A=11 9
3 C=-7 D=13 E=1 A=2 B=-4 5
4 D=1 E=6 A=1 B=-2 C=-3 3
5 E=-3 A=5 B=-5 C=4 D=6 7
y
..k
-18 18 -4 5 9 y
...
=10
Los totales para los
tratamientos son:

A y
.1.
= 18
B y
.2.
= -24
C y
.3.
= -13
D y
.4.
= 24
E y
.5.
= 5
Ing. Felipe Llaugel
DISEOS CUADRADOS LATINOS
Las ecuaciones de estos resultados son:
t SS
680
25
6760
2
10 ( )
.
lote SS

2
2 2 2 2
2
14
9 5 3 7
10
5 25
680
( )
.
operador SS




2
2
2
2 2
2
18
18
4
5 9
10
5 25
1500
( ) ( )
.
tratamiento SS




2 2 2
2 2
2
18 24 13
24 5
10
5 25
3300
( ) ( ) ( )
.
e t lote operario tratamiento SS SS SS SS SS
676 0 68 0 150 0 330 0 128 0 . . . . .
Ing. Felipe Llaugel
DISEOS CUADRADOS LATINOS
TABLA ANOVA PARA DISEOS CUADRADOS LATINOS
Fuente de
Variacin
Suma de
Cuadrados
Grados de
Libertad
Cuadrado
Medio
F
0
Formulacion 330.0 4 82.5 7.73
Lotes 68.0 4 17.0
Operarios 150.0 4 37.5
Error 128.0 12 10.67
Total 676.0 24
Buscando en la tabla del estadstico F para = 0.05 y 4 y 12 grados
de libertad, tenemos que F
0.05,4,12
= 3.26, lo que indica que las
formulaciones difieren ya que se rechaza la hiptesis nula de que
sean iguales porque F
0
> F
0.05,4,12
.
Ing. Felipe Llaugel
Se dice que el diseo es balanceado, cuando la
combinacin de tratamientos en cada bloque es igual, es
decir, cada par de tratamientos ocurre junto el mismo
numero de veces en cada bloque.
DISEO 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 tamao
fsico del bloque. En casos como estos recurrimos al uso
de diseos con bloques incompletos.
Ejercicio 4.2 con MINITAB (1 de 3)
Ejercicio 4.2 con MINITAB (2 de 3)
Ejercicio 4.2 con MINITAB (3 de 3)
Ing. Felipe Llaugel
DISEO CON BLOQUES INCOMPLETOS BALANCEADOS
RESULTADOS DE EXPERIMENTO DE CATALTICO
LOTE DE MATERIA PRIMA
Cataltico 1 2 3 4 y
i.
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
...
Ing. Felipe Llaugel
DISEO CON BLOQUES INCOMPLETOS BALANCEADOS
TABLA ANOVA PARA DISEOS
DE BLOQUES INCOMPLETOS BALANCEADOS
Fuente de
Variacin
Suma de
Cuadrados
Grados de
Libertad
Cuadrados
medios
F
0
Tratamientos
(ajustados)
k
a
i
Q
2

a - 1
tratamiento ajustado SS
a
(
1
tratamiento ajustado
e
MS
MS
(
Bloque
.
..
j
y
y
k N
2
2

b - 1
Bloque SS
b 1
Error e SS N -a - b + 1
e SS
N a b 1
Total
ij
y
y
N
2
2

..
N - 1
Ing. Felipe Llaugel
DISEO CON BLOQUES INCOMPLETOS BALANCEADOS
Donde:
T
ij
j i
SS
y
y
N


2
2
..
T tratamiento ajustado bloque E SS SS SS SS

( )
bloque
j
j
b
SS
y
y
k N

.
..
2
1
2
tratamiento ajustado
i
i
a
SS
Q
k
a
( )

2
1

i i
ij
j
j
b
Q y
n
y
k

. .
1
1
i=1,2,...,a
Ing. Felipe Llaugel
DISEO CON BLOQUES INCOMPLETOS BALANCEADOS
Sustituyendo en las frmulas tenemos:
T SS
63156
12
8100
2
870
, .
( )
bloque SS



2 2 2 2 2
221 207 224 218 870
3 12
5500 .
1
218
1
3
221 224 218 9 3
Q
( ) /
2
214
1
3
207 224 218 7 3
Q
( ) /
3
216
1
3
221 207 224 4 3
Q
( ) /
4
222
1
3
221 207 218 20 3
Q
( ) /

tratamiento ajustado SS ( )
( / ) ( / ) ( / ) ( / )
( )( )
.

3
2 4
22 75
2 2 2 2
9 3 7 3 4 3 20 3
SS
E
= 81.00 - 22.75 - 55.00 = 3.25
Ing. Felipe Llaugel
DISEO CON BLOQUES INCOMPLETOS BALANCEADOS
Sustituyendo en la tabla ANOVA tenemos:
Fuente de
Variacin
Suma de
Cuadrados
Grados de
Libertad
Cuadrado
medio
F
0
Tratamient
o (ajustado)
22.75 3 7.58 11.66
Bloque 55.00 3 -
Error 3.25 5 0.65
Total 81.00 11
Del anlisis de varianza mostrado en la tabla anterior, vemos
que para un nivel de significacin de 0.05, F
0.05,3.5
= 5.41, lo
que indica que hay diferencia notable entre los catalticos.
Ejercicio 4.3 con MINITAB (1 de 2)
Ejercicio 4.3 con MINITAB (2 de 2)