RCBD avec NHANES
Rappelez-vous que notre expérience avec blocs comprend un traitement où le médecin demande au patient de réduire les gras ou les calories dans son alimentation, et nous testons l'effet de cela sur le poids (bmxwt). Nous prévoyons de bloquer par sexe, variable qui se trouve dans l'ensemble de données NHANES sous riagendr. N'oubliez pas que le blocage vise à créer des groupes expérimentaux aussi semblables 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 la personne.
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.
Cette activité fait partie du cours
Plan d'expériences avec R
Instructions de l’exercice
- Utilisez
aov()pour créernhanes_rcbd. Rappelez-vous que le traitement est stocké dansmcq365det que vous testez l'issuebmxwt, 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 ce code d’exemple.
# 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 = ___(___))