CSPATIALREG Procedure

Spatial Durbin Models

Unlike a SAR model, a spatial Durbin model (SDM) can account for exogenous interaction effects in addition to the endogenous interaction effects. Let y Subscript i denote the observation that is associated with a 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 of your choice. Further assume that bold x Subscript i is a p times 1 vector that denotes values of p regressors that are recorded for the spatial unit bold s Subscript i and assume that bold z Subscript i is a q times 1 vector that denotes values of q regressors that are measured at unit bold s Subscript i.

The SDM can be described in vector form as (LeSage and Pace 2009)

bold y equals rho bold upper W bold y plus bold upper X bold-italic beta plus bold upper W bold upper Z bold-italic theta 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, 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 is an n times p matrix where each row consists of bold x prime Subscript i, bold upper Z is an n times q matrix where each row consists of bold z prime Subscript i, and bold-italic beta and bold-italic theta are p times 1 and q times 1 parameter vectors, respectively.

By letting bold 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, you can rewrite the SDM as

bold y equals rho bold upper W bold y plus bold upper X overTilde bold-italic beta overTilde plus bold-italic epsilon

The log-likelihood function for the SDM is

script upper L equals minus StartFraction n Over 2 EndFraction ln left-parenthesis 2 pi sigma squared right-parenthesis minus StartFraction left-parenthesis bold upper A bold y minus bold upper X overTilde bold-italic beta overTilde right-parenthesis prime left-parenthesis bold upper A bold y minus bold upper X overTilde bold-italic beta overTilde right-parenthesis Over 2 sigma squared EndFraction plus ln StartAbsoluteValue bold upper A EndAbsoluteValue

where bold upper A equals bold upper I Subscript n Baseline minus rho bold upper W.

For the SDM, the gradients are

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

Both the SDM and the SAR model account for endogenous interaction effects. However, in some cases there might be an interaction among error terms. In such cases, you might consider a spatial error model, which addresses spatial interaction among error terms.

Last updated: July 09, 2026