我正在拟合零通胀Poisson回归模型,它是摘要输出:
Call:
zeroinfl(formula = m ~ lcounts | 1, dist = "poisson")
Pearson residuals:
Min 1Q Median 3Q Max
-22.090 -13.297 -4.686 26.634 115088.430
Count model coefficients (poisson with log link):
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.788e+00 5.901e-04 6418.7 <2e-16 ***
lcountsf 5.306e-01 7.644e-05 6941.9 <2e-16 ***
lcountsm 1.095e-02 4.031e-05 271.6 <2e-16 ***
lcountso 1.999e-01 8.525e-05 2345.0 <2e-16 ***
lcountss 7.289e-02 4.432e-05 1644.8 <2e-16 ***
lcountst 3.916e-02 4.258e-05 919.6 <2e-16 ***
Zero-inflation model coefficients (binomial with logit link):
Estimate Std. Error z value Pr(>|z|)
(Intercept) -6.6645 0.4011 -16.62 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Number of iterations in BFGS optimization: 14
Log-likelihood: -1.965e+07 on 7 Df
我在徘徊如何使用这种模型进行预测?
正常的glm
就是predict(model, newdata)
。
有专门的命令可以做到吗?