我计算了以下模型:
a.Pfad<-lmer(FlowZ~MTsqZ+(1|VP04_01), data=MT)
与MTsqZ
是一个二次项。
我正在寻找一种绘制这个模型的方法:我想知道变量之间是否存在倒U形关系。做这个的最好方式是什么?由于我对R中的任何图表都没有经验,我真的很期待你的帮助。
您可以查看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")