SPATIALREG Procedure

Spatial Autoregressive Confused Models

The spatial autoregressive confused (SAC) model, like the SARMA model, can accommodate spatial dependence in both the dependent variable and error terms. However, the covariance structure for the error terms in a SAC model is different from that of the SARMA model. 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 1 and bold upper W 2 be two spatial weights matrices. 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 SAC model can be described in vector form by using the following two-stage formulation (LeSage and Pace 2009),

bold y equals rho bold upper W 1 bold y plus bold upper X bold-italic beta plus bold u
bold u equals lamda bold upper W 2 bold-italic u plus 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 prime. Moreover, bold upper X is an n times p matrix that has bold x prime Subscript i in each row. In addition, bold-italic beta is a p times 1 parameter vector.

The log-likelihood function for the SAC 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 left-parenthesis bold upper A bold y minus bold upper X bold-italic beta right-parenthesis right-bracket prime left-bracket bold upper B left-parenthesis bold upper A bold y minus bold upper X bold-italic beta right-parenthesis right-bracket Over 2 sigma squared EndFraction plus ln StartAbsoluteValue bold upper A EndAbsoluteValue plus ln StartAbsoluteValue bold upper B EndAbsoluteValue

where bold upper A equals bold upper I Subscript n Baseline minus rho bold upper W 1 and bold upper B equals bold upper I Subscript n Baseline minus lamda bold upper W 2.

For the SAC model, the gradients are

StartFraction partial-differential script upper L Over partial-differential bold-italic beta EndFraction equals StartFraction left-parenthesis bold upper B bold upper X right-parenthesis prime left-bracket bold upper B left-parenthesis bold upper A 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 rho EndFraction equals StartFraction left-parenthesis bold upper B bold upper W 1 bold y right-parenthesis prime bold upper B left-parenthesis bold upper A bold y minus bold upper X bold-italic beta right-parenthesis Over sigma squared EndFraction minus normal t normal r left-parenthesis bold upper A Superscript negative 1 Baseline bold upper W 1 right-parenthesis
StartFraction partial-differential script upper L Over partial-differential lamda EndFraction equals StartFraction 1 Over sigma squared EndFraction left-bracket bold upper W 2 left-parenthesis bold upper A bold y minus bold upper X bold-italic beta right-parenthesis right-bracket prime left-bracket bold upper B left-parenthesis bold upper A bold y minus bold upper X bold-italic beta right-parenthesis right-bracket minus normal t normal r left-parenthesis bold upper B Superscript negative 1 Baseline bold upper W 2 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 left-parenthesis bold upper A bold y minus bold upper X bold-italic beta right-parenthesis right-bracket prime left-bracket bold upper B left-parenthesis bold upper A 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