一般化線形モデルにおけるゼロ膨張ポアソン分布<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Fit Modelプラットフォームで一般化線形モデルにゼロ膨張ポアソン分布(誤差構造)を指定する方法はありますか?正規分布、二項分布、ポアソン分布、指数分布を指定できます。しかし、私が必要としているのはゼロ膨張ポアソン分布です。 任意の助けは大歓迎です。 ありがとう。
一般化線形モデルにおけるRe:ゼロ膨張ポアソン分布Clay BarkerはJMPファイル交換に関する良いZIPoissonの例を持っています。
//上で定義した開始値を使用して、ZIPの負の対数尤度を設定します
Exform = expr(
<<式の設定(パラメータ(
{p = expr(pstart)、λ= expr(lstart)}、
(:Y == 0の場合
-Log(p +(1 - p)* Exp( - λ))、
-ログ(
(1 - p)*((Exp( - λ)*λ^:Y)/
階乗(:Y))
)
)
)))
//負の対数尤度列に開始値を代入します
eval(eval Expr(exform));
//負の対数尤度を最小化するために非線形を使用します。これは尤度を最大化することと同じです。
nl =非線形(
損失(:NegLogLik)、
損失は負の対数尤度(1)です。
ニュートン、
終了
;
一般化線形モデルにおけるRe:ゼロ膨張ポアソン分布どのタイプのZIモデルにも自動的にフィットする方法がわかりません。 私はいつも2つの部品モデル、1つの部品は0または0以外をモデル化し、もう1つは0以外の部分をモデル化してからモデルを結合しました。