BIBD - fonction rénale du chat
Pour bien comprendre ce à quoi ressemble un BIBD, construisons à partir de zéro un jeu de données qui en contient un.
Supposons que nous voulions comparer quatre nourritures humides différentes dans l'alimentation des chats pour voir leur effet sur la fonction rénale. La nourriture pour chats étant coûteuse, nous ne testerons que 3 nourritures par bloc pour économiser. Le facteur de blocage est la color du chat, car ce n'est pas l'objet de notre expérience. Le résultat sera le taux de créatinine sanguine mesuré, un indicateur de la fonction rénale et des dysfonctions, autant chez les chats que chez les humains.
Cette activité fait partie du cours
Plan d'expériences avec R
Instructions de l’exercice
- 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 pourrez voir l'ordre dans lequel les traitements alimentaires sont utilisés dans chaque bloc.
- Créez
cat_modelavecaov()selon la description de l'expérience ci-dessus et examinez les résultats avecsummary(). Le type de nourriture humide a-t-il un effet sur les niveaux de créatinine ?
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# 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 = ___)
___