kde2dnew.fortran: A 2-d normal kernel estimator

In etasFLP: Mixed FLP and ML Estimation of ETAS Space-Time Point Processes for Earthquake Description

A 2-d normal kernel estimator

Description

A simple and quick 2-d weighted normal kernel estimator, with fixed bandwidth and relative integral.

Usage

kde2dnew.fortran(
#        parallel=FALSE,
         xkern,ykern,gx,gy,h,
         factor.xy=1,eps=0,w=replicate(length(xkern),1),
hvarx=replicate(length(xkern),1),hvary=replicate(length(xkern),1)
)

kde2d.integral(xkern,ykern,gx=xkern,gy=ykern,eps=0,factor.xy=1,
           h  = c(   bwd.nrd(xkern,w),bwd.nrd(ykern,w)),w=replicate(length(xkern),1),
           hvarx=replicate(length(xkern),1),hvary=replicate(length(xkern),1)
           )

Arguments

xkern	x-values of kernel points of length n (n=length(xkern)).
ykern	y-values of kernel points of length n.
gx	x-values of the points where densities must be estimated.
gy	y-values of the points where densities must be estimated.
h	bandwidths: a length 2 numerical vector.
eps	enlargment factor for the region of interest.
factor.xy	expansion factor for bandwidths (density will be smoother if factor.xy>1).
w	vector of weights to give to observed points (length n).

[]

hvarx	Longitude bandwidths adjustement used in the kernel estimator of background seismicity. The length must be equal to the number of events of the catalog after event selection (can be less than nrow(cat.orig)). Default value = replicate(length(xkern),1)
hvary	Longitude bandwidths adjustement used in the kernel estimator of background seismicity. The length must be equal to the number of events of the catalog after event selection (can be less than nrow(cat.orig)).Default value = replicate(length(xkern),1)

Details

A standard bivariate normal kernel estimator.

Value

grid values and estimated densities.

Author(s)

Marcello Chiodi.

References

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer. Wand, M.P and Jones, M.C. (1995). Kernel Smoothing. London: Chapman & Hall/CRC.

etasFLP documentation built on Sept. 14, 2023, 5:06 p.m.