SPATIALREG Procedure

Spatial Moving Average Models

The spatial moving average (SMA) model accounts for spatial dependence among the error terms; thus it is similar to the SEM model but with a different autocorrelation structure. The SMA model is used for modeling local autocorrelation. Let y Subscript i denote the observation associated with the spatial unit bold s Subscript i for i equals 1 comma 2 comma ellipsis comma n. For these spatial units, let bold upper W be an n times n spatial weights matrix. Further, let bold x Subscript i be a p times 1 vector that denotes values of p regressors recorded at unit bold s Subscript i.

The SMA model can be described in vector form by using the following two-stage formulation,

bold y equals bold upper X bold-italic beta plus bold u
bold u equals left-parenthesis bold upper I Subscript n Baseline minus lamda bold upper W right-parenthesis bold-italic epsilon

where bold y equals left-parenthesis y 1 comma y 2 comma ellipsis comma y Subscript n Baseline right-parenthesis prime and bold-italic epsilon equals left-parenthesis epsilon 1 comma epsilon 2 comma ellipsis comma epsilon Subscript n Baseline right-parenthesis prime, with epsilon Subscript i Baseline tilde Overscript normal i normal i normal d Endscripts upper N left-parenthesis 0 comma sigma squared right-parenthesis. Moreover, bold upper X is an n times p matrix that has bold x prime Subscript i in each row, and bold upper Z is an n times q matrix that has of bold z prime Subscript iin each row. In addition, bold-italic beta is a p times 1 parameter vector.

The log-likelihood function for the SMA model is

script upper L equals minus StartFraction n Over 2 EndFraction ln left-parenthesis 2 pi sigma squared right-parenthesis minus StartFraction left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X bold-italic beta right-parenthesis right-bracket prime left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X bold-italic beta right-parenthesis right-bracket Over 2 sigma squared EndFraction minus ln StartAbsoluteValue bold upper B EndAbsoluteValue

where bold upper B equals bold upper I Subscript n Baseline minus lamda bold upper W.

For the SMA model, the gradients are

StartFraction partial-differential script upper L Over partial-differential bold-italic beta EndFraction equals StartFraction left-parenthesis bold upper B Superscript negative 1 Baseline bold upper X right-parenthesis prime left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X bold-italic beta right-parenthesis right-bracket Over sigma squared EndFraction
StartFraction partial-differential script upper L Over partial-differential lamda EndFraction equals minus StartFraction 1 Over sigma squared EndFraction left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X bold-italic beta right-parenthesis right-bracket prime left-bracket bold upper B Superscript negative 1 Baseline bold upper W right-bracket left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X bold-italic beta right-parenthesis right-bracket plus normal t normal r left-parenthesis bold upper B Superscript negative 1 Baseline bold upper W right-parenthesis
StartFraction partial-differential script upper L Over partial-differential sigma squared EndFraction equals minus StartFraction n Over 2 sigma squared EndFraction plus StartFraction left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X bold-italic beta right-parenthesis right-bracket prime left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X bold-italic beta right-parenthesis right-bracket Over 2 sigma Superscript 4 Baseline EndFraction
Last updated: June 19, 2025