零膨胀模型R的二进制部分的响应曲线

Question

我已针对所有预测值绘制了每个预测变量的响应曲线，以确定每个预测变量如何影响我的计数。但是，我还想绘制零膨胀模型的二进制部分，以了解零膨胀模型的二进制部分中的预测变量如何帮助解释错误零的可能性。我正在尝试获得类似于下面链接页面底部的图，但是在该示例中，它们不提供可重现的代码。

https://fukamilab.github.io/BIO202/04-C-zero-data.html#sketch_fitted_and_predicted_values

我在下面添加了一些代码，这些代码具有零膨胀模型和使用的预测变量。然后，我使用predict函数来预测更大的栅格网格（new.data）的估计，并且我想查看这些预测值与我在整个栅格网格中使用的预测变量之间的响应。

mod1 = zeroinfl(Response~x1+x2|x1,link ="logit",data=data, 
                          dist="negbin")
modpred=predict(mod1, new.data, se.fit=T, type = "response")

response1 <- ggplot(data, aes(x = x1, y = modpred)) + geom_point()+
  +geom_smooth(data = data, aes(x = x1, y = modpred))