廣義線性模型中的零膨脹泊松分佈<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />有沒有辦法在Fit Model平台中為廣義線性模型指定零膨脹泊松分佈(誤差結構)?我看到我可以指定正態,二項式,泊松和指數分佈;但我需要的是零膨脹的泊松分佈...... 任何幫助非常感謝。 謝謝。
Re:廣義線性模型中的零膨脹泊松分佈Clay Barker在JMP文件交換中有一個很好的ZIPoisson示例..
//使用上面定義的起始值設置ZIP負數對數似然
Exform = expr(
ll <<設定公式(參數(
{p = expr(pstart),lambda = expr(lstart)},
如果(:Y == 0,
-Log(p +(1-p)* Exp(-lambda)),
-Log(
(1 - p)*((Exp(-lambda)* lambda ^:Y)/
析因(:Y))
)
)
)));
//現在將起始值插入到負對數似然列中
eval(eval Expr(exform));
//使用非線性來最小化負對數似然,這相當於最大化可能性
nl =非線性(
損失(:NegLogLik),
損失是負LogLikelihood(1),
牛頓,
完
);
Re:廣義線性模型中的零膨脹泊松分佈我不知道如何自動適應任何類型的ZI模型。 我總是適合兩個部分模型,一部分模型0或不模型0,第二部分模型非0部分,然後組合模型。