如何使用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