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

【问题标题】：How to run a truncated and inflated Poisson model in R?如何在 R 中运行截断和膨胀的泊松模型？
【发布时间】：2020-07-19 09:37:06
【问题描述】：

我的数据不包含任何零。我的结果 y 的最小值是 1，这就是被夸大的值。我的目标是使用 R 运行截断和膨胀的泊松回归模型。

我已经知道如何区分每个回归零截断和零膨胀的方式。我想知道如何将这两个条件组合成一个模型。

感谢您的帮助。

【问题讨论】：

欢迎堆栈溢出。请提供尽可能多的详细信息，以最大限度地提高您获得帮助的机会。你能进一步解释一下吗？参考可能会有所帮助。我会推荐tweedie distribution 和tweedie。 R 会给你一种实现它的方法。在crossvalidate 上您可能会有更好的机会。

【解决方案1】：

对于零膨胀模型或零障碍模型，标准方法是使用pscl 包。我还写了一个适合这种模型的包here，但它还没有成熟和完全测试。除非您拥有大量数据，否则我仍然建议您使用更灵活、更强大且有文档的pscl。

对于零截断模型，您可以查看VGML::vglm 函数。您可能会找到有用的信息here。

请注意，您没有进行相同的分布假设，因此您不需要相同的估计数据。鉴于您的数据集的描述，我认为您正在寻找一个零截断模型（因为您没有观察到零）。使用零膨胀模型，您可以将观察到的模式分解为由选择模型生成的零和由计数数据模型生成的其他零。这看起来与您的数据集不一致。

【讨论】：

感谢您的回答。但实际上我想获得零截断和一个膨胀模型的信息。在上传到问题之前，我已经检查了各种软件包，例如 pscl、vglm、glmmTMB 等。但我找不到这个模型的方法（“零截断和膨胀模型）。我希望你有其他想法。