如何使用glmmADMB获得零膨胀NB模型的零膨胀参数（pz）？

Question

我使用R中的glmmADMB软件包运行了一个零膨胀负二项模型。据我所知，pz参数是零膨胀参数，它由软件包拟合到你运行的模型 - pz搜索最适合您数据的值，包从pz = 0.2开始搜索。这是默认值，可以更改。

运行模型后，是否有人知道如何找到为数据选择的pz值？

Answer 1

可以从模型对象获得零膨胀估计（连同其标准偏差）。请参阅下面的glmmADMB包中的内置数据：

library(glmmADMB)
# munge data
Owls = transform(Owls, Nest = reorder(Nest, NegPerChick), 
                 logBroodSize = log(BroodSize), NCalls = SiblingNegotiation)
# fit model
fit_zinb = glmmadmb(NCalls ~ (FoodTreatment + ArrivalTime) * SexParent +
                             offset(logBroodSize) + (1 | Nest), 
                    data = Owls, zeroInflation = TRUE, 
                    family = "nbinom")
# overall summary, check for match
summary(fit_zinb)
# zero-inflation estimate
fit_zinb$pz
# zero-inflation standard deviation
fit_zinb$sd_pz