1

我今天的第二个问题：我想使用这个visreg包来绘制我的逻辑回归模型。只要我不使用属性“by”，它就像一个魅力，但是当我想使用它时，我得到一个错误。我用来创建模型的代码如下：

m3<- glm(alive ~ seatbelt*dvcat + sex + ageOFocc + airbag, family = binomial, data = nassCDS, start=)
summary(m3)

如果我然后使用：

visreg(m3, "seatbelt", scale = "response")

我得到以下结果 在此处输入图像描述

这很好。但是，如果我现在添加“by”属性，则会出现错误：

visreg(m3, "seatbelt", by="dvcat", scale ="response")

在此处输入图像描述

我用谷歌搜索，据我了解，该函数找不到绘制模型的数据。但是我在哪里可以提供数据？我已经尝试过“data=”属性，但它对我不起作用（或者我做错了）。没有控制台输出，我只能提供图表本身的消息。有人可以帮助我吗？亲切的问候，简 :)

编辑：我使用了来自 vincent arel-bundocks github 的“nassCDS”数据卫星，您可以在此处找到：https : //vincentarelbundock.github.io/Rdatasets/datasets.html能够使用逻辑回归。因此，我使用了 dplyr 包和以下代码：

nassCDS <- nassCDS %>% 
  mutate(dead1 = as.integer(dead)) %>%
  mutate(alive = sjmisc::rec(dead1, rec = "2=0; 1=1")) %>%
select(seatbelt, dead, alive, dvcat, sex, ageOFocc, everything()) %>% 
  select(-dead1)

此外，我将列安全气囊和安全带更改为数字，正如另一位 stackoverflow 用户所建议的那样。

4

0 回答 0