Join Count Statistics

Join count statistics are a method of spatial analysis used to assess the degree of association, in particular the autocorrelation, of categorical variables distributed over a spatial map. They were originally introduced by Australian statistician P. A. P. Moran. Join count statistics have found widespread use in econometrics, remote sensing and ecology.Dale MR, Fortin MJ. Spatial analysis: a guide for ecologists. Cambridge University Press; 2014 Sep 11. Join count statistics can be computed in a number of software packages including PASSaGE, GeoDA, PySAL and spdep.

Binary data

Given binary data $x_i \in \$ distributed over $N$ spatial sites, where the neighbour relations between regions $i$ and $j$ are encoded in the spatial weight matrix
:$w_ = \begin 1 \qquad &i\textj\\ 0 &\text \end$
the join count statistics are defined as

:$J = J_ + J_ + J_$
Where
:$J_ = \frac\sum_ w_ x_i x_j$
:$J_ = \frac\sum_ w_ (x_i-x_j)^2$
:$J_ = \frac\sum_ w_ (1-x_i) (1-x_j)$
:$J = \frac\sum_ w_$
The $B,W$ subscripts refer to 'black'=1 and 'white'=0 sites. The relation $J = J_ + J_ + J_$ implies only three of the four numbers are independent. Generally speaking, large values of $J_$ and $J_$ relative to $J_$ imply autocorrelation and relatively large values of $J_$ imply anti-correlation.

To assess the statistical significance of these statistics, the expectation under various null models has been computed.Sokal RR, Oden NL. Spatial autocorrelation in biology: 1. Methodology. Biological journal of the Linnean Society. 1978 Jun 1;10(2):199-228. For example, if the null hypothesis is that each sample is chosen at random according to a Bernoulli process with probability
:$p = \frac = \frac$
then Cliff and Ord show that
:$E(J_) = \frac S_0 p^2$
:
var(J_) = \frac ( S_1(1-p) + S_2p

:$E(J_) = S_0 p(1-p)$
:
var(J_) = \frac 4 S_1 + S_2(1-4p(1-p))

where
:$S_0 = \sum_ w_$
:$S_1 = \frac\sum_(w_ + w_)^2$
:$S_2 = \sum_( \sum_j w_ + \sum_j w_)^2$
However in practice an approach based on random permutations is preferred, since it requires fewer assumptions.

Local join count statistic

Anselin and Li introducedAnselin L, Li X. Operational local join count statistics for cluster detection. Journal of geographical systems. 2019 Jun 1;21:189-210. the idea of the local join count statistic, following Anselin's general idea of a Local Indicator of Spatial Association (LISA). Local Join Count is defined by e.g.
:$J_ = x_i \sum_j w_ x_j$
with similar definitions for $BW$ and $WW$. This is equivalent to the Getis–Ord statistics computed with binary data. Some analytic results for the expectation of the local statistics are available based on the hypergeometric distribution but due to the multiple comparisons problem a permutation based approach is again preferred in practice.

Extension to multiple categories

When there are $k \geq 2$ categories join count statistics have been generalised
:$J_ = \frac \sum_ I_r(x_i) I_s(x_j)$
Where $I_r(x_i) = \delta_$ is an indicator function for the variable $x_i$ belonging to the category $r$. Analytic results are availableEpperson, B.K., 2003. Covariances among join-count spatial autocorrelation measures. Theoretical Population Biology, 64(1), pp.81-87. or a permutation approach can be used to test for significance as in the binary case.

Spatial analysis
Covariance and correlation

References

{{reflist