我正在使用计数模型对零截断过程进行建模,并试图确定数据是否过度分散。泊松分布的方差等于其均值, 负二项式模型通过估计过度离散参数。这种方法有两种方案, 和 我经常看到第二个版本。我估计了三个具有相同协变量但具有三个不同链接函数(、、)的模型。
NB1
NB2
Poisson
NB1
NB2
Model logLik alpha p(alpha)
Poisson -7942 1 NA
NB1 -6399 1.001 0
NB2 -7944 403.4 0
显然,过度分散参数很重要,因此我应该使用其中一个NB
模型。但是,当模型中有附加参数时,对数似然怎么可能NB2
低于呢?Poisson
当我估计所有三个模型时,我得到一个警告,我的方差-协方差矩阵可能不是半正定的。但是我没有高度共线的变量,并且所有相关矩阵都有严格的正特征值,所以我不明白为什么我会收到这个警告。这可能与可能性问题有关吗?
我正在使用同名包glmmadmb
中的函数。R