glmmADMBを使用してゼロ膨張したNBモデル​​のゼロ膨張パラメータ（pz）をどのように取得しますか？

Question

私は、RのglmmADMBパッケージを使ってゼロ膨張の負の二項モデルを実行しました。私が理解するところでは、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