시작하기무료로 시작하기

귀무가설 검정

귀무가설 검정은 p-값을 사용해 어떤 변수가 0과 유의미하게 다른지 판단해요. 최근에는 귀무가설 검정과 p-값의 남용·과용 문제로 인해 American Statistical Association이 p-값 사용에 관한 성명을 발표하기도 했습니다.

이와 같은 비판과 기타 수치적 어려움 때문에, lme4 패키지의 제작자인 Doug Bates는 자신의 패키지에 p-값을 포함하지 않았어요. 그럼에도 불구하고 p-값을 추정해야 할 필요가 있을 수 있습니다. 이를 위해 lmerTest 패키지를 포함한 여러 패키지가 제공됩니다.

lmerTest는 lme4 패키지와 동일한 lmer() 문법을 사용하지만, 출력은 일부 다릅니다. 이번 연습에서는 lmerTest와 lme4를 사용해 각각 lmer() 모델을 적합해 보겠습니다.

이 연습은 강의의 일부입니다

R로 배우는 계층 및 혼합 효과 모형

강의 보기

실습형 인터랙티브 연습

이 예제를 이 샘플 코드를 완성하여 풀어보세요.

# Load lmerTest
library(___)

# Fit a lmer use lme4
out_lme4 <- 
	lme4::___(Crime ~ Year2 + (1 + Year2 | County), 
              data = md_crime)

# Fit a lmer use lmerTest
out_lmerTest <- 
	lmerTest::___(Crime ~ Year2 + (1 + Year2 | County), 
                  data = md_crime)

# Look at the summaries
summary(out_lme4)
summary(out_lmerTest)
코드 편집 및 실행