不再可以将 vglm 用于 R 中的欠分散计数数据?
No longer can use vglm for underdispersed count data in R?
我在别处读到,VGAM 包可用于通过 genpoisson 系列对欠分散的计数数据进行建模。但是,当我查找 genpoisson0、genpoisson1 和 genpoisson2 的帮助文件时,他们都说如下:
“理论上,允许 λ 参数为负数以处理欠分散,但不再支持该参数,因此 0 < λ < 1。”
“理论上 \varphi 参数可能被允许小于 1 以处理欠分散,但这不受支持。”
“理论上 α 参数可能被允许为负数以处理欠分散,但这是不受支持的。”
我在哪里可以处理欠分散?
- 您可以使用准似然方法,例如
family=quasipoisson
在 glm()
中(在基础 R 中)
glmmTMB
包支持 COM-Poisson (family = compois
) 和广义泊松 (family = genpois
) 条件分布。
我不清楚 VGAM
中不再支持分散不足的广义泊松分布的原因(简要地)讨论 here 是否也适用于 glmmTMB
中的实现。 ..
Brooks et al (2019) 中对 glmmTMB
parameterizations/implementations 的 COM-Poisson 和广义泊松进行了一些讨论。
Brooks、Mollie E.、Kasper Kristensen、Maria Rosa Darrigo、Paulo Rubim、María Uriarte、Emilio Bruna 和 Benjamin M. Bolker。 “年生育率模式的统计模型。”生态 100,没有。 7 (2019): e02706。 https://doi.org/10.1002/ecy.2706.
我在别处读到,VGAM 包可用于通过 genpoisson 系列对欠分散的计数数据进行建模。但是,当我查找 genpoisson0、genpoisson1 和 genpoisson2 的帮助文件时,他们都说如下:
“理论上,允许 λ 参数为负数以处理欠分散,但不再支持该参数,因此 0 < λ < 1。”
“理论上 \varphi 参数可能被允许小于 1 以处理欠分散,但这不受支持。”
“理论上 α 参数可能被允许为负数以处理欠分散,但这是不受支持的。”
我在哪里可以处理欠分散?
- 您可以使用准似然方法,例如
family=quasipoisson
在glm()
中(在基础 R 中) glmmTMB
包支持 COM-Poisson (family = compois
) 和广义泊松 (family = genpois
) 条件分布。
我不清楚 VGAM
中不再支持分散不足的广义泊松分布的原因(简要地)讨论 here 是否也适用于 glmmTMB
中的实现。 ..
Brooks et al (2019) 中对 glmmTMB
parameterizations/implementations 的 COM-Poisson 和广义泊松进行了一些讨论。
Brooks、Mollie E.、Kasper Kristensen、Maria Rosa Darrigo、Paulo Rubim、María Uriarte、Emilio Bruna 和 Benjamin M. Bolker。 “年生育率模式的统计模型。”生态 100,没有。 7 (2019): e02706。 https://doi.org/10.1002/ecy.2706.