在R中的多级模型中绘制二次效应

Question

我计算了以下模型：

a.Pfad<-lmer(FlowZ~MTsqZ+(1|VP04_01), data=MT)

MTsqZ是二次项。

我正在寻找一种绘制这个模型的方法：我想知道变量之间是否存在倒U形关系。 做这个的最好方式是什么？由于我对R中的任何图表都没有经验，我真的很期待你的帮助。

Answer 1

您可以查看sjPlot包来绘制lmer型号的效果尺寸。您想要的功能是sjp.lmer，如果您想绘制MTsqZ的固定效果斜率，您可以输入：

library(sjPlot)
sjp.lmer(a.Pfad, type = "fe.slope", vars = "MTsqZ")

附加评论：从它的外观来看，您没有将MT和FlowZ之间的关系建模为二次方，因为您只包含变量的平方版本。

编辑：这就是如何将关系建模为二次方：

a.Pfad<-lmer(FlowZ~MT+I(MT^2)+(1|VP04_01), data=MT)

然后可以使用type = "poly"绘制关系：

sjp.lmer(a.Pfad, type = "poly", poly.term = "MT")