如何在R中运行截断和膨胀的泊松模型？

install.packages("truncdist")
install.packages("pscl")

library(truncdist)
library(pscl)

trunc_poisson <- function(x, lambda) {
  return(dtrunc(x, spec="pois", a=0, b=Inf, mean=lambda))
}

inflated_poisson <- function(y, formula) {
  return(zeroinfl(y ~ formula, dist="poisson"))
}

# 创建截断泊松分布对象
trunc_dist <- trunc_poisson(x, lambda)

# 创建膨胀泊松模型对象
inflated_model <- inflated_poisson(y, formula)

# 模型拟合
fit <- fit(inflated_model)

# 预测
predict <- predict(fit, newdata=data)