BIBD - função renal de gatos
Para garantir que você realmente entendeu como é um BIBD, vamos montar do zero um conjunto de dados que contenha um BIBD.
Suponha que a gente queira testar a diferença entre quatro tipos de ração úmida na dieta de gatos e seu efeito na função renal. No entanto, ração para gatos é cara, então vamos testar apenas 3 alimentos por bloco para economizar. O fator de bloqueio é a color do gato, já que não estamos interessados nisso como parte do experimento. O desfecho será o nível de creatinina no sangue, um indicador de função e disfunção renal tanto em gatos quanto em humanos.
Este exercício faz parte do curso
Planejamento de Experimentos em R
Instruções do exercício
- A função personalizada
lambda()já foi carregada para você. Calculelambdacomt= 4,k= 3 er= 3 para verificar se um BIBD é possível. - Execute o código para montar o conjunto de dados. Você pode ver a ordem em que os tratamentos de alimento são usados em cada bloco.
- Crie
cat_modelcomaov()de acordo com a descrição do experimento acima e examine os resultados comsummary(). O tipo de ração úmida faz diferença nos níveis de creatinina?
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# 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 = ___)
___