NHANES RCBD
Rappelez-vous que notre expérience avec blocs comportait un traitement où le médecin demande au patient de réduire les graisses ou les calories dans son alimentation, et nous testons l'effet que cela a sur le poids (bmxwt). Nous prévoyons de bloquer par sexe, qui dans le jeu de données NHANES est stocké dans riagendr. Souvenez-vous que le blocage vise à créer des groupes expérimentaux aussi similaires que possible. Bloquer cette expérience par sexe signifie que si nous observons un effet du traitement sur bmxwt, il est plus probable que l'effet soit réellement dû au traitement plutôt qu'au sexe de l'individu.
Dans votre code R, vous indiquez une expérience avec blocs en utilisant une formule de la forme : outcome ~ treatment + blocking_factor dans la fonction de modélisation appropriée.
nhanes_final est disponible.
Cet exercice fait partie du cours
Plan d’expériences en R
Instructions
- Utilisez
aov()pour créernhanes_rcbd. Rappelez-vous que le traitement est stocké dansmcq365det que vous testez la variable de résultatbmxwt, avec le facteur de blocageriagendr. - Examinez les résultats de
nhanes_rcbdavecsummary(). - Utilisez les fonctions de
dplyrpour examiner les poids moyens selonmcq365detriagendr.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# Use aov() to create nhanes_rcbd
___
# Check results of nhanes_rcbd with summary()
___
# Print mean weights by mcq365d and riagendr
___ %>%
group_by(___, ___) %>%
summarize(mean_wt = ___(___))