<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Gibt es eine Möglichkeit, eine Null-aufgeblasene Poisson-Verteilung (Fehlerstruktur) für ein verallgemeinertes lineares Modell in der Fit Model-Plattform anzugeben? Ich sehe, ich kann normale, Binomial-, Poisson- und Exponentialverteilungen angeben; aber was ich brauche ist eine null-aufgeblasene poisson-verteilung ... Jede Hilfe sehr geschätzt. Vielen Dank.
Clay Barker hat ein schönes ZIPoisson-Beispiel für den JMP-Dateiaustausch.
// Die ZIP-Negativprotokollwahrscheinlichkeit mit den oben definierten Startwerten einrichten
Exform = Ausdruck (
ll << Formel setzen (Parameter (
{p = Ausdruck (pstart), Lambda = Ausdruck (lstart)},
Wenn (: Y == 0,
-Log (p + (1 - p) * Exp (-Lambda)),
-Log(
(1 - p) * ((Exp (-lambda) * lambda ^: Y) /
Fakultät (: Y))
)
)
)));
// Stecken Sie jetzt die Startwerte in die Spalte für negative Log-Likelihood
eval (eval Expr (exform));
// Verwenden Sie nichtlinear, um die negative Log-Wahrscheinlichkeit zu minimieren. Dies entspricht der Maximierung der Wahrscheinlichkeit
nl = nichtlinear (
Verlust (: NegLogLik),
Verlust ist Neg LogLikelihood (1),
Newton,
Fertig
);
Ich weiß nicht, wie man ein ZI-Modell eines beliebigen Typs automatisch anpasst. Ich habe immer zwei Teilmodelle montiert, ein Teil für Modell 0 oder nicht 0 und ein zweites Modell, das nicht 0 ist, und dann die Modelle kombiniert.