grecolatinos

Capítulo 8
Diseños en cuadrados greco-latinos
8.1. Introducción
El modelo en cuadrado greco-latino se puede considerar como una extensión del
cuadrado latino en el que se incluye una tercera variable de control o variable de bloque.
En este modelo, como en el diseño en cuadrado latino, todos los factores deben tener el
mismo número de niveles K y el número de observaciones necesarias siguesiendo K2. Este
diseño es, por tanto, una fracción del diseño completo en bloques aleatorizados con un
factor principal y 3 factores secundarios que requeriría K4 observaciones.
Los cuadrados greco-latinos se obtienen por superposición de dos cuadrados latinos del
mismo orden y ortogonales entre sí, uno de los cuadrados con letras latinas el otro con
letras griegas. Dos cuadrados reciben elnombre de ortogonales si, al superponerlos, cada
letra latina y griega aparecen juntas una sola vez en el cuadrado resultante.
En el Apéndice C se muestra una tabla de cuadrados latinos que dan lugar, por superposición
de dos de ellos, a cuadrados greco-latinos. Notamos que no es posible formar
cuadrados greco-latinos de orden 6.
La Tabla 5-8 ilustra un cuadrado greco-latino para K = 4
Tabla5-8.
Cuadrado greco-latino
A α B β C γ D δ
D γ C δ B α A β
B δ A γ D β C α
C β D α A δ B γ
1
2 Diseños en cuadrados greco-latinos
8.2. Planteamiento del modelo
En un diseño en cuadrado greco-latino la variable respuesta yij(hp) viene descrita
por la siguiente ecuación
yij(hp) = μ + τi + βj + γh + δp + ǫij(hp)


i = 1, 2 . . . ,K
j = 1, 2 . . . ,K
h = 1, 2 . . . ,K
p = 1, 2. . . ,K
, (8.1)
donde
μ es un efecto constante, común a todas las unidades.
τi es el efecto producido por el i-ésimo nivel del factor fila. Dichos efectos están
sujetos a la restricción

i τi = 0.
βj es el efecto producido por el j-ésimo nivel del factor columna. Dichos efectos están
sujetos a la restricción

j βj = 0.
γh es el efecto producido por el h-ésimo nivel del factor letralatina. Dichos efectos
están sujetos a la restricción

h γh = 0.
δp es el efecto producido por el p-ésimo nivel del factor letra griega. Dichos efectos
están sujetos a la restricción

p δp = 0.
ǫij(hp) son variables aleatorias independientes con distribución N(0, σ).
La notación yij(hp) indica que los niveles i y j determinan los niveles h y p para
un cuadrado greco-latinoespecificado. Es decir, los subíndices h y p toman valores que
dependen de la celdilla (i, j).
Se utiliza la siguiente notación:
N = K2 es el número total de observaciones.
El total y el promedio de todas las observaciones
y.... =

i

j
yij(..) ¯y.... =
y....
K2
El total y el promedio para cada fila
yi... =
K
j=1
yij(..) ¯yi... =
yi...
K
(8.2)
8.3 Estimación de los parámetros delmodelo 3
El total y el promedio para cada columna
y.j.. =
K
i=1
yij(..) ¯y.j.. =
y.j..
K
(8.3)
El total y el promedio para cada letra latina
y..h. =

i,j
yij(h.) ¯y..h. =
y..h.
K (8.4)
El total y el promedio para cada letra griega
y...p =

i,j
yij(.p) ¯y...p =
y...p
K (8.5)
y..h. se obtiene sumando las K observaciones en las que la letra latina se ha fijado al
nivel h.
y...pse obtiene sumando las K observaciones en las que la letra griega se ha fijado al
nivel p.
Comentario 8.1
Uno de los inconvenientes del cuadrado greco-latino, al igual que el cuadrado latino, es
que requiere el mismo número de niveles para los cuatro factores que intervienen. Además
no hay cuadrados greco-latinos de dimensión 6.
8.3. Estimación de los parámetros del modelo
Siguiendo elmismo proceso que en los diseños anteriores se obtienen los siguientes
estimadores máximos verosímiles de los parámetros del modelo
μ =

i

j
yij(..)
K2 = ¯y.... , (8.6)
τi =
1
K
K
j=1
yij(..) − μ = ¯yi... − ¯y.... , (8.7)
βj =
1
K
K
i=1
yij(..) − μ = ¯y.j.. − ¯y.... , (8.8)
4 Diseños en cuadrados greco-latinos
γh =
1
K

i,j
yij(h.) − μ = ¯y..h. − ¯y.... ,...