BIBD - fonction rénale du chat
Pour bien comprendre à quoi ressemble un BIBD, construisons un jeu de données de A à Z.
Supposons que nous voulions tester la différence entre quatre aliments humides différents dans l’alimentation des chats sur leur fonction rénale. Or, la nourriture pour chats est coûteuse, donc nous ne testerons que 3 aliments par bloc pour limiter les dépenses. Le facteur de blocage est la color du chat, car cet aspect ne nous intéresse pas dans l’expérience. Le critère de jugement sera le taux de créatinine sanguine mesuré, un indicateur de la fonction rénale et de la dysfonction, chez les chats comme chez les humains.
Cet exercice fait partie du cours
Plan d’expériences en R
Instructions
- La fonction personnalisée
lambda()a été chargée pour vous. Calculezlambdaavect= 4,k= 3 etr= 3 pour vérifier qu’un BIBD est possible. - Exécutez le code pour assembler le jeu de données. Vous pouvez voir l’ordre dans lequel les traitements alimentaires sont utilisés dans chaque bloc.
- Créez
cat_modelavecaov()conformément à la description de l’expérience ci-dessus et examinez les résultats avecsummary(). Le type d’aliment humide a-t-il un effet sur les niveaux de créatinine ?
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# 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 = ___)
___