NHANES RCBD
Weet je nog dat ons geblokkeerde experiment een behandeling omvatte waarbij de arts de patiënt vraagt om het vet of de calorieën in hun dieet te verminderen, en dat we testen wat voor effect dit heeft op het gewicht (bmxwt). We willen blokkeren op geslacht, wat in de NHANES-gegevensset is opgeslagen als riagendr. Onthoud dat blokkeren wordt gedaan om experimentele groepen te creëren die zo vergelijkbaar mogelijk zijn. Dit experiment blokkeren op geslacht betekent dat, als we een effect van de behandeling op bmxwt zien, het waarschijnlijker is dat dit effect daadwerkelijk door de behandeling komt en niet door het geslacht van de persoon.
In je R-code geef je een geblokkeerd experiment aan met een formule als: outcome ~ treatment + blocking_factor in de juiste modelleringsfunctie.
nhanes_final is beschikbaar.
Deze oefening maakt deel uit van de cursus
Experimenteel ontwerp in R
Oefeninstructies
- Gebruik
aov()omnhanes_rcbdte maken. Onthoud dat de behandeling is opgeslagen inmcq365den dat je de uitkomstbmxwttest, met de blokkeringsfactorriagendr. - Bekijk de resultaten van
nhanes_rcbdmetsummary(). - Gebruik functies uit
dplyrom de gemiddelde gewichten permcq365denriagendrte bekijken.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# 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 = ___(___))