原假设检验
原假设检验使用p 值来检验某个变量是否与 0 存在"显著"差异。近期,由于对原假设检验和 p 值的滥用与过度使用,美国统计协会发布了关于p 值使用的声明。
由于上述批评及其他数值计算方面的挑战,Doug Bates(lme4 包的作者)没有在其包中提供 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)