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

Question

I would like to plot the estimates of the fixed effects of an lme4::lmer model. Currently, the sjPlot::plot_model output includes two additional plots that I do not want: SD (Observations) and SD (Intercept). How can I remove these two additional plots?

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

Essentially, I would like an output that looks like what is produced when plotting the fixed effects of an glm model, e.g.:

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

I believe that this what the sjPlot::plot_model function used to produce by default when plotting lmer models before because this is what my code used to produce. Thank you for your help!

sjp · Accepted Answer · 2021-11-23 13:29:33Z

2

Have you tried reinstalling all the packages and restarting your RStudio session? It's possible that you overwrote some sort of default because I can't reproduce the images that you have attached. I'm using the most recent version of sjPlot from CRAN (2.8.9).

This is what I get when I run your code, slightly modified to call the 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

And then I run the second part:

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

answered Nov 23, 2021 at 13:29

sjp

9101 gold badge6 silver badges11 bronze badges

Sign up to request clarification or add additional context in comments.

2 Comments

Sauterelle Over a year ago

Thank you very much! I had already tried restarting my R session several time but reinstalling the sjPlotpackage did the trick!

sjp Over a year ago

You're welcome, glad it's sorted!

Collectives™ on Stack Overflow

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

1 Answer 1

2 Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

2 Comments

Your Answer

Sign up or log in

Post as a guest

Related