提取系数
我们常常想知道模型估计得到的系数取值。虽然 summary() 函数会提供模型输出,但我们有时也希望直接访问这些结果。
固定效应的估计可直接使用 fixef() 函数提取。随机效应的估计可直接使用 ranef() 函数提取。我们还可以用 confint() 函数为固定效应提取置信区间。
broom.mixed 包还为 lmer() 模型提供了用于提取结果的 tidy 方法,即 tidy() 函数。不过,由于混合效应模型本身较为复杂,这些结果通常也更复杂,并不如许多 tidy 输出那样"整洁"。
已经为您拟合了一个 lmer 模型,并保存为 out。
本练习是课程的一部分
R 中的分层与混合效应模型
练习说明
- 使用保存的模型
out,通过fixef()提取固定效应系数。 - 使用保存的模型
out,通过ranef()提取随机效应系数。 - 使用保存的模型
out,通过confint()估计 95% 置信区间。 - 使用
tidy()配合out和conf.int = TRUE,用一条"整洁"的命令复现前面三步的代码调用。
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Extract the fixed-effect coefficients
___
# Extract the random-effect coefficients
___
# Estimate the confidence intervals
___
# Use the tidy function
___