5.5 Diseño en Cuadrados Grecolatinos
5.5 Diseño en Cuadrados Grecolatinos
5.5 Diseño en Cuadrados Grecolatinos
5 DISEO EN
CUADRADOS
GRECOLATINOS
Unidad V
INTRODUCCIN
En los arreglos por bloques, se pueden analizar 4 factores,
introduciendo un cuarto factor o bloque en un diseo
cuadrado latino, siguiendo las mismas reglas utilizadas
para introducir un tercer factor en un diseo cuadrado de
dos factores.
A este cuarto factor o bloque se le denomina componente
griego, ya que se utilizan letras griegas para identificar sus
niveles.
A la adicin de un diseo cuadrado latino y un cuarto
factor, se le llama Diseo Cuadrado Grecolatino.
INTRODUCCIN
Un cuadrado grecolatino es una extensin del cuadrado
latino en el que se incluye una tercera variable de
control o variable de bloque.
Todos los factores deben tener el mismo nmero de
niveles K y el nmero de observaciones necesarias sigue
siendo K2.
Este diseo es, por tanto, una fraccin del diseo completo
en bloques aleatorizados con un factor principal y 3
factores secundarios que requerira K4 observaciones.
DEFINICIN
ORIGEN
Nacen de una investigacin
llevada a cabo por el
matemtico y fsico
Leonhard Paul Euler,
quien fuera muy distinguido
durante el siglo XVIII por sus
estudios realizados.
Entre sus muchos estudios
se interes en los
cuadrados mgicos y realiz
un trabajo algebraico ms
que numrico con estos
para ello introdujo las letras
latinas en vez de nmeros,
a estos les llam
ORIGEN
Un cuadrado grecolatino es el resultado de la superposicin
de dos cuadrados latinos, el cuadrado resultante no tiene
ninguna de las letras o signo repetido en fila o columna,
ello es posible si ambos cuadrados latinos son ortogonales
entre s.
Cuadrado grecolatino
Cuadr
o
1
2
3
4
A
B
C
D
B
A
D
C
C
D
A
B
Alfa
Beta
D
Gamma
C
B Delta
4
para un
cuadrado grecolatino especificado. Es decir, los subndices l toman
valores que dependen de la celdilla .
es un efecto constante, comn a todas las unidades.
es el efecto producido por el i-simo nivel del factor fila. Dichos efectos estn
sujetos a la restriccin
es el efecto producido por el j-simo nivel del factor columna. Dichos efectos
estn sujetos a la restriccin
es el efecto producido por el h-simo nivel del factor letra latina. Dichos
efectos estn sujetos a la restriccin
es el efecto producido por el p-simo nivel del factor letra griega. Dichos
efectos estn sujetos a la restriccin
son variables aleatorias independientes con distribucin N (0, ).
Descomposicin de la variabilidad
Simblicamente se puede escribir:
SCT = SCF + SCC + SCL + SCG + SCR
SCT suma total de cuadrados.
SCF suma de cuadrados debida al efecto fila.
SCC suma de cuadrados debida al efecto columna.
SCL suma de cuadrados debida a las letras latinas.
SCG suma de cuadrados debida a las letras griegas.
SCR suma de cuadrados del error.
Anlisis de varianza
Anlisis de varianza
Por lo tanto, la suma de cuadrados debida al factor letra griega
puede calcularse usando los totales de la letra griega. El error
experimental se reduce en esta cantidad.
Las hiptesis nulas de igualdad entre los renglones, entre las
columnas, entre los tratamientos de la letra latina y entre los
tratamientos de la letra griega pueden probarse dividiendo la
media de cuadrados correspondiente entre la media de
cuadrados del error.
La regin de rechazo es el extremo superior de la distribucin
F p 1, (p 3) (p 1).
Suma de Cuadrados
Grados de Libertad
Fuente de Variacin
Tratamiento de Letra
Latina
p 1
Tratamiento de Letra
Griega
p 1
Renglones
p 1
Columnas
Columnas
p
p 1
1
Error
Error
(p
(p 3)(p
3)(p 1)
1)
p2 1
Anlisis de varianza
Total
Ejemplo 1
En la obtencin de un determinado producto qumico se est
interesado en comparar 4 procedimientos.
Se supone que en dicha obtencin tambin puede influir la
temperatura, presin y tipo de catalizador empleado,
decidindose realizar un experimento en cuadrado grecolatino.
Para ello, se consideran 4 niveles de cada uno de estos
factores. La tabla adjunta muestra el cuadrado greco-latino
que resulta elegido y las cantidades de producto obtenidas.
Ejemplo 1
En dicha tabla:
Las filas representan el factor principal, procedimientos.
Las columnas representan el factor temperatura.
Las letras latinas(A, B, C, D) representan el factor presin.
Las letras griegas(, , , ) representan el factor tipo de
catalizador.
Ejemplo 1
Temperaturas
Procedimient
os
T1
T2
T3
T4
P1
12
13
Dy
13
P2
By
10
15
11
P3
Ay
P4
11
10
Cy
N=16
n=4
EJEMPLO 1
Temperaturas
Procedimi
entos
P1
P1
P2
P2
P3
P3
P4
P4
T1
C
C
B
By
y
D
D
A
A
T2
T3
12
12
10
10
7
7
11
11
B
B
C
C
Ay
Ay
D
D
29
29
5
5
6
6
841
T4
13
13
15
15
5
5
10
10
A
A
D
D
B
B
Cy
Cy
37
37
1369
13
13
11
11
43
43
42
42
1849
1849
1764
1764
5
5
8
8
Dy
Dy
A
A
C
C
B
B
7
7
9
9
24
24
38
38
576
576
1444
1444
41
41
40
40
1681
1600
5491
Ejemplo 1
Temperaturas
Procedimi
entos
P1
P2
P3
P4
T1
T2
T3
T4
C 5
B 12
13
Dy 13
By
C 10
D 15
A 11
D 7
Ay
B 5
C 7
A 11
D 10
Cy
PARA LATINO
A=
11+5+13+11
40
1600
B=
6+12+5+9
32
1024
C=
5+10+8+7
30
900
D=
7+10+15+13
45
2025
147
5549
SUMAS
Ejemplo 1
Temperaturas
Procedimi
entos
P1
P2
P3
P4
T1
T2
T3
T4
PARA GRECOLATINO
11+12+15+7
45
2025
5+10+5+11
31
961
C 5
B 12
13
Dy 13
6+5+8+13
32
1024
By
C 10
D 15
A 11
7+10+13+9
39
1521
D 7
Ay
B 5
C 7
147
5531
A 11
D 10
Cy
SUMAS
Ejemplo 1
Datos al cuadrado
1
C 5
25
B 12
144 A 13
169 Dy 13
169
507
By 6
36
C 10
100 D 15
225 A 11
121
482
D 7
49
Ay 5
A 11
12
1
D 10
231
369
B 5
25
C 7
49
148
100 Cy 8
64
B 9
81
366
25
483
420
1503
Ejemplo 1
Ejemplo 1
Ejemplo 1
PARA LATINO
A=
11+5+13+11
40
1600
B=
6+12+5+9
32
1024
C=
5+10+8+7
30
900
D=
7+10+15+13
45
2025
147
5549
SUMAS
Ejemplo 1
PARA GRECOLATINO
11+12+15+7
45
2025
5+10+5+11
31
961
6+5+8+13
32
1024
7+10+13+9
39
1521
147
5531
SUMAS
Ejemplo 1
Coeficiente de
variacin
Filas
Filas
Grados de
libertad (Y)
n-124-1=3
4-1=3
Ejemplon-1
(procedimiento
(procedimiento
))
Columnas
Columnas
(temperatura)
(temperatura)
Tratamiento
Tratamiento
latino
latino
Tratamiento
greco
ERROR
TOTAL
n-1
n-1
4-1=3
4-1=3
n-1
n-1
4-1=3
4-1=3
n-1
4-1=3
(n-3)(n-1) (4-3)
(4-1)
(1)(3)=3
n2-1
42-1
16-1=15
Suma de
cuadrados
Cuadrados
medios
57.6875
57.6875
CMF=57.6875/3
CMF=57.6875/3
=19.2291
=19.2291
22.1875
22.1875
CMC=22.1875/3
CMC=22.1875/3
=7.3958
=7.3958
36.6875
36.6875
32.1875
3.6875
152.4375
CML=36.6875/3
CML=36.6875/3
=12.2291
=12.2291
CMG=32.1875/3
=10.7291
CME=3.6875/3
=1.2291
F Calculada
Ejemplo 1
F de tablas
v1 ; v2 con
Hiptesis a probar:
Ejemplo 1
Para procedimientos (filas)
Ejemplo 1
Para tratamiento latino (presin)
Ejemplo 1
CONCLUSIN:
Se concluye que se aceptan las hiptesis de igualdad de efectos de
columnas y de letra griega y se rechazan las hiptesis de igualdad de
efecto de filas y de letra latina. Es decir, son significativos los efectos
de los procedimientos y presin, pero no lo son los efectos de la
temperatura y catalizador.
Ejemplo 2
Un experimentador opina que las lneas de ensamblaje son
fuentes de variacin al momento de reproducir la frmula
para la elaboracin de dinamita. Para comprobarlo disea
un arreglo de cuadrado Grecolatino.
Supngase que en el experimento para comparar las
frmulas para la dinamita, tiene importancia el factor
adicional lnea de montaje. Considrese que existen cinco
lneas de montaje representadas por las letras griegas, ,
, y .
Ejemplo 2
El diseo de cuadrado grecolatino 5 x 5 que resulta se
muestra en la Tabla.
Operadores
Lotes
Materia
Prima
1
2
3
4
5
1
A
B
Cy
D
E
2
24
17
18
26
22
By
C
D
E
A
3
20
24
38
31
30
C
D
E
Ay
B
4
19
30
26
26
20
D
Ey
A
B
C
5
24
27
27
23
29
E
A
B
C
Dy
24
36
21
22
31
N=25
n=5
Ejemplo 2
Operadores
Lotes
Materia
Prima
24
By
20
19
24
24
111
12 321
B
Cy
Cy
D
D
17
18
18
26
26
C
D
D
E
E
24
38
38
31
31
D
E
E
Ay
Ay
30
26
26
26
26
Ey
A
A
B
B
27
27
27
23
23
A
B
B
C
C
36
21
21
22
22
134
130
130
128
128
17
16
16
16
16
E
E
22
22
A
A
30
30
B
B
20
20
C
C
29
29
Dy
Dy
31
31
132
132
107
107
143
143
121
121
130
130
134
134
17 424
17 424
80985
11449
20449
14641
1690
0
3
3
4
4
5
5
1795 81 395
6
956
900
900
384
384
Ejemplo 2
Operadores
Lotes
Materi
a
Prima
PARA LATINO
1
A=
24+30+26+27+3
143
B=
6
17+20+20+23+2
101
10 201
1
18+24+19+29+2
112
122 544
D=
2
26+38+30+24+3
149
22 201
1
22+31+26+27+2
130
16 900
635
82295
24
By
20
19
D 24 E
24
17
24 D 30
Ey 27 A
36 C=
Cy
18
38
26
A 27 B
21
26
31
Ay
26
B 23 C
22
22
30
20
C 29 Dy
31 E=
2
20 449
SUMAS
Ejemplo 2
Operadores
Lotes
Materi
a
Prima
PARA GRECOLATINO
1
24
By
20
19
D 24 E
24
17
24 D 30
Ey 27 A
36
Cy
18
38
26
A 27 B
21
26
31
Ay
26
B 23 C
22
22
30
20
C 29 Dy
31
24+31+30+29+2
135
2
18225
1
17+30+26+24+2
119
14161
2
18+20+26+27+3
122
14884
1
26+24+20+27+2
121
14641
4
22+38+19+23+3
138
19044
635
80955
SUMAS
Ejemplo 2
Datos al cuadrado
A
24 576 By
20 400 C
19 361 D 24 576 E
5
24
576
2,489
B 17 289 C
24 556 D 30 900 Ey
27 729 A
36 1296
3,790
Cy
18 324 D
144
38
E
4
26 676 A
27 729 B
21
441
3,614
D 26 676 E
31 961 Ay
26 676 B
23 529 C
22
484
3,326
22 484 A 30 900 B
20 400 C
29 841 Dy 31
961
3,856
2,349
4,281
3,013
3,404
3,758
16,805
Ejemplo 2
Ejemplo 2
Ejemplo 2
PARA LATINO
A=
24+30+26+27+3
143
2
20 449
B=
6
17+20+20+23+2
101
10 201
C=
1
18+24+19+29+2
112
122 544
D=
2
26+38+30+24+3
149
22 201
E=
1
22+31+26+27+2
130
16 900
Ejemplo 2
PARA GRECOLATINO
24+31+30+29+21
135
17+30+26+24+22
119
14161
18+20+26+27+31
122
14884
26+24+20+27+24
121
14641
22+38+19+23+36
138
19044
635
80955
SUMAS
2
18225
Ejemplo 2
Coeficiente de
variacin
Filas (lotes)
Columnas
(operador)
Tratamiento
latino
Tratamiento
greco
ERROR
TOTAL
Grados de
libertad (Y)
n-1
5-1=4
n-1
5-1=4
n-1
5-1=4
n-1
5-1=4
(n-3)(n-1) (5-3)
(5-1)
(2)(4)=8
n2-1
52-1
25-1=24
Suma de
cuadrados
Cuadrados
medios
F Calculada
68
CMF=68/4=17
150
CMC=150/4=37.5
330
CML=330/4=82.5
62
CMG=62/4=15.5
66
CME=66/8=8.25
676
160.75
Ejemplo 2
F de tablas
v1 ; v2 con
Hiptesis a probar:
Ejemplo 2
Para lote (filas)
Ejemplo 2
Para tratamiento latino
10 > 3.84
Ejemplo 2
CONCLUSIN:
Al parecer hay varianza en los Lotes de materia prima y a un
ms en los operadores lo que significa que hay una gran
fuente de variacin para los tratamientos.
EJEMPLO 3
En el caso de ejemplo 3.1, donde se comparan cuatro
mtodos de ensamble y se tiene el factor de bloque
operador, se podran tener dos factores de bloque
adicionales: orden en el que se hace el ensamble y lugar
donde se hace. De acuerdo con esto, el diseo en cuadro
grecolatino se observa en la tabla
EJEMPLO 3
EJEMPLO 3
Ejercicio 1
Se compara el rendimiento de tres procesos de fabricacin
(A, B y C) en tres condiciones experimentales (, y ) tres
das distintos con tres procedimientos de medicin. El
diseo y los resultados obtenidos se indican en el cuadro.
Da 1
Da 2
Da 3
Mtodo 1
10
13
11
Mtodo 2
10
Mtodo 3
11
10
Ejercicio 2
Una compaa distribuidora ubicada en los suburbios est
interesada en estudiar la diferencia en costos (tiempo y
gasolina) entre las cuatro rutas (A, B, C, D) que llegan a la
zona comercial, ms importante para ellos, en el otro
extremo de la ciudad. Deciden correr un experimento en
cuadro grecolatino controlando los factores de bloque
chofer, marca de vehculo (, , , ) y da de la semana. El
experimento se repite en dos semanas diferentes, en
las cuales no hay das festivos ni quincenas. Los costos
observados en pesos se muestran en la siguiente tabla:
NOTA: se puede realizar un promedio para las
dos semanas
Ejercicio 2
Ejercicio 3
El rendimiento de un proceso qumico se midi utilizando:
cinco lotes de materia prima,
cinco concentraciones del cido,
cinco tiempos de procesamiento (A, B, C, D y E) y
cinco concentraciones del catalizador (, , , , ).
Se us el cuadrado grecolatino siguiente. Analizar los datos
de este experimento (utilizar = 0,05) y sacar conclusiones.
Ejercicio 3
Referencias bibliogrficas
Universidad nacional del santa Facultad de ingeniera E.A.P.
Ingenieria agroindustrial(Noviembre 2013), Diseo de
cuadrado latino y grecolatino, pp. 21-31, Fecha de consulta
11 de noviembre de 2015 Disponible desde http://bit.ly/1PUvf9P
Reyes Moreno, K. A.(28 de julio de 2015) DISEOS EN
CUADRADOS GRECO-LATINOS(presentacin de PowerPoint)
Fecha de consulta 11 de noviembre de 2015 URL: http://
bit.ly/1I5h6Wc
Gutirrez Pulido, M. De la Vara Salazar, R.. (2008). Anlisis y
diseo de experimentos. Mxico: McGraw-Hill.
Referencias bibliogrficas
Universidad de Castilla~La mancha; Lpez Fidalgo, J(2009),
Modelos lineales 3 Diplomatura de Estadstica Ejercicios;
pp. 4-5; Fecha de consulta: 22 de noviembre de 2015; URL:
http://bit.ly/1PITI3D
Jimnez Gonzlez, R. (Agosto de 2012) Estadstica
Inferencial II p. 107. Fecha de consulta 10 de noviembre de
2015; URL http://bit.ly/1PUww0F