A statistical method for the detection of geographic clustering

Kernel-based, smoothed estimates of spatial variables are useful in exploratory analyses because they yield a clear visual image of geographic variability in the underlying variable. In this paper I suggest an approach for assessing the significance of peaks in the surface that result from the application of the smoothing kernel. The approach may also be thought of as a method for assessing the maximum among a set of suitably defined local statistics. Local statistics for data on a regular grid of cells are first defined by using a Gaussian kernel. Results from integral geometry are then used to find the probability that the maximum local statistic (M) exceeds a given critical value (M). Approximations are provided that make implementation of the approach straightforward. Future work will address several other issues associated with local statistics that have been defined in this way, including edge effects, and the effects of global spatial autocorrelation on the choice of critical value.