est_intens_spde {highriskzone}

Estimates the intensity of the point pattern by using the SPDE method from r-INLA.

Description

Estimates the intensity of the point pattern by using the SPDE method from r-INLA.

Usage

est_intens_spde(
  coords,
  win = NULL,
  npixel = 50,
  fine_mesh = FALSE,
  mesh = NULL,
  weights = NULL,
  alpha = 2,
  ...
)

Arguments

coords

ppp object or matrix with x and y coordinates of the observed bombs

win

observation window, either of class owin or a matrix with the x and y coordinates of the boundary, not neccessary if coords is a ppp object

npixel

number of pixel per dimension (see spatstat.options)

fine_mesh

logical, if FALSE a coarse mesh will be created, if TRUE a fine mesh will be created, only used if argument mesh is NULL

mesh

(optional) a predefined mesh for the spde model

weights

(optional) integration weights for the spde model, only used if argument mesh is NULL

alpha

(optional) alpha value for the spde model, only used if argument spde is NULL

...

additional arguments for the construction of the spde model (see INLA/inla.spde2.matern documentation)

Value

A list of

intensest

Pixel image with the estimated intensities of the random field.

mesh

The mesh.

Examples

## Not run: 
data(craterA)
est_spde <- est_intens_spde(coords=craterA)
image.plot(list(x=est_spde$intensest$xcol, y=est_spde$intensest$yrow, 
                z=log(t(est_spde$intensest$v))), main="estimated logarithmic intensity")
points(craterA)

## End(Not run)