est_intensity_np: Estimate intensity surface by using non-parametric method.

In wangyuchen/sppmix: 2D Normal Mixture Models for Spatial Point Process

Description

Using Epanechnikov kernel to estimate intensity surface.

Usage

est_intensity_np(pattern, win = pattern$window, h, L = 10,
  kernel = c("Epanechnikov"), edgecorrect = TRUE, truncate = TRUE)

Arguments

pattern	A two-dimesional spatial point pattern with class ppp.
win	Object of class owin.
h	Kernel bandwith. h should be a postive number.
L	Number of grid for x and y axis.
kernel	Kernel used to estimate intensity surface. Currently, only support Epanechnikov kernel.
edgecorrect	Logical flag indicating whether to use edge-correction in estimating intensity surface. The default is TRUE.
truncate	Logical flag indicating whether truncation is used for pattern. The default is TRUE.

Value

An object of class im.

Examples

pattern1 <- rsppmix(demo_mix, lambda = 30, win = square(5))

# estimate and plot the estimated intensity surface
surf1 <- est_intensity_np(pattern1, h=0.05, L=100)
plot(surf1)

wangyuchen/sppmix documentation built on May 4, 2019, 12:58 a.m.