How to model a spatial point process with INLA?

  • We fit the intensity by modelling the number of points per unit area as a Poisson distribution.
  • The area offset for the mesh vertices (integration points) is determined by a dual mesh: polygons around the primary mesh vertices. In effect, the offset serves as a weight.
    • The response at the mesh vertices is set to an initial observation of 0 points and offset according to the dual mesh.
    • The response at the event locations is set as 1 with an offset of 0.