---
title: Local Moran’s I (LISA)
url: https://www.emergentmind.com/topics/local-moran-s-i-lisa
type: topic
---

# Local Moran’s I (LISA)

Local Moran’s I (LISA) provides a formal, spatially explicit decomposition of spatial autocorrelation by quantifying the association between a value at a geographic location and the values observed at its spatial neighbors. As a local indicator of spatial association, it extends the global Moran’s I statistic, producing a separate measure for each spatial unit (“location”) that can detect local clusters (“hot-spots” and “cold-spots”) as well as local spatial outliers. The rigorous mathematical foundation, multiple forms of weights normalization, inferential permutation testing, and recent advances in visualization and normalization schemes underpin its broad uptake in geospatial modeling, spatial econometrics, epidemiology, and other fields.

## 1. Mathematical Definitions and Normalization Schemes

The classical (univariate) Local Moran’s I for an attribute vector $x = (x_1,\ldots,x_n)$ over $n$ spatial units is defined, after mean-removal and (optionally) standardization, using a spatial weights matrix $W = (w_{ij})$ to encode the neighborhood structure. Mathematically:

\[
I_i = \sum_{j=1}^n w_{ij}\,(x_i - \bar{x})\,(x_j - \bar{x})
\]
where $\bar{x}$ is the global mean. This is formally a local “gamma index” of matrix similarity, and up to proportionality can be seen as the $i$th row "projection" of the global Moran’s I calculation [2110.10622].

Normalization of both weights and variables fundamentally alters the inferential and interpretive properties. Three prominent variants have been formally analyzed [2202.11207]:

| Index | Weights Matrix | Variable Type | Formula | Proportionality to Global $I$? |
|-------|---------------|--------------|---------|-------------------------------|
| MI₁   | un-normalized $v_{ij}$   | mean-centered $y_i$  | $I_i^{(1)} = \frac{\sum_j v_{ij} y_i y_j}{o^2}$ | Yes (factor $S_0$) |
| MI₂   | row-normalized $w_{ij}^*$| standardized $z_i$   | $I_i^{(2)} = z_i \sum_j w_{ij}^* z_j$ | No |
| MI₃   | global-normalized $w_{ij}$ | standardized $z_i$  | $I_i^{(3)} = z_i \sum_j w_{ij} z_j$   | Yes ($=I$) |

Here, $y_i = x_i - \bar{x}$, $z_i = y_i/o$ (where $o^2$ is the population variance), $v_{ij}$ is the raw (un-normalized) adjacency, $w_{ij}^* = v_{ij}/S_i$ with $S_i = \sum_j v_{ij}$, and $w_{ij} = v_{ij}/S_0$ with $S_0 = \sum_{ij} v_{ij}$.

Of these, MI₃ (using global-normalized weights and standardization by population variance) uniquely satisfies Anselin’s second requirement that “the sum of LISAs for all observations is proportional to a global indicator of spatial association,” with constant $\gamma=1$ [2202.11207].

## 2. Panel and Multivariate Extensions

For spatial panel data—multiple outcomes or residuals $y_{t,i}$ measured for each region $i$ across time points $t=1,\ldots,T$—Local Moran’s I generalizes to a multivariate inner product form. Let $y_i \in \mathbb{R}^T$ denote the $T$-dimensional vector for region $i$, $\bar{y}$ the across-region mean vector, and $M$ a positive-definite $T\times T$ matrix (often the identity). The panel-data LISA is:

\[
I_i = \sum_{j=1}^n w_{ij}\, \langle y_i - \bar{y},\ y_j - \bar{y} \rangle_M
= \sum_{j=1}^n w_{ij}\, (y_i - \bar{y})^\mathsf{T} M (y_j - \bar{y})
\]

Matrix notation enables the simultaneous calculation of all $I_i$ entries as the diagonal of $[(Y - \bar{Y})^\mathsf{T} M (Y - \bar{Y}) W]$, where $Y$ is the $T \times n$ data matrix [2110.10622].

A plausible implication is that this generalization retains theoretical properties (e.g., local/global decomposability under global-normalized weights) under appropriate centering and normalization, but with the local statistic now quantifying multivariate (“vector”) similarity.

## 3. Computational Algorithms and Inference

LISA values are not interpreted without significance testing. A standard workflow is:

1. Standardize variables (usually via z-scores or mean-centering and population variance scaling).
2. Construct an appropriate weights matrix (e.g., binary adjacency, inverse distance, with possible normalizations).
3. For each region $i$, compute $I_i$.
4. Assess significance via permutation: with region $i$’s value fixed, permute neighbor labels $M$ times, recompute $I_i$ for each, and estimate a pseudo $p$-value
\[
\text{pseudo-}p_i = \frac{\#\{|I^{(k)}_i| \ge |I_i|\} + 1}{M+1}
\]
5. Apply multiple-testing correction, e.g., Benjamini–Hochberg FDR, to control for family-wise error [2408.02418, 2110.10622].

For spatial panel models, analytic (non-sampling) bounds based on concentration inequalities can provide fast $O(\text{number of edges})$ computation of $p$-values, avoiding the cost of $(n-1)!$ enumerations or Monte Carlo approximation:
\[
\Pr\{\gamma_i(\pi)-m_i\,\bar{\lambda}_{-i} \ge \gamma_i\} \le \frac{1}{\sqrt{\pi} \Gamma(\cdots; 1/2) + O(n^{-4})}
\]
where $\Gamma$ is the upper incomplete gamma, and the constants are calculated from the weight structure and similarity function [2110.10622].

## 4. Interpretation and Classification of Local Patterns

The sign and magnitude of each $I_i$ are interpreted in conjunction with the value at $i$ and the spatial lag $\text{lag}_i = \sum_j w_{ij} z_j$. Classical quadrant-based classification yields four cluster/outlier types:

- **High–High (HH)**: $I_i > 0$, $z_i > 0$, $\text{lag}_i > 0$—localized hot-spot (high value surrounded by highs)
- **Low–Low (LL)**: $I_i > 0$, $z_i < 0$, $\text{lag}_i < 0$—cold-spot (low value among lows)
- **High–Low (HL)**: $I_i < 0$, $z_i > 0$, $\text{lag}_i < 0$—high-value outlier amidst lows
- **Low–High (LH)**: $I_i < 0$, $z_i < 0$, $\text{lag}_i > 0$—low outlier among highs

A significant positive $I_i$ indicates the attribute at $i$ moves with its neighbors; a significant negative $I_i$ marks a spatial outlier. Final map labeling should always consider both the sign of $I_i$ and the accompanying $z_i$/$\text{lag}_i$ values [2408.02418].

## 5. Visualization Techniques

Advanced interactive visualizations have been proposed to facilitate exploration and interpretation of Local Moran’s I statistics. Three systems highlighted include [2408.02418]:

- **Moran Dual-Density Plot**: Shows the distribution of $z$-scores for the whole study region (with neighbors marked), the spatial lag, the computed $I_i$, and the permutation-based significance.
- **Moran Network Scatterplot**: Displays the Moran scatterplot ($z_i$ vs. $\text{lag}_i$), with neighbor relations rendered interactively (edges colored and scaled by weight, neighbor cluster labels).
- **Spatial Lag Radial Plot**: Radially arranges neighbors around $i$ by bearing, with distance set by $z_j$, and marker size by $w_{ij}$; a dashed circle indicates $\text{lag}_i$.

Brushing and linking allow coordinated updates across views, with map-based, scatterplot, and spatial context jointly updating as the user examines different locations, promoting holistic understanding of spatial patterns.

## 6. Comparative Properties and Best Practices

Key methodological considerations and comparisons include:

- Global-normalized weights and population variance standardization (MI₃) uniquely ensure that LISA values sum exactly to the global Moran’s I [2202.11207].
- Row normalization is often used in practical software and visualization for interpretability, but such MI₂ LISA values do not decompose the global measure.
- Choice of weights (contiguity, distance, or inverse distance) and normalization method critically impact detected patterns.
- Multiple testing is inherent: reporting cluster significance demands Type I error control; FDR procedures are standard [2110.10622, 2408.02418].
- Caution is warranted for sparse neighborhoods: unstable significance inference and wide null distributions can result for locations with very few neighbors.
- For panel models, LISA can be directly applied to model residuals, supporting model checking for spatial dependence in spatial econometric regression [2110.10622].

## 7. Practical Applications and Case Studies

Local Moran’s I (LISA) is applied to diverse spatial problems. Notable recent implementations include:

- **Electoral Panel Model**: Over 3,100 U.S. counties across five presidential elections (2000–2016). Covariates included population density, median income, and racial composition. LISA identified statistically significant Republican “corridors” and Democratic clusters, guiding substantive political analysis [2110.10622].
- **Cancer Mortality Studies**: Age-adjusted cancer deaths in U.S. counties (2011–2020), with HH and LL clusters mapped at fine spatial scale. Visualization dashboards (available as open-source moranplot library) bring transparency and contextualization to spatial “hot-spot” findings [2408.02418].
- **Urban Demography**: Empirical verification on Beijing–Tianjin–Hebei regional data demonstrates the exactness of MI₃ sum-to-global proportionality and illustrates how alternate weightings affect the cluster identification [2202.11207].

These studies exemplify the broad applicability and rigorous interpretability afforded by formally constructed LISA statistics, in both univariate and panel settings.

Source: https://www.emergentmind.com/topics/local-moran-s-i-lisa