Я построил кривые отклика для каждого из моих предикторов по всем прогнозируемым значениям, чтобы определить, как каждый предиктор влияет на мои показатели.Однако я также хочу построить бинарную часть моей модели с нулевым раздувом, чтобы увидеть, как предикторы в двоичной части модели с нулевым раздувом помогают объяснить вероятность ложных нулей.Я пытаюсь получить график, похожий на график внизу страницы по ссылке ниже, однако в этом примере они не предоставляют воспроизводимый код.
https://fukamilab.github.io/BIO202/04-C-zero-data.html#sketch_fitted_and_predicted_values
Я включил ниже некоторый код, где у меня есть модель с нулевым раздуванием и используемые предикторы.Затем я использую функцию прогнозирования, чтобы прогнозировать оценки для гораздо большей растровой сетки (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))