Мы пытаемся смоделировать переменную count с избыточными нулями, используя пуассон с нулевым надуванием (как реализовано в пакете pscl).Вот (упрощенный) вывод, показывающий как категориальные, так и непрерывные пояснительные переменные:
library(pscl)
> m1 <- zeroinfl(y ~ treatment + some_covar, data = d, dist =
"poisson")
> summary(m1)
Count model coefficients (poisson with log link):
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.189253 0.102256 31.189 < 2e-16 ***
treatmentB -0.282478 0.107965 -2.616 0.00889 **
treatmentC 0.227633 0.103605 2.197 0.02801 *
some_covar 0.002190 0.002329 0.940 0.34706
Zero-inflation model coefficients (binomial with logit link):
Estimate Std. Error z value Pr(>|z|)
(Intercept) 0.67251 0.74961 0.897 0.3696
treatmentB -1.72728 0.89931 -1.921 0.0548 .
treatmentC -0.31761 0.77668 -0.409 0.6826
some_covar -0.03736 0.02684 -1.392 0.1640
резюме дало нам несколько хороших ответов, но мы ищем ANOVA-подобную таблицу.Итак, вопрос: можно ли использовать car :: Anova для получения такой таблицы?
> Anova(m1)
Analysis of Deviance Table (Type II tests)
Response: y
Df Chisq Pr(>Chisq)
treatment 2 30.7830 2.068e-07 ***
some_covar 1 0.8842 0.3471
Кажется, это работает нормально, но я не совсем уверен, является ли правильный подход, так как документация отсутствует (кажется, рассматривает только часть «модели подсчета»?).Вы рекомендуете следовать этому подходу или есть лучший способ?