Как вы получаете параметр нулевой инфляции (pz) нулевой вздутой модели NB с использованием glmmADMB?

Question

Я запустил отрицательную биномиальную модель с нулевым раздувом, используя пакет glmmADMB в R. Из того, что я понимаю, параметр 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