Question

我正在尝试从线性混合效果模型中提取每个组的斜率并进行截距。该模型是使用lme4库中的lmer构建的，我可以使用jtools库中的interact_plot查看每个组的结果。如何获得每条线的斜率和截距？

我知道我可以使用summary（）或summ（）来查看固定效果的估算值和随机效果的方差，但是看不到随机效果的估算值。因此，我无法准确计算模型的斜率和截距。

>library(lme4)
> cond_waterxsilver <- lmer(LnAg ~ LnVolume + (LnVolume | FilterID) + SilverType + WaterType + SilverType*WaterType + SilverType*LnVolume +  WaterType*LnVolume, data=capwater_removed.data)


> library(jtools)
> interact_plot(cond_ranin_waterxsilver, pred = LnVolume, modx = WaterType, mod2 = SilverType)

我只是想从模型中获取斜率和六条线的截距（两种不同的WaterType和三种不同的SilverType）。 jtools或其他软件包中是否有工具可以帮助我从模型中提取斜率和截距？

Answer 1

我是这个软件包的开发者！

简短说明：我将要提到的这个功能和另一个功能刚刚移到了一个名为interactions的新程序包中，该程序包正在添加到CRAN中。假设您尚未更新到jtools的最新版本（2.0.0；几天前才发布），则这些功能在jtools包中仍然可用。如果您确实更新到jtools 2.0.0，则需要遵循this link以获得有关如何下载interactions到CRAN的说明。

您的问题应该有一个简单的答案。 sim_slopes（“简单坡度”的缩写）功能可以为您提供所需的内容。

sim_slopes(cond_ranin_waterxsilver, pred = LnVolume, modx = WaterType, mod2 = SilverType, cond.int = TRUE)

这将打印出条件斜率和截距（仅当cond.int = TRUE时才打印截距。

如果需要使用这些值进行编程，则可以保存sim_slopes对象。

ss <- sim_slopes(cond_ranin_waterxsilver, pred = LnVolume, modx = WaterType, mod2 = SilverType, cond.int = TRUE)

ss$slopes # Matrix of slopes with test statistics, etc.
ss$ints # Matrix of intercepts with test statistics, etc.

绘制线性固定效果模型时，如何从jtools中的interact_plots中提取不同组的斜率和截距值

1 个答案: