r - Johnson-Neyman 情节适度调解

Question

我有一个相对简单的调节中介模型，看起来像这样（所有连续变量）：

当我过去进行调节分析时，我使用交互包创建了 Johnson-Neyman 图，以显示调节变量的值，在该值处，预测变量的斜率在不显着和显着之间变化。

我现在正在尝试创建与当前模型类似的东西。本文中的图 3表明可以在 mod-med 模型中创建这样的图，但我不确定他们是如何做到的。在过去的情节中，我指定了调节器和预测器并引用了线性模型（lm），但我认为这次的方法不正确......

任何指导都会非常有帮助。