SPATIALREG Procedure

Parameter Space of Spatial Coefficients

For all models except linear regression models in PROC SPATIALREG, the autoregressive parameters rho and lamda are often assumed to satisfy some assumptions to ensure consistency of the maximum likelihood estimator (Elhorst 2013). For SAR and SDM models, the Jacobian term involves the log-determinant of bold upper I minus rho bold upper W 1, and the parameter space of rho is often specified such that bold upper I minus rho bold upper W 1 is nonsingular. For SEM, SMA, SDM, and SDMA models, the Jacobian term involves the log-determinant of bold upper I minus lamda bold upper W 1, and the parameter space of lamda is often specified such that bold upper I minus lamda bold upper W 1 is nonsingular. For SAC, SDAC, SARMA, and SDARMA models, the Jacobian term involves the log-determinants of both bold upper I minus rho bold upper W 1 and bold upper I minus lamda bold upper W 1. As a result, the parameter space of rho and lamda is often specified such that both bold upper I minus rho bold upper W 1 and bold upper I minus lamda bold upper W 1 are nonsingular.

In the SPATIALREG procedure, the parameter space of spatial coefficients rho and lamda depends on the spatial weights matrix bold upper W that you choose. For bold upper W, the parameter space of the autoregressive parameters rho and lamda in PROC SPATIALREG is determined as follows:

  1. For a symmetric bold upper W, the nonsingularity condition requires rho element-of left-parenthesis omega Subscript min Superscript negative 1 Baseline comma omega Subscript max Superscript negative 1 Baseline right-parenthesis and lamda element-of left-parenthesis omega Subscript min Superscript negative 1 Baseline comma omega Subscript max Superscript negative 1 Baseline right-parenthesis. Here omega Subscript min and omega Subscript max denote the smallest (that is, most negative) and largest real eigenvalues of bold upper W, respectively.

  2. If bold upper W is symmetric and subsequently row-standardized, the nonsingularity condition requires rho element-of left-parenthesis omega Subscript min Superscript asterisk Baseline Superscript negative 1 Baseline comma 1 right-parenthesis and lamda element-of left-parenthesis omega Subscript min Superscript asterisk Baseline Superscript negative 1 Baseline comma 1 right-parenthesis. Here omega Subscript min Superscript asterisk denotes the smallest purely real eigenvalue of the row-standardized bold upper W.

  3. If bold upper W is asymmetric and subsequently row-standardized, the nonsingularity condition requires rho element-of left-parenthesis r Subscript min Superscript asterisk Baseline Superscript negative 1 Baseline comma 1 right-parenthesis and lamda element-of left-parenthesis r Subscript min Superscript asterisk Baseline Superscript negative 1 Baseline comma 1 right-parenthesis. Here r Subscript min Superscript asterisk denotes the smallest purely real eigenvalue of the row-standardized bold upper W.

  4. When Taylor approximation or Chebyshev approximation is used for SAR and SDM models, bold upper W is required to be row-standardized. In these cases, the restriction on the autoregressive coefficient rho is rho element-of left-parenthesis negative 1 comma 1 right-parenthesis.

    For the interpretation of the spatial coefficient rho, a significant and positive (or negative) value of rho in the SAR, SDM, SAC, SDAC, SARMA, and SDARMA models would suggest positive (or negative) spatial lag dependence. Similarly, a significant and positive (or negative) value of lamda in the SEM, SDEM, SAC, and SDAC models would suggest positive (or negative) spatial error dependence, which is contrary to the SMA, SDMA, SARMA, and SDARMA models. For the SMA, SDMA, SARMA, and SDARMA models, a significant and positive (or negative) value of lamda would suggest negative (or positive) spatial error dependence.

Last updated: June 19, 2025