NHANES RCBD

Připomeňme si, že naš blokový experiment zahrnoval léčbu, při níž lékař požádá pacienta, aby omezil tuk nebo kalorie ve stravě – a my testujeme, jaký vliv to má na hmotnost (bmxwt). Plánujeme blokovat podle pohlaví, které je v datasetu NHANES uloženo jako riagendr. Blokování slouží k vytvoření experimentálních skupin, které jsou si co nejvíce podobné. Blokováním tohoto experimentu podle pohlaví dosáhneme toho, že pokud pozorujeme vliv léčby na bmxwt, je pravděpodobnější, že tento efekt skutečně pochází z léčby, a ne z pohlaví konkrétního jedince.

V R kódu se blokový experiment vyjadřuje vzorcem ve tvaru: outcome ~ treatment + blocking_factor v příslušné modelovací funkci.

K dispozici máš nhanes_final.

Toto cvičení je součástí kurzu

Experimental Design in R

Zobrazit kurz

Pokyny k cvičení

Pomocí aov() vytvoř objekt nhanes_rcbd. Léčba je uložena v mcq365d, testovaný výsledek je bmxwt a blokující faktor je riagendr.
Prozkoumej výsledky nhanes_rcbd pomocí funkce summary().
Pomocí funkcí z balíčku dplyr prozkoumej průměrné hmotnosti podle mcq365d a riagendr.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

# 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 = ___(___))

Upravit a spustit kód