독성학 데이터
한 독성학자는 어떤 화학 물질에 노출될수록 투여량(dose)이 증가하면 사망률이 높아지는지 궁금해합니다. 그녀는 시험 생물을 서로 다른 투여량에 노출시키고, 각 시험 수조에서 사망한 개체 수를 기록했습니다. 각 투여량은 서로 다른 세 날에 반복 측정되었습니다.
먼저 GLM을 만들어 보세요. 단순한 모델부터 구축하면 디버깅이 쉬워집니다. 또한 이후에 GLMER 결과와 비교할 예정입니다. 이 모델에는 dose와 replicate를 고정효과로 포함합니다.
다음으로 GLMER을 만들어 보세요. 이번에는 dose는 고정효과, replicate는 랜덤 효과 절편으로 설정합니다.
마지막으로 각 모델에서 coef()로 계수 추정치를 확인하세요. 모델마다 절편 추정치가 표시되는 방식이 어떻게 다른지 살펴보세요.
이 연습은 강의의 일부입니다
R로 배우는 계층 및 혼합 효과 모형
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
# Fit glm_out and look at its coefficient estimates
glm_out <- glm(___ ~ ___ + ___,
family = "___", data = df)
coef(___)