过度离散的次数模型:拟泊松、负二项、零膨胀和 hurdle

出险次数(或门诊次数)的方差比均值大很多、0 特别多时,泊松模型错在哪里?

SOA SRMSOA PACAS MAS-I
对过的标准答案
3 条
覆盖的考试公式
1 个
参数
8 个
开源库
2 个

01当场运行

改一个参数,当场重算

正在用默认参数运行,计算量大的方法最多要 20 秒。

04出处和说明

公式从哪里来,有什么限制

列出参考的教材、论文和考试公式条目。方法的限制和实现上的说明也写在这里。

参考

  • Cameron、Trivedi:Regression Analysis of Count Data第 3–4 章
  • Frees:Regression Modeling with Actuarial and Financial Applications第 12 章 Count dependent variables
  • Cameron、Trivedi (1990) Regression-based tests for overdispersion in the Poisson model

覆盖的考试公式

  • REG-029过度离散:拟似然、负二项、零膨胀与 hurdle 模型SRMPAMAS-I

公式条目按精算考试的公式表整理,编号和中文名是平台自己的。

这个方法算什么

对同一份次数数据拟合泊松、拟泊松、负二项 NB2、零膨胀泊松(ZIP)和 hurdle 泊松,比较对数似然、AIC、BIC 和预测的 0 次的个数;做 Pearson 离散度、Cameron-Trivedi 检验和负二项对泊松的边界似然比检验;列出选定模型的系数和各次数的实际/预测频数。

说明和限制

  • 拟泊松的系数和泊松相同,标准误乘 √φ̂。它没有似然,不参加 AIC 的比较。
  • 负二项 NB2 对 ln α 做一维最大似然,给定 α 时用 IRLS 解 β。α → 0 时退化为泊松,所以似然比检验的原假设在边界上,p 值取 χ²₁ 尾概率的一半。
  • ZIP 和 hurdle 的区别:ZIP 的 0 来自两处(结构性的 0 和泊松的 0),hurdle 的 0 只来自跨栏的那一部分。截距模型时两者的似然和预测的 0 的个数相同。
  • Cameron-Trivedi 检验:((y − μ̂)² − y)/μ̂ 对 μ̂ 做不带截距的回归,系数的 t 检验 α = 0。
  • RAND 数据取原数据每 10 行的 1 行(2,019 人),控制运行时间。

数据

看参数里选的数据:选平台的数据集时是真实数据或公开样例,填自己的数时只用输入的数。上面运行的输出里写明这一次用的是哪种。

开源库和版本

  • scipy1.17.1BSD分布拟合、统计检验
  • numpy2.4.6BSD-3-Clause数值计算

05同族的方法

广义线性模型,还有 21 个方法

同一个方法族里的其他方法。在智能体里,用一句话就能调用这些方法。

在智能体里这样问“出险次数(或门诊次数)的方差比均值大很多、0 特别多时,泊松模型错在哪里?拟泊松、负二项、零膨胀泊松、hurdle 模型哪个更合适?”