在 R (plot_model) 中编辑 sjplot 上的 y 轴

Question

我在创建的这个 sjplot 上的 y 轴有问题。我不确定为什么这些值是这样排列的（见图），并且想知道是否有人可以帮助我，例如将我的 y 轴设置为从 0 开始。

library(jtools)
library(carData)
library(effects)
library(sjPlot)
mod <- glmer(Golden.Trevally ~ Maturity.Status + Behavioural.Activity + (1 | ID.Number), family = "binomial", data = mydf2)
summary(mod)
plot_model(mod, "pred", title="")

Answer 1

到目前为止，回答这个问题最难的部分是重新创建数据以使其可重现。但是，以下内容非常接近：

library(jtools)
library(carData)
library(effects)
library(sjPlot)
library(lme4)

set.seed(69)

Behavioural.Activity <- factor(sample(c("Cleaning", "Courtship", 
                                        "Cruising", "Feeding"),
                                        size = 10000, 
                                        replace = TRUE))
Maturity.Status <- factor(sample(LETTERS[1:3], 10000, TRUE))
ID.Number <- factor(sample(500, 10000, TRUE))
Golden.Trevally <- rbinom(10000, 1, prob =
                          (c(6, 4, 7, 3)/600)[as.numeric(Behavioural.Activity)] *
                          c(0.8, 1, 1.2)[as.numeric(Maturity.Status)] *
                          (as.numeric(ID.Number) / 1000 + 0.75))

mydf2 <- data.frame(ID.Number, Golden.Trevally, 
                    Behavioural.Activity, Maturity.Status)

mod <- glmer(Golden.Trevally ~ Maturity.Status + Behavioural.Activity + (1 | ID.Number), 
             family = "binomial", data = mydf2)

my_sjplot <- plot_model(mod, "pred", title = "")

my_sjplot$Behavioural.Activity

这里的解决方案是实现plot_model返回的对象是一个包含两个ggplot对象的列表。您看到的是Behavioural.Activity。它看起来是这样的，因为它有一个scale_y_continuous，它的标签功能是将中断标记到最接近的百分比。你可以简单地用你自己的规模来覆盖这个规模：

my_sjplot$Behavioural.Activity +
   scale_y_continuous(limits = c(0, 0.01), 
                      labels = scales::percent_format(accuracy = 0.01))