Data toksikologi
Seorang ahli toksikologi ingin mengetahui apakah paparan terhadap suatu bahan kimia meningkatkan mortalitas seiring kenaikan dosisnya. Ia memaparkan organisme uji pada berbagai dosis bahan kimia dan mencatat jumlah yang mati di setiap tangki uji. Setiap dosis diulang pada tiga hari yang berbeda.
Pertama, bangun sebuah GLM. Membangun model sederhana memudahkan proses debug. Selain itu, Anda nanti akan membandingkan hasilnya dengan GLMER. Model ini memasukkan dose dan replicate sebagai efek tetap (fixed-effects).
Kedua, bangun sebuah GLMER. Kali ini, dose adalah efek tetap dan replicate adalah intersep efek acak (random-effects intercept).
Terakhir, periksa taksiran koefisien dengan coef() dari tiap model. Perhatikan bagaimana taksiran intersep ditampilkan berbeda pada setiap model.
Latihan ini adalah bagian dari kursus
Model Hierarki dan Mixed Effects di R
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# Fit glm_out and look at its coefficient estimates
glm_out <- glm(___ ~ ___ + ___,
family = "___", data = df)
coef(___)