首页 > 解决方案 > 在 ggplot 中绘制混合模型

问题描述

我构建了一个混合模型,现在使用sjp.lmer,我正在根据我的每个预测器绘制预测值。拟合模型的名称是mdl.lmer,并且names.vec是预测变量名称的向量。我想做一个这样的循环,以便所有图都出现在一个 pdf 中

pdf("my.pdf")
for(i in seq_along(names.vec)){

    clim <- names.vec[i]
    sjp.lmer(mdl.lmer, type = "pred.fe", vars = clim, show.scatter = T)
    ggsave(plot = last_plot())
}

问题是 sjp.lmer 线制作了一个绘图,但我无法将其分配给一个对象,然后将该对象打印到一个打开的设备。通常我这样做:

pdf("my.pdf")
for(i in seq_along(names.vec)){

  clim <- names.vec[i]

  p <- ggplot(dat, aes(x = clim, y = response)) + geom_line()
  print(p)
}

此命令将所有绘图保存在打开的 pdf 中。但是,对于混合模型案例,我无法实现这一点。

标签: rggplot2lme4mixed-models

解决方案


推荐阅读