我使用R中的glmmADMB软件包运行了一个零膨胀负二项模型。据我所知,pz参数是零膨胀参数,它由软件包拟合到你运行的模型 - pz搜索最适合您数据的值,包从pz = 0.2开始搜索。这是默认值,可以更改。
运行模型后,是否有人知道如何找到为数据选择的pz值?
答案 0 :(得分:0)
可以从模型对象获得零膨胀估计(连同其标准偏差)。请参阅下面的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