Dane toksykologiczne
Toksykolog zastanawia się, czy ekspozycja na substancję chemiczną zwiększa śmiertelność wraz ze wzrostem dawki. Organizmy testowe były narażane na różne dawki substancji, a liczba padłych osobników była rejestrowana w każdym zbiorniku testowym. Każdą dawkę powtarzano przez trzy różne dni.
Na początku zbuduj model GLM. Prosty model ułatwia debugowanie. Będziesz też później porównywać jego wyniki z modelem GLMER. W tym modelu dose i replicate pełnią rolę efektów stałych.
Następnie zbuduj model GLMER. Tym razem dose jest efektem stałym, a replicate – losowym wyrazem wolnym.
Na koniec zbadaj oszacowania współczynników za pomocą coef() dla obu modeli. Zwróć uwagę, jak różnie wyświetlany jest wyraz wolny w każdym z nich.
To ćwiczenie jest częścią kursu
Hierarchiczne i mieszane modele efektów w R
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# Fit glm_out and look at its coefficient estimates
glm_out <- glm(___ ~ ___ + ___,
family = "___", data = df)
coef(___)