如何在R中运行截断后膨胀的Poisson模型？

Question

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

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

感谢您的帮助。

Answer 1

对于零充气模型或零障碍模型，标准方法是使用pscl包。我还写了一个适合这种型号here的包装，但尚未成熟并经过充分测试。除非您有大量数据，否则我仍然建议您使用更灵活，更可靠和更详尽的pscl。

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

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