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