如何使用 sjPlot::plot_model 仅显示 lmer 模型的固定效应估计

How to only show fixed effect estimates of lmer model using sjPlot::plot_model

我想绘制 lme4::lmer 模型的固定效应估计值。目前,sjPlot::plot_model 输出包括两个我不想要的附加图:SD(观察)和 SD(截距)。我怎样才能删除这两个额外的地块?

library(lme4)
data(iris)
lmerMod <- lmer(Sepal.Length ~ Sepal.Width * Petal.Length + (1|Species), data=iris)

library(sjPlot)
plot1 <- plot_model(lmerMod,
           type = "est",
           show.intercept = TRUE,
           title = "Estimates of fixed effects")

plot1

本质上,我想要一个看起来像绘制 glm 模型的固定效应时产生的输出,例如:

glmMod <- glm(Sepal.Length ~ Sepal.Width * Petal.Length, data=iris)
plot2 <- plot_model(glmMod,
           type = "est",
           show.intercept = TRUE,
           title = "Estimates of fixed effects")

plot2

我相信这是 sjPlot::plot_model 函数在之前绘制 lmer 模型时默认生成的内容,因为这是我的代码用来生成的内容。感谢您的帮助!

您是否尝试过重新安装所有软件包并重新启动 RStudio 会话?您可能覆盖了某种默认值,因为我无法复制您附加的图像。我正在使用来自 CRAN (2.8.9) 的最新版本的 sjPlot。

这是我 运行 你的代码时得到的,稍微修改以调用 plot to plot window:

library(lme4)
data(iris)
lmerMod <- lmer(Sepal.Length ~ Sepal.Width * Petal.Length + (1|Species), data=iris)

library(sjPlot)
plot1 <- plot_model(lmerMod,
           type = "est",
           show.intercept = TRUE,
           title = "Estimates of fixed effects")
plot1

然后我运行第二部分:

glmMod <- glm(Sepal.Length ~ Sepal.Width * Petal.Length, data=iris)
plot2 <- plot_model(glmMod,
                    type = "est",
                    show.intercept = TRUE,
                    title = "Estimates of fixed effects")
plot2