现在才遇到这个。为了避免混淆,我是原始问题中引用的 Wilson of Wilson(2015),它询问 Poisson 和截断 Poisson 模型是否是嵌套的、非嵌套的等。稍微简化一下,如果较大的模型,较小的模型将嵌套在较大的模型中如果其参数的子集固定在规定的值,则模型会缩减为较小的模型;如果两个模型在各自参数的子集固定为特定值时都归结为同一模型,则它们是重叠的,如果无论参数如何固定,一个都不能归结为另一个,则它们是非嵌套的。根据这个定义,截断泊松和标准泊松是非嵌套的。然而,这一点似乎被许多人忽视了,Vuong 的分配理论指的是严格嵌套,严格非嵌套,并且严格重叠。“严格”指的是在嵌套等的基本定义中添加了六个限制。这些限制并不完全简单,但它们确实意味着 Vuong 关于对数似然比分布的结果不适用于以下情况模型/分布嵌套在参数空间的边界(就像泊松/零膨胀泊松与零膨胀参数的恒等链接的情况一样),或者当一个模型倾向于另一个模型时,当参数趋于无穷时,如当使用 logit 链接对零膨胀参数建模时,泊松/零膨胀泊松就是这种情况。Vuong 没有提出关于在这些情况下对数似然比分布的理论。不幸的是,这里
以下 R 代码将模拟泊松分布和截断泊松对数似然比。它需要VGAM软件包。
n<-30
lambda1<-1
H<-rep(999,10000)
for(i in 1:10000){
print(i)
y<-rpospois(n, lambda1)
fit1 <- vglm(y ~ 1, pospoisson)
fit2<-glm(y~1, family=poisson(link="log"))
H[i]<-logLik(fit1)-logLik(fit2)
}
hist(H,col="lemonchiffon")