2/4/2021
Desarrollado para:
1.1. Controlar efecto de condiciones ajenas al experimento
1.2. Combinar varios factores
Arreglo en filas y columnas de factores experimentales
2.1. Cada factor aparece 1 vez por fila y columna
2.2. Factor con \(N\) niveles \(\rightarrow\) cuadro \(N \times N\)
Control de efectos ajenos
1.1. 1 Factor experimental con 3 niveles:
\[ \mathrm{Factor\ A} = \{a, b, c \} \] \[ \mathrm{Factor\ B} = \{\alpha, \beta, \gamma \}\]
\(a, \alpha\) | \(b, \beta\) | \(c, \gamma\) |
---|---|---|
\(c, \beta\) | \(a, \gamma\) | \(b, \alpha\) |
\(b, \gamma\) | \(c, \alpha\) | \(a, \beta\) |
Sólo para estudios simétricos (factores con mismo número de niveles)
1.1. Cuadrados superpuestos
Cada combinación aparece 1 vez
Estadísticamente:
3.1. Sólo se estiman efectos aditivos (no hay grados de libertad para interacciones)
3.2. Se necesitan más cuadrados para estimar interacciones
3.2.1. Deben ser aleatorizados
Complejidad crece a ritmo \(N^2\)
Pocos tratamientos = Pocos grados de libertad
Se recomienda:
3.1. \(4 \geq N \leq 8\)
Si interacciones biológicas son posibles \(\rightarrow\) EVITAR
Cuando \(N\) tratamientos es manejable es un diseño robusto
Pueden existir réplicas –cuadrado latino extendido
2.1. Réplicas suficientes \(\rightarrow\) estimación de interacciones
Puede ser aleatorizado
Cada bloque es una réplica
\(N\) tratamientos = \(N\) réplicas
Para análisis de 1 sólo factor en R:
aov(y ~ x, datos) lm(y ~ x, datos)
aov(y ~ x1 + x2, datos) lm(y ~ x1 + x2, datos)
Donde el número de tratamientos de \(x_1\) es el mismo que de \(x_2\)
Cada cuadrado aleatorizado = bloque
Todos los tratamientos están contenidos en cada bloque
Factor de agrupamiento es el bloque
3.1. Efecto aleatorio por bloque
3.2. Bloques necesarios para interacciones = Tratamientos + 1
El análisis:
aov(y = x1 * x2 + Error(bloque/(x1*x2)), datos)
Frecuentemente en diseños CL se considera que
1.1. Hay dos factores de agrupación –filas y columnas 1.2. Es posible modelarlos con efectos aleatorios
library(lme4) lmer(y ~ x1 + x2 + (1|filas) + (1|columnas), datos)
El modelo es aditivo porque no podemos ajustar la interacción \(x_1 : x_2\)
Si no hay un gradiente identificado que tomar en cuenta \(\rightarrow\) EVITAR CL
Identificar el gradiente puede incluso ayudar a incluir en análisis y aislar varianza que explica dicho gradiente (modelo ANCOVA en Unidad 5)