est_intens: Estimates the intensity of the point pattern.

In highriskzone: Determining and Evaluating High-Risk Zones

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)

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

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.

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)

highriskzone documentation built on Aug. 29, 2023, 5:10 p.m.