est_intens_spde: Estimates the intensity of the point pattern by using the...

In highriskzone: Determining and Evaluating High-Risk Zones

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)

highriskzone documentation built on Aug. 29, 2023, 3:01 p.m.