CSPATIALREG Procedure

Spatial Durbin Moving Average Models

The term spatial Durbin moving average (SDMA) model is used to refer to the SMA model that accommodates exogenous interaction effects. Let y Subscript i denote the observation that is 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 that are recorded at unit bold s Subscript i, and let bold z Subscript i be a q times 1 vector that denotes values of q covariates that are measured at unit bold s Subscript i.

The SDMA model can be described in vector form as

where bold y equals left-parenthesis y 1 comma y 2 comma ellipsis comma y Subscript n Baseline right-parenthesis prime; 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; bold upper X and bold upper Z are n times p and n times q matrices, where each row consists of bold x prime Subscript i and bold z prime Subscript i, respectively; bold-italic beta is a p times 1 parameter vector; and bold-italic theta is a q times 1 parameter vector.

By letting upper X overTilde equals left-bracket bold upper X bold upper W bold upper Z right-bracket and bold-italic beta overTilde equals left-parenthesis bold-italic beta prime bold-italic theta Superscript prime Baseline right-parenthesis prime, the SDMA model can be written as

bold y equals bold upper X overTilde bold-italic beta overTilde plus bold upper B bold-italic epsilon

The log-likelihood function for the SDMA 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 overTilde bold-italic beta overTilde right-parenthesis right-bracket prime left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X overTilde bold-italic beta overTilde 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 and StartAbsoluteValue bold upper B EndAbsoluteValue denotes the determinant of matrix bold upper B.

For the SDMA model, the gradients are

StartFraction partial-differential script upper L Over partial-differential bold-italic beta overTilde EndFraction equals StartFraction left-parenthesis bold upper B Superscript negative 1 Baseline bold upper X overTilde right-parenthesis prime left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X overTilde bold-italic beta overTilde 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 overTilde bold-italic beta overTilde 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 overTilde bold-italic beta overTilde 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 overTilde bold-italic beta overTilde right-parenthesis right-bracket prime left-bracket bold upper B Superscript negative 1 Baseline left-parenthesis bold y minus bold upper X overTilde bold-italic beta overTilde right-parenthesis right-bracket Over 2 sigma Superscript 4 Baseline EndFraction
Last updated: July 09, 2026