<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Existe-t-il un moyen de spécifier une distribution de Poisson (structure d'erreur) à zéro gonflé dans un modèle linéaire généralisé dans la plateforme de modèle d'ajustement? Je vois que je peux spécifier des distributions normales, binomiales, de poisson et exponentielles; mais ce dont j'ai besoin, c'est d'une distribution de poisson à gonflage nul ... Toute aide grandement appréciée. Merci.
Clay Barker a un bel exemple ZIPoisson sur l'échange de fichiers JMP.
// configure la probabilité de journal négatif ZIP avec les valeurs de départ définies ci-dessus
Exform = expr (
ll << définir la formule (Parameter (
{p = expr (pstart), lambda = expr (lstart)},
Si (: Y == 0,
-Log (p + (1 - p) * Exp (-lambda)),
-Bûche(
(1 - p) * ((Exp (-lambda) * lambda ^: Y) /
Factorielle (: Y))
)
)
)));
// branche maintenant les valeurs de départ dans la colonne de probabilité de journal négatif
eval (eval Expr (exform));
// utilise non linéaire pour minimiser la probabilité de log négative, ce qui équivaut à maximiser la probabilité
nl = non linéaire (
Perte (: NegLogLik),
Loss is Neg LogLikelihood (1),
Newton,
terminer
)
Je ne connais pas de moyen d'adapter un modèle ZI de tout type automatiquement. J'ai toujours adapté deux modèles de pièces, une partie au modèle 0 ou non 0 et une seconde pour modéliser la partie non 0, puis j'ai combiné les modèles.