R code for fitting a multitype NHPP model to grouped temporal data

The following R code demonstrates how to fit a multitype NonHomogeneous Poisson Process (NHPP) model to grouped temporal data. Multitype Poisson processes belong to the family of marked point processes.
Note that grouped temporal data occurs when only the number of recurrences within a given time interval are known.

The R code may be used to fit 1) a multitype NHPP model with a loglinear intensity function, with the intensity at time t defined by  e^{\gamma_{0} + \gamma_{1}t}, or 2) a multitype NHPP model with a power law intensity function, with the intensity at time t defined by  \gamma_{0}*t^{\gamma_{1}}.

Furthermore, the R code also shows how to fit a multitype Homogeneous Poisson Process (HPP) model to grouped temporal data.

NHPP multitype grouped

Continue reading R code for fitting a multitype NHPP model to grouped temporal data

R code for fitting a multitype NHPP model to temporal data

The following R code demonstrates how to fit a multitype NonHomogeneous Poisson Process (NHPP) model to temporal data. Multitype Poisson processes belong to the family of marked point processes.

The R code may be used to fit 1) a multitype NHPP model with a loglinear intensity function, with the intensity at time t defined by  e^{\gamma_{0} + \gamma_{1}t}, or 2) a multitype NHPP model with a power law intensity function, with the intensity at time t defined by  \gamma_{0}*t^{\gamma_{1}}.

Furthermore, the R code also shows how to fit a multitype Homogeneous Poisson Process (HPP) model to temporal data.

NHPP multitype

Continue reading R code for fitting a multitype NHPP model to temporal data

R code for fitting a nonhomogeneous temporal Poisson process model using the spatstat package

In my previous blog post I showed how to fit power and loglinear nonhomogeneous Poisson process (NHPP) models to temporal data.

The following R code fits these same two models, but this time using the spatstat package.

NHPP spatstat

Continue reading R code for fitting a nonhomogeneous temporal Poisson process model using the spatstat package

R code for fitting a nonhomogeneous temporal Poisson process model

The following R code demonstrates how to fit a nonhomogeneous Poisson process (NHPP) model to temporal data.

The R code may be used to fit 1) a NHPP model with a loglinear intensity function, with the intensity at time t defined by e^{\gamma_{0} + \gamma_{1}t} , or 2) a NHPP model with a power law intensity function, with the intensity at time t defined by \frac{\beta}{\eta}(\frac{t}{\eta})^{\beta-1}. The latter NHPP model is sometimes referred to as the Crow-AMSAA model.

NHPP loglinear

Continue reading R code for fitting a nonhomogeneous temporal Poisson process model

R code for fitting a loglinear-NHPP model to grouped temporal data

The following R code fits a NonHomogeneous Poisson Process (NHPP) model to grouped temporal data. Grouped temporal data occurs when only the number of recurrences within a given time interval are known.
The R code fits a NHPP model with a loglinear intensity function, where the intensity at time t is given by e^{\gamma_{0} + \gamma_{1}t} .
NHPP loglikelihood Continue reading R code for fitting a loglinear-NHPP model to grouped temporal data