现在就遇到这个。为避免混淆,我是原始问题中引用的Wilson(2015)的Wilson,该问题询问Poisson模型和截短的Poisson模型是嵌套的,非嵌套的等。稍微简化一下,如果较大的模型嵌套在较大的模型中,如果模型的参数子集固定为指定值,则模型缩小为较小的模型;如果两个模型在各自参数的子集固定为特定值时都归结为同一模型,则它们是重叠的,如果参数如何固定,一个模型都不能归结为另一个模型,则这两个模型是非嵌套的。根据此定义,截短的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")
vuong
中的软件包pscl
中的函数,该函数说这是针对非嵌套模型的。我只是谷歌搜索并找到包含参数“嵌套”的vuongtest
包nonnest2
中的函数。是吗