零膨胀模型的“最小计数不是零”错误

Question

y为每2分钟内火车站站台的乘客人数，A1至A17为大堂17个研究区域的乘客人数。通过改变Xs，时间间隔已经得到了考虑。因为有时候，在大厅的学习区没有人等待，所以会出现过多的零。我计划使用零膨胀模型。我已经尝试了如图所示的代码，但它说“最小计数不是零”这意味着什么，我该如何解决它？我做了泊松，这没关系，但零膨胀不起作用。

    > setwd('C:/Users/zuzymelody/Desktop')
> try<-read.csv('0inflated_2mins27peak.csv',header=TRUE)
> attach(try)
> names(try)
 [1] "y"   "A1"  "A2"  "A3"  "A4"  "A5"  "A6"  "A7"  "A8"  "A9"  "A10" "A11"
[13] "A12" "A13" "A14" "A15" "A16" "A17"

    > model1<-glm(y~A1+A2+A3+A4+A5+A6+A7+A8+A9+A10+A11+A12+A13+A14+A15+A16+A17,family="poisson")
    > summary(model1)

Call:
glm(formula = y ~ A1 + A2 + A3 + A4 + A5 + A6 + A7 + A8 + A9 + 
    A10 + A11 + A12 + A13 + A14 + A15 + A16 + A17, family = "poisson")

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-7.8598  -3.4571  -0.3663   2.1867  12.5183  

Coefficients:
             Estimate Std. Error z value Pr(>|z|)    
(Intercept)  6.102009   0.164497  37.095  < 2e-16 ***
A1          -0.017555   0.003665  -4.790 1.66e-06 ***
A2          -0.026101   0.017569  -1.486 0.137371    
A3          -0.179988   0.014976 -12.018  < 2e-16 ***
A4          -0.032584   0.007735  -4.213 2.52e-05 ***
A5          -0.019908   0.007014  -2.839 0.004532 ** 
A6          -0.044144   0.010266  -4.300 1.71e-05 ***
A7           0.049829   0.006518   7.645 2.09e-14 ***
A8          -0.080712   0.009819  -8.220  < 2e-16 ***
A9           0.007390   0.007105   1.040 0.298273    
A10          0.041116   0.004085  10.065  < 2e-16 ***
A11         -0.041420   0.008418  -4.921 8.62e-07 ***
A12         -0.008241   0.007304  -1.128 0.259171    
A13         -0.033161   0.008966  -3.699 0.000217 ***
A14          0.020818   0.005250   3.965 7.34e-05 ***
A15         -0.002995   0.006125  -0.489 0.624887    
A16         -0.061997   0.017122  -3.621 0.000294 ***
A17         -0.025025   0.008391  -2.982 0.002860 ** 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for poisson family taken to be 1)

    Null deviance: 1137.71  on 29  degrees of freedom
Residual deviance:  599.74  on 12  degrees of freedom
AIC: 840.1

Number of Fisher Scoring iterations: 5







     >with(model1, cbind(res.deviance = deviance, df = df.residual,
      p = pchisq(deviance, df.residual, lower.tail=FALSE)))

 res.deviance df             p
[1,]     599.7445 12 1.202013e-120

> require( pscl )
> Zip<-zeroinfl(model1,link="logit",dist="poisson")
**Error in zeroinfl(model1, link = "logit", dist = "poisson") : 
  invalid dependent variable, minimum count is not zero**

dput（try）结构（列表（y=c（156L、74L、221L、207L、168L、36L、128L、208L、99L、117L、228L、211L、341L、173L、196L、310L、112L、203L、104L、183L、325L、143L、218L、166L、218L、127L、136L、38L、102L、34L），A1=c（24L、24L、24L、19L、20L、9L、14L、23L、15L、23L、14L、16L、15L、25L、19L、24L、26L、25L、26L、26L、22L、14L、13L、15L、9L、9L、12L、15L、15L、18L），A2=c（2L、4L、0L、1L、1L、2L、1L、2L、2L、2L、2L、0L、1L、1L、1L、1L、3L、2L、2L、3L、2L、3L、3L、3L、3L、3L、3L、3L、3L、3L、2L、3L、3L、3L、3L、3L、3L、3L、3L、3L、3L、3L、3L、3L、3L、3L、3L、3，A3=c（2L、2L、0L、1L、1L、9L、3L、0L、0L、0L、1L、1L、0L、0L、1L、2L、3L、1L、0L、1L、1L、0L、1L、1L、0L、0L、0L、2L）、A4=c（15L、11L、6L、7L、10L、5L、5L、5L、7L、9L、4L、6L、6L、13L、9L、13L、10L、10L、6L、6L、6L、6L、6L、6L、7L、6L、6L、6L、6L、6L、6L、7L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6，A5=c（13L、10L、6L、6L、11L、19L、13L、14L、7L、7L、6L、8L、10L、5L、7L、9L、9L、9L、9L、9L、9L、13L、8L、8L、8L、6L、8L、9L、14L、9L、4L、7L、7L、12L、11L、11L、12L、8L、6L、7L、7L、8L、8L、8L、8L、8L、8L、8L、6L、14L、9L、9L、9L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L，（21L，16L，16L，13L，13L，13L，13L，13L，13L，12L，12L，12L，12L，12L，12L，12L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，7L，6L，6L，6L，6L，6L，6L，6L，4L，6L，4L，6L，6L，4L，4L，4L，4L，5L，4L，5L，5L，5L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，6L，7L，7L，6L，7L，7L，7L，6L，7L，7L，7L，7L，7L，（8L、10L、6L、6L、10L、6L、6L、6L、6L、6L、0L、6L、6L、6L、6L、2L、2L、3L、3L、6L、6L、6L、6L、6L、6L、2L、2L、2L、2L、2L、6L、2L、2L、2L、6L、2L、2L、2L、2L、2L、3L、3L、3L、3L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、6L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、6L、6L、6L、6L、2L、6L、6L、2L、2L、2L、6L、6L、6L、，A12=c（12L、14L、14L、14L、17L、17L、17L、10L、14L、14L、17L、10L、14L、14L、13L、13L、13L、19L、19L、19L、7L、7L、7L、5L、5L、19L、19L、19L、19L、7L、5L、5L、5L、5L、5L、5L、5L、5L、5L、5L、5L、5L、5L、1L、1L、1L、1L、1L、1L、1L、1L、1L、1L、1L、1L、1L、1L、4L、4L、1L、1L、1L、4L、1L、4L、4L、1L、4L、1L、4L、4L、4L、4L、1L、4L、4L、4L、4L、4L，A14=c（14L、13L、16L、16L、11L、16L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、8L、8L、8L、8L、8L、6L、6L、6L、6L、6L、6L、8L、6L、6L、6L、6L、6L、6L、6L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L、11L，A16=c（2L、1L、3L、3L、3L、3L、3L、3L、3L、1L、1L、1L、2L、3L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、2L、3L、2L、2L、2L、3L、2L、3L、3L、2L、3L、2L、3L、2L、3L、3L、2L、3L、2L、3L、3L、2L、3L、3L、3L、2L、2L、3L、2L、3L、2L、3L、3L、2L、3L、3L、2L、3L、2L、2L、2L、2L、3L、2L、2L、2L、3L、3L，A17=c（10L、13L、13L、2L、5L、1L、3L、3L、5L、4L、4L、6L、4L、6L、3L、2L、2L、2L、2L、7L、8L、3L、7L、5L、6L、3L、4L、3L））。名称=c（“y”、“A1”、“A2”、“A3”、“A4”、“A5”、“A6”、“A7”、“A8”、“A9”、“A10”、“A11”、“A12”、“A13”、“A14”、“A15”、“A16”、“A17”），class=“data.frame”，行。名称=c（NA，-30L））

以上是可复制的示例。抱歉，这是我第一次在这里发帖，我不太懂规则

Answer 1

数据框的因变量$y$中不包含零值：

min（mydata$y）[1]34

你需要至少有一个$y=0$。