使用包visreg绘制物流模型不工作未找到数据

使用包visreg绘制物流模型不工作未找到数据,r,plot,graph,visreg,R,Plot,Graph,Visreg,我今天的第二个问题:我想用visreg软件包来绘制我的逻辑回归模型。只要我不使用属性,它就像一个符咒,但当我想使用它,我会得到一个错误。我用于创建模型的代码如下所示: m3<- glm(alive ~ seatbelt*dvcat + sex + ageOFocc + airbag, family = binomial, data = nassCDS, start=) summary(m3) 我得到以下结果 这很好。但如果现在添加by属性,则会出现错误: visreg(m3, &quo

我今天的第二个问题:我想用visreg软件包来绘制我的逻辑回归模型。只要我不使用属性,它就像一个符咒,但当我想使用它,我会得到一个错误。我用于创建模型的代码如下所示:

m3<- glm(alive ~ seatbelt*dvcat + sex + ageOFocc + airbag, family = binomial, data = nassCDS, start=)
summary(m3)
我得到以下结果

这很好。但如果现在添加by属性,则会出现错误:

visreg(m3, "seatbelt", by="dvcat", scale ="response")
我在谷歌上搜索过,据我所知,函数无法找到绘制模型的数据。但是我在哪里可以提供数据呢?我已经尝试了data=属性,但它对我不起作用,或者我做错了。没有控制台输出,我只能提供图形本身上的消息。有人能帮我吗?问候你,简:

编辑:我使用了vincent arel bundocks github提供的nassCDS数据,您可以在这里找到:我刚刚通过列dead插入了列alive,因此我可以使用逻辑回归。因此,我使用了带有以下代码的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用户的建议,我将安全气囊和安全带列更改为nummeric。

NASCDS是某个软件包附带的默认数据集吗?如果没有,则通过复制并粘贴DPUTNASCDS的输出来编辑问题。数据集是以前默认包的一部分,但现在不再是了。我想,函数dput的输出太长,无法复制它?我从加拿大一位教授的github获得了数据集:您只需搜索nassCDS数据集不包含活动列,而是包含死列,并且在您的代码中缺少start=。@BappaDas我编辑了上面的问题,提供了我用于创建活动列的代码。nassCDS是随某些包提供的默认数据集吗?如果没有,则通过复制并粘贴DPUTNASCDS的输出来编辑问题。数据集是以前默认包的一部分,但现在不再是了。我想,函数dput的输出太长,无法复制它?我从加拿大一位教授的github获得了数据集:您只需搜索NASSCDS数据集不包含活动列,而是包含死列,并且在您的代码中缺少start=。@BappaDas我编辑了上面的问题,提供了用于创建活动列的代码。
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)