Calculates an estimate of the K function (M.A. Mugglestone, S.A. Harding, B.Y.Y. Lee, P.J. Diggle & B.S. Rowlingson).
|What to print (
||Vertical coordinates of each spatial point pattern; no default – this parameter must be set|
||Horizontal coordinates of each spatial point pattern; no default – this parameter must be set|
||Vertical coordinates of each polygon; no default – this parameter must be set|
||Horizontal coordinates of each polygon; no default – this parameter must be set|
||Vectors of distances to use; no default – this parameter must be set|
||Variates to receive the estimated values of the K function|
The K function, or reduced second-order moment function, relates to the distribution of the inter-event distances between all ordered pairs of events in a spatial point pattern (see Diggle 1983). The function is formally defined as the expected number of further events within distance s of an arbitrary event, divided by the overall density of events per unit area. An approximately unbiased estimator of K which incorporates corrections for edge effects can be obtained using the method of Ripley (1976). This estimator, denoted by K^(s), is essentially an empirical distribution function (EDF) of weighted inter-event distances. The weight associated with an inter-event distance s derived from events at positions (x1, y1) and (x2, y2) is the reciprocal of the conditional probability that any event separated from the point (x1, y1) by the distance s will fall in the study region and so be observed.
The term complete spatial randomness (CSR) is used to represent the hypothesis that the overall density of events in a spatial point pattern is constant throughout the study region, and that the events are distributed independently and uniformly. Under CSR, the expected number of further events which lie within a distance s of an arbitrary event in the pattern is simply the area of a circle of radius s, multiplied by the overall density of events. Thus, the K function for a completely random pattern is given by
K(s) = π × s2 .
The K function for a clustered (regular) pattern will tend to be larger (smaller) than the values given by the above expression, at least for small distances.
The procedure KHAT calculates Ripley’s (1976) estimator for K given the coordinates of a spatial point pattern (specified by the parameters
Y), the coordinates of a polygon containing the points (specified by the parameters
YPOLYGON) and a vector of distances at which to calculate the EDF of K (specified by the parameter
S). The output of the procedure is a vector of estimates of K corresponding to the distances in
S. The estimated K function can be saved using the parameter
Printed output is controlled using the
summary prints the distances at which the K function is estimated and the estimates themselves under the headings
PTCHECKXY is called to check that
Y have identical restrictions. A similar check is made on
YPOLYGON. The procedure then calls
PTCLOSEPOLYGON to close the polygon specified by
SORT function is then used to create a variate containing the distances in
S arranged in ascending order. (The original variate is left unchanged.) The procedure then calls a procedure
PTPASS to call a Fortran program to calculate an edge-corrected estimate of the K function.
S may be restricted, as long as
X has the same restriction as
XPOLYGON has the same restriction as
YPOLYGON. Only the subset of values specified by each restriction will be included in the calculations.
Diggle, P.J. (1983). Statistical Analysis of Spatial Point Patterns. Academic Press, London.
Ripley, B.D. (1976). The second-order analysis of stationary point processes. Journal of Applied Probability, 13, 255-266.
Commands for: Spatial statistics.
CAPTION 'KHAT example'; STYLE=meta VARIATE pinex,piney READ [SETNVALUES=yes] pinex,piney 0.09 0.91 0.02 0.71 0.03 0.62 0.18 0.61 0.03 0.52 0.02 0.41 0.16 0.35 0.13 0.33 0.13 0.27 0.03 0.21 0.13 0.14 0.08 0.11 0.02 0.02 0.18 0.98 0.31 0.89 0.22 0.58 0.13 0.52 0.21 0.38 0.23 0.27 0.23 0.11 0.41 0.98 0.44 0.97 0.42 0.93 0.42 0.48 0.43 0.36 0.59 0.92 0.63 0.92 0.63 0.88 0.66 0.88 0.58 0.83 0.53 0.69 0.52 0.68 0.49 0.58 0.52 0.48 0.52 0.09 0.58 0.06 0.68 0.66 0.68 0.63 0.67 0.53 0.67 0.48 0.67 0.41 0.68 0.34 0.66 0.24 0.73 0.27 0.74 0.11 0.78 0.06 0.79 0.02 0.86 0.03 0.84 0.88 0.94 0.89 0.95 0.83 0.79 0.79 0.84 0.71 0.83 0.68 0.86 0.65 0.79 0.61 0.93 0.48 0.83 0.42 0.93 0.31 0.93 0.23 0.97 0.64 0.96 0.64 0.96 0.61 0.96 0.57 0.97 0.38 : VARIATE xpoly; VALUES=!(0,1,1,0) & ypoly; VALUES=!(0,0,1,1) & s; VALUES=!(0.01,0.02...0.1) KHAT [PRINT=*] Y=piney; X=pinex; YPOLYGON=ypoly; XPOLYGON=xpoly;\ S=s; KVALUES=kpines CALCULATE kcsr = CONSTANTS('pi')*(s**2) PRINT s,kpines,kcsr