highriskzone (version 1.4.9)

est_intens: Estimates the intensity of the point pattern.

Description

Estimates the intensity of the point pattern by a kernel method (See density.ppp).

Usage

est_intens(ppdata, covmatrix = NULL, weights = NULL)

Value

A list of

intensest

Estimated intensity (object of class "im", see density.ppp).

covmatrix

Covariance matrix. If covmatrix = NULL, the matrix is estimated by Hscv.

Arguments

ppdata

data of class ppp

covmatrix

(Optional) Covariance matrix of the kernel of a normal distribution

weights

(Optional) vector of weights attached to each observation

See Also

density.ppp, Hscv, eval.im

Examples

data(craterA)
#change npixel = 50 to 1000 to get a nicer picture
spatstat.geom::spatstat.options(npixel=50)
# use only ten observations for fast computation
thin.craterA <- craterA[1:10]
int <- est_intens(thin.craterA)
# Plot estimated intensity
plot(int$intensest, main = "pixel image of intensity")
plot(craterA$window, main = "contour plot of intensity")
contour(int$intensest, add =TRUE)