逸脱はGLMの概念であり、ZIPおよびZINBモデルはGLMではなく、GLMである分布の有限混合として定式化されるため、EMアルゴリズムを介して簡単に解くことができます。
これらのノートは、逸脱の理論を簡潔に説明しています。これらのメモを読むと、ポアソン回帰の飽和モデルに対数尤度があるという証拠が表示されます
ℓ(λs)=∑i=1,∀yi≠0n[yilog(yi)−yi−log(yi!)]
これは、プラグインがを推定した結果です。yi=λ^i
数学がより単純で、ZINBでも同様の結果が得られるため、ZIPの可能性について説明します。残念ながらZIPについては、ポアソンのような単純な関係はありません。番目の観測は、対数尤度がありますi
ℓi(ϕ,λ)=Zilog(ϕ+(1−ϕ)e−λ)+(1−Zi)[−λ+yilog(λ)−log(yi!)].
あなたがWRT偏微分を取る必要があると思い、これを解決することが観察されていないの両方と、0に数式を設定し、解決のためにと。ここでの問題は値です。これらはまたはに入る可能性があり、を観察せずに観測値を入れることはできません。ただし、値がわかっている場合は、欠損データがないため、ZIPモデルは必要ありません。観測されたデータは、EM形式の「完全なデータ」の可能性に対応しています。Ziλϕλϕyi=0λ^ϕ^Ziyi=0Zi
合理的かもしれない1つのアプローチは、を削除して期待値に置き換える、完全なデータログ尤度の期待値を使用することです。最新の更新でEMアルゴリズムが計算する内容(Eステップ)の一部。逸脱からこのアプローチを研究した文献はありません。ZiE(ℓi(ϕ,λ))Ziexpected
また、この質問が最初に尋ねられたので、この投稿に回答しました。ただし、同じトピックに関する別の質問があり、Gordon Smythによる素晴らしいコメントがここにあります。
ゼロ膨張複合ポアソンモデルの逸脱、連続データ(R)
で、同じ応答について言及しました(これは、私がコメントした詳細です)言う)加えて、他の投稿へのコメントで彼らがあなたが読みたいと思うかもしれない論文を述べた。(免責事項、私は参照された論文を読んでいません)