将 R 中的 2 个预测概率图与 SjPlot 包相结合

Question

假设我有这个代码：

dt <- data.frame(x = (1:9), z = c(6,5,3,7,8,3,1,6,7), y = c(1,0,1,1,0,0,0,1,0))
model1 <- glm(y ~ x + z, family = binomial(link = "logit"),
              data = dt)
model2 <- glm(y ~ x + z, family = binomial(link = "logit"),
              data = dt)
summary(model1)
summary(model2)
library(sjPlot)
m <- plot_model(model1, title = "Predicted probabilities", type = "pred", terms = "x")
n <- plot_model(model2, title = "Predicted probabilities", type = "pred", terms = "x")

我想合并这两个图表，这样我就可以将它们保存在一个文件中。在这两个图中，您都可以从 sjPlot 包中找到预测概率。

我怎样才能做到这一点？

我已经尝试过par(mfrow=c(2,2))，但它根本不起作用。

任何帮助表示赞赏！

Answer 1

由于您尝试了par() 函数，我假设您想要并排创建两个图表。您可以使用库gridExtra 的grid.arrange() 函数实现此目的，如下所示：

使用的代码：

library(sjPlot)
library(gridExtra)

dt <- data.frame(x = (1:9), z = c(6,5,3,7,8,3,1,6,7), y = c(1,0,1,1,0,0,0,1,0))

model1 <- 0
model2 <- 0

model1 <- glm(y ~ x + z, family = binomial(link = "logit"),
              data = dt)
model2 <- glm(y ~ x^4 + z, family = binomial(link = "logit"),
              data = dt)
summary(model1)
summary(model2)

m <- plot_model(model1, title = "Predicted probabilities", type = "pred", terms = "x")
n <- plot_model(model2, title = "Predicted probabilities", type = "pred", terms = "x")

grid.arrange(m, n, ncol = 1, heights = c(1, 1))

如果您希望它们并排而不是一个在另一个之上，只需将 ncol 值更改为 2