Kiểm định giả thuyết không
Kiểm định giả thuyết không sử dụng p-value để kiểm tra xem một biến có khác đáng kể so với 0 hay không. Gần đây, việc lạm dụng và sử dụng quá mức kiểm định giả thuyết không và p-value đã khiến Hiệp hội Thống kê Hoa Kỳ đưa ra tuyên bố về việc sử dụng p-value.
Do những phê bình như vậy cùng các thách thức tính toán khác, Doug Bates (tác giả gói lme4) không đưa p-value vào gói của mình. Tuy vậy, bạn vẫn có thể muốn hoặc cần ước lượng p-value. Để đáp ứng nhu cầu này, có một số gói, bao gồm lmerTest.
lmerTest dùng cùng cú pháp lmer() như gói lme4, nhưng cho ra các đầu ra khác. Trong bài tập này, bạn sẽ khớp một mô hình lmer() bằng cả lmerTest và lme4.
Bài tập này là một phần của khóa học
Mô hình phân cấp và hiệu ứng hỗn hợp trong R
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# 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)