---
title: Getis-Ord Gi* Statistic
url: https://www.emergentmind.com/topics/getis-ord-gi-statistic
type: topic
---

# Getis-Ord Gi* Statistic

The Getis-Ord Gi* statistic is a local indicator of spatial association (LISA) specifically designed to quantify the presence and intensity of high-value (“hotspot”) or low-value (“coldspot”) spatial clusters within a spatial dataset. Originally introduced in the spatial statistics and spatial econometrics literature, Gi* is widely used across fields such as geostatistics, epidemiology, environmental science, crime analysis, and remote sensing to detect spatial heterogeneity, evaluate spatial association, and map spatially significant anomaly regions at granular scales. Rigorous developments and recent methodological advances have strengthened its interpretation, computational efficiency, and theoretical connection to global autocorrelation indices such as Moran’s I.

## 1. Mathematical Basis and Definition

The Getis-Ord Gi* statistic evaluates, for each spatial unit $i$, whether the attribute values in its local neighborhood are significantly higher or lower than would be expected under global spatial randomness. The classical form of Gi*, as standardized to a $z$-score, is:

\[
G_i^* = \frac{ \sum_{j=1}^n w_{ij}x_j - \bar x \sum_{j=1}^n w_{ij} }{ S \sqrt{ \frac{n\sum_{j=1}^n w_{ij}^2 - \left( \sum_{j=1}^n w_{ij} \right)^2 }{n-1} } }
\]

where:
- $x_j$ is the observed value at location $j$ (e.g., count/rate/intensity).
- $w_{ij}$ is the spatial weight encoding the strength of connection between $i$ and $j$.
- $n$ is the total number of spatial units.
- $\bar x$ is the global mean.
- $S$ is the global (sample) standard deviation.

This standardized statistic allows direct use of the standard normal distribution for inference under the null hypothesis of spatial randomness. High positive $G_i^*$ indicates a hotspot, while large negative $G_i^*$ marks a coldspot [2508.12188, 2506.03356, 2012.08647].

Both the “self” and “neighbor” contributions are included, distinguishing Gi* from its exclusionary counterpart $G_i$, though both can be rewritten in quadratic or matrix notation for large-scale or eigendecomposition-based analyses [2508.19908, 1803.06139].

## 2. Spatial Weights and Neighborhood Construction

The choice and construction of the spatial weights matrix $W=[w_{ij}]$ are central to Gi*’s local character. Different study domains employ different schemes:
- **Binary contiguity:** $w_{ij}=1$ if $i$ and $j$ are contiguous (e.g., Queen’s case—sharing a side or vertex on a political map/grid); zero otherwise [2508.12188, 2506.03356].
- **Row-standardization:** Each row is normalized such that $\sum_j w_{ij}=1$, equalizing influence from varying neighbor counts [2508.12188].
- **Distance bands:** $w_{ij}=1$ if the distance $d_{ij}\leq D$, for some fixed threshold; $w_{ij}=0$ otherwise [2601.12548].
- **Continuous kernels:** $w_{ij}$ decay inversely or exponentially with $d_{ij}$, relevant for gravity-model or interaction potential reformulations [2508.19908, 1803.06139].

All these choices affect sensitivity: larger or more numerous neighbors (wider spatial reach) dilute local extreme values; tight, sparse weighting accentuates highly localized clusters. Specialized domains (e.g., fine-grained urban grids, state-level epidemiological analyses) require weight construction consistent with the scale and the underlying spatial process [2506.03356, 2508.12188].

## 3. Computation, Standardization, and Inference

### Computation Workflow

The procedure consists of:
1. Selection or calculation of local variable(s) $x_j$ (counts, rates, attributes), often standardized or normalized in advance (or, for rates, directly used).
2. Construction of spatial weights $w_{ij}$ via the chosen adjacency, contiguity, or distance-based scheme.
3. For each unit $i$:
   a. Compute the weighted sum in the numerator.
   b. Calculate the mean and standard deviation over all $x_j$.
   c. Standardize using the denominator to obtain $G_i^*$.
4. Assess significance by comparing $G_i^*$ against critical $z$-values for the desired confidence level.

Pseudocode implementations and full variable flow have been detailed for spatial point pattern analysis [2601.12548], raster/grid analysis [2506.03356], and even integration within convolutional neural network layers [1912.10667].

### Inference and Multiple Testing

The classic hypothesis testing approach is analytic, interpreting $G_i^*$ as a $z$-score under approximate normality. Hotspots correspond to $G_i^*$ above a chosen percentile, often $|G_i^*|>1.96$ for a $p<0.05$ two-sided test [2508.12188, 2601.12548]. Alternative inferential strategies include:
- **Permutation-based inference:** Randomly reassign $x_j$ among locations to generate the null distribution of $G_i^*$ [2012.08647].
- **Computation-free nonparametric testing:** Closed-form analytic bounds for $p$-values using Khintchine-type inequalities, requiring only $O(n)$ complexity and suitable even for non-Gaussian, small-$n$ scenarios [2012.08647].

Formal multiple-testing correction is rarely applied in published work; thus, some detected clusters may be spurious in analyses reporting results for many regions or features [2508.12188, 2506.03356].

## 4. Theoretical Framework and Relationship to Global Indices

Recent mathematical work situates Gi* within a unified spatial statistics framework:
- **Quadratic form and potential theory:** The global Getis-Ord index is a quadratic form $G=\mathbf{p}^T W \mathbf{p}$, and local Gi* is $G_i=(W\mathbf{p})_i$, where $\mathbf{p}$ is the normalized attribute vector [2508.19908, 1803.06139].
- **Gravity-model equivalence:** With suitably defined kernels $v_{ij}=f(r_{ij})$, local Gi* is directly proportional to classical gravity-model potential at $i$. Hence, Gi* can be interpreted as measuring local “spatial interaction intensity” [2508.19908, 1803.06139].
- **Decomposition of Moran’s I:** Moran’s I can be decomposed into the global Getis-Ord index, the sum of local Gi*, a size-correlation function, and the number of elements, reflecting both global structure and local clustering [2508.19908].
- **Scatterplot diagnostics:** Gi* values versus unitized variables can be visualized in scatterplots analogous to Moran’s I plots, partitioning observations into high-high, high-low, low-high, and low-low quadrants [1803.06139].

This framework clarifies that local clustering (hotspots/coldspots) is not merely a secondary effect, but a constituent part of global spatial autocorrelation structure.

## 5. Applications Across Spatial Domains

Gi* is applicable in any setting with indexed spatial data and an interest in mapping local association. Notable domains and implementations include:

| Domain                   | Attribute/Unit               | Spatial Weight           |
|--------------------------|------------------------------|-------------------------|
| Chronic disease mapping  | State-level mortality rates  | Queen’s contiguity (row-standardized) [2508.12188]           |
| Urban accident analysis  | Grid cell crash counts       | Queen contiguity (binary) [2506.03356]                      |
| Traffic collision severity | Pointwise collision events (weights = severity) | Fixed distance band (binary) [2601.12548]                |
| Electoral spatial analysis| Vote share per precinct     | Adjacency graph (binary) [2012.08647]                      |
| Urban systems science    | City populations             | Distance-decay kernels [2508.19908, 1803.06139]             |
| Deep learning for remote sensing | CNN activations (feature maps) | Local windows with distance-based weights [1912.10667]      |

The Gi* statistic’s flexibility allows adjustment for varying data resolutions, attribute types (counts, rates, intensities), and neighbor definitions, facilitating both classical (e.g., epidemiological, urban planning) and novel (e.g., spatial pooling in neural nets) applications.

## 6. Interpretation, Visualization, and Caveats

Interpretation is standardized: significant, large positive Gi* indicates a spatial “hotspot” (local cluster of high values), while large negative Gi* marks a “coldspot” (local low-value cluster). Visual mappings typically employ warm/cool color scales to represent hotspots/coldspots, neutral for non-significant areas [2508.12188, 2506.03356].

Caveats include:
- **Scale sensitivity and MAUP:** The result depends on the spatial resolution and weight construction; changing grid size or adjacency alters hotspot detection [2506.03356].
- **Aggregation bias:** Coarse units may mask within-unit heterogeneity; observed hotspots may be artefacts of zonation.
- **Null distribution validity:** Analytic normality is approximate, especially for skewed or heavy-tailed $x_j$ or small $n$; permutation or computation-free bounds may be preferable [2012.08647].
- **Multiple testing:** Numerous parallel Gi* tests inflate type-I error rate; the absence of systematic correction is a limitation [2508.12188].
- **Purely descriptive nature:** Gi* detects spatial association but does not infer causality or adjust for confounding/covariate effects.

Mapping conventions, critical value selection, and cluster frequency tallies provide operational guidance for result reporting and policy-oriented spatial targeting [2508.12188, 2601.12548].

## 7. Extensions and Computational Developments

Extensions to the Gi* statistic have appeared in domains with complex spatial dependencies or high computational demands:
- **Computation-free nonparametric inference:** Analytic bounds substituting for Monte-Carlo permutation preserve statistical validity at vastly lower computational burden, enabling multi-scale or network-scale testing [2012.08647].
- **Integration in machine learning:** The Gi* statistic as a pooling mechanism in CNNs improves generalization for geospatial segmentation by enforcing spatially meaningful activation selection [1912.10667].
- **Structural decomposition:** Gi* serves as a building block for decomposing and interpreting global spatial autocorrelation indices (e.g., explicit mathematical linkages with Moran’s I), unifying local and global perspectives [2508.19908, 1803.06139].

The methodological versatility and deep theoretical foundation of Getis-Ord Gi* ensure its continuing prominence in spatial data science, spatial epidemiology, urban analytics, and spatially-informed machine learning.

Source: https://www.emergentmind.com/topics/getis-ord-gi-statistic