BIBD - función renal de los gatos
Para asegurarnos de que entendemos cómo es un BIBD, vamos a construir desde cero un conjunto de datos que contenga un BIBD.
Imagina que queremos probar la diferencia entre cuatro tipos de comida húmeda en la dieta de los gatos sobre su función renal. Sin embargo, la comida para gatos es cara, así que solo probaremos 3 comidas por bloque para ahorrar dinero. El factor de bloqueo es el color del gato, ya que no nos interesa como parte del experimento. El resultado será el nivel de creatinina en sangre medido, un indicador de función y disfunción renal tanto en gatos como en humanos.
Este ejercicio forma parte del curso
Diseño experimental en R
Instrucciones del ejercicio
- La función personalizada
lambda()ya está cargada. Calculalambdacont= 4,k= 3 yr= 3 para comprobar que un BIBD es posible. - Ejecuta el código para construir el conjunto de datos. Puedes ver el orden en el que se usan los tratamientos de comida en cada bloque.
- Crea
cat_modelconaov()según la descripción del experimento anterior y examina los resultados consummary(). ¿El tipo de comida húmeda marca una diferencia en los niveles de creatinina?
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Calculate lambda
___
# Build the data.frame
creatinine <- c(1.98, 1.97, 2.35, 2.09, 1.87, 1.95, 2.08, 2.01, 1.84, 2.06, 1.97, 2.22)
food <- as.factor(c("A", "C", "D", "A", "B", "C", "B", "C", "D", "A", "B", "D"))
color <- as.factor(rep(c("Black", "White", "Orange", "Spotted"), each = 3))
cat_experiment <- as.data.frame(cbind(creatinine, food, color))
# Create cat_model and examine with summary()
___ <- ___ (___ ~ ___ + color, data = ___)
___