RCBD na danych NHANES

W naszym zablokowanym eksperymencie badamy efekt interwencji, w której lekarz zaleca pacjentowi ograniczenie tłuszczu lub kalorii w diecie – i sprawdzamy, czy ma to wpływ na wagę (bmxwt). Jako czynnik blokujący wybieramy płeć, która w zbiorze danych NHANES jest zapisana w zmiennej riagendr. Blokowanie służy tworzeniu grup eksperymentalnych jak najbardziej zbliżonych do siebie pod względem wybranych cech. Blokowanie tego eksperymentu według płci oznacza, że jeśli zaobserwujemy wpływ interwencji na bmxwt, z większym prawdopodobieństwem będzie on wynikał z samego zabiegu, a nie z płci uczestnika.

W kodzie R zablokowany eksperyment zapisujesz przy użyciu formuły w postaci: wynik ~ leczenie + czynnik_blokujący w odpowiedniej funkcji modelującej.

Dostępny jest obiekt nhanes_final.

To ćwiczenie jest częścią kursu

Projektowanie eksperymentów w R

Zobacz kurs

Instrukcje do ćwiczenia

Użyj aov(), aby utworzyć obiekt nhanes_rcbd. Pamiętaj, że zmienną z interwencją jest mcq365d, zmienną wynikową – bmxwt, a czynnikiem blokującym – riagendr.
Sprawdź wyniki nhanes_rcbd za pomocą funkcji summary().
Użyj funkcji z pakietu dplyr, aby zbadać średnie wartości wagi według zmiennych mcq365d i riagendr.

Interaktywne ćwiczenie praktyczne

Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.

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

Edytuj i uruchom kod