Simulación de procesos de Poisson no homogéneos y autoexcitados

Fueron mis modelos favoritos un tiempo, cuando modelaba visitas y revisitas de usuarios a cierto malhadado portal.

Si las visitas fuesen aleatorias (en cierto sentido), tendrían un aspecto no muy distinto del que se obtiene haciendo

Es decir,

o bien una distribución uniforme en el tiempo. Pero bien puede ocurrir que una visita incremente la probabilidad de otra inmediatamente después, por lo que las visitas tenderían a arracimarse en determinados momentos. Con el paquete IHSEP de R pueden simularse (y ajustarse) este tipo de modelos. Por ejemplo,

proporciona

que es más realista. La magia es obra de la función g, que incrementa la probabilidad de nuevos eventos después de que ocurra alguno, aunque la estela que dejan decae exponencialmente (por construcción).

Notas:

  • Los procesos anteriores tienen una base homogénea (la función de intensidad base se ha definido igual a .1), pero podría utilizarse una función de t arbitraria que los convertiría en procesos no homogéneos.
  • El paquete también proporciona mecanismos para estimar tanto la función base como la de autoexcitación. Que no sé qué tal funcionarán. Y si lo harán en tiempo y forma.

Comenta

Your email address will not be published.

Puedes usar estas etiquetas y atributos de HTML: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code class="" title="" data-url=""> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong> <pre class="" title="" data-url=""> <span class="" title="" data-url="">

Este sitio usa Akismet para reducir el spam. Aprende cómo se procesan los datos de tus comentarios.