零膨胀负二项分布

Question

我需要 R 中以下任务的帮助： 从零膨胀负二项分布生成 10000 个样本（每个样本有 30 个观测值）。 我找到了函数

rzinbinom()

，但无法计算出要输入的参数值；我们没有被分配任务。我们可以输入默认值来使函数正常工作吗？

Answer 1

零膨胀负二项式有三个参数：0 的概率，然后是两个“经典”负二项式参数。有大约十几种方法来参数化标准负二项式，这会导致混乱。在 R 中，这两者通常是均值和

size

，后者相当于 Klugman-Panger-Willmot 中的“r”，或者如果将 NB 视为具有 gamma 先验的泊松分布，则相当于混合 gamma 的 alpha。

这意味着，即使知道每个样本的预期值（我认为“每个样本有 30 个观察值”）也不足以准确确定要使用哪个 ZINBINOM。不同的 pi 和大小确实有无数的可能性，平均为 30。要么你的老师希望你选择任何一组有效的参数，要么你的指示是不够的。