我有一个平均值 = 3.2 的计数数据集，并且有点零膨胀。

X1      X2      X3  Y
Food3   Low     13  2
Food3   High    27  1
Food2   Low     13  1
Food1   Medium  27  1
Food1   High    20  8
Food3   Low     20  1
Food1   High    13  5
Food2   Medium  13  4
Food1   Low     13  0
Food2   High    20  6
Food1   Medium  13  2
Food1   Low     13  1
Food1   Low     13  1
Food3   Low     13  1
Food2   Medium  13  5
Food1   Medium  27  0
Food3   Low     13  2
Food1   Medium  20  3
Food3   Medium  13  7
Food1   Low     20  1
Food3   Medium  13  5

我用泊松族拟合了 GLM 模型：

model1 <- glm(formula=Y~X1+X2+X3+X1:X2+X1:X3+X2:X3, 
              family=poisson(link="log"), data=Df)

summary(model1)输出显示有点过度分散，我还尝试拟合 glm.nb() 负二项式 GLM。

但是这个模型的问题是，对于泊松 GLM 和负二项式都有一些负面预测。这些是怎么来的，我应该如何解决这个问题？