<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />¿Hay una manera de especificar una distribución de poisson inflada a cero (estructura de error) a un modelo lineal generalizado en la plataforma del modelo de ajuste? Veo que puedo especificar distribuciones normales, binomiales, de poisson y exponenciales; pero lo que necesito es una distribución de poisson inflada en cero ... Cualquier ayuda muy apreciada. Gracias.
Clay Barker tiene un buen ejemplo de ZIPoisson en el intercambio de archivos JMP.
// configura la probabilidad de registro negativo ZIP con los valores iniciales definidos anteriormente
Exform = expr (
ll << configurar fórmula (Parámetro (
{p = expr (pstart), lambda = expr (lstart)},
Si (: Y == 0,
-Log (p + (1 - p) * Exp (-lambda)),
-Iniciar sesión(
(1 - p) * ((Exp (-lambda) * lambda ^: Y) /
Factorial (: y)
)
)
)));
// ahora inserte los valores iniciales en la columna de probabilidad de registro negativa
eval (eval Expr (exform));
// usar no lineal para minimizar la probabilidad de registro negativa, que es equivalente a maximizar la probabilidad
nl = No lineal (
Pérdida (: NegLogLik),
La pérdida es Neg LogLikelihood (1),
Newton,
Terminar
);
No conozco una manera de ajustar automáticamente un modelo ZI de ningún tipo. Siempre he ajustado dos modelos de partes, una parte para modelar 0 o no 0 y una segunda para modelar la parte no 0, luego combiné los modelos.