spaMM (version 1.7.2)

adjacency: Fitting autoregressive models

Description

This describes the autoregressive models (AR) considered by spaMM and how to fit them. This documentation is mainly a stub, as more implementation of some these models (in particular, of simultaneous AR, i.e., SAR) is to be expected in the future. The only model currently implemented is the adjacency model (a conditional AR, i.e., CAR).

Details

The most widely applicable fitting strategy in spaMM is to fit possibly complex spatial models by numerical maximization over correlation parameters. However, more efficient algorithms are available for CAR and SAR models, and have been widely used in particular in the econometric literature (e.g., LeSage and Pace 2009). Currently such an algorithm is implemented only for one CAR model, the adjacency model. The taxonomy of AR models is not well settled, but adjacency should retain its current meaning as a specific CAR model in future versions of spaMM. In the adjacency model, the covariance matrix of random effects u can be described as $\lambda$(I$-\rho$ W$)^{-1}$ where W is the (symmetric) adjacency matrix. If you use corrHLfit, the default fitting method for the adjacency model is numerical maximization of the likelihood (or restricted likelihood) as function of the correlation parameter $\rho$. However, a call to HLCor fits more efficiently the variance and correlation parameters of random effects (see examples). The ML fits by both methods should be practically equivalent. The REML fits should slightly differ from each other, due to the fact that the REML approximation for GLMMs does not maximize a single likelihood function. Finally, a call to corrHLfit with the additional argument init.HLfit=list(rho=0) should be equivalent in speed and result to the HLCor call. The efficient method uses the spectral decomposition W=VDV' where D is a diagonal matrix of eigenvalues $d_i$, the covariance of V'u is $\lambda$(I$-\rho$ D$)^{-1}$, which is a diagonal matrix with elements $\lambda_i$=$\lambda$/(1$-\rho d_i$). Hence $1/\lambda_i$ is in the linear predictor form $\alpha$+$\beta d_i$ This can be used to fit $\lambda$ and $\rho$ efficiently. If the efficient fitting method is used, the results are reported as the coefficients $\alpha$ ((Intercept)) and $\beta$ (adjd) of the predictor for $1/\lambda_i$, in addition to the resulting values of $\rho$ and of the common $\lambda$ factor.

References

LeSage, J., Pace, R.K. (2009) Introduction to Spatial Econometrics. Chapman & Hall/CRC.

Examples

data(scotlip)
# CAR by Laplace with 'outer' estimation of rho
if (spaMM.getOption("example_maxtime")>5) {          
  corrHLfit(cases~I(prop.ag/10) +adjacency(1|gridcode)+offset(log(scotlip$expec)),
          adjMatrix=Nmatrix,family=poisson(),data=scotlip,HLmethod="ML") ## ~ 5 s.
}

# CAR by Laplace with 'inner' estimation of rho
HLCor(cases~I(prop.ag/10) +adjacency(1|gridcode)+offset(log(scotlip$expec)),
          adjMatrix=Nmatrix,family=poisson(),data=scotlip,HLmethod="ML") ## < 1 s.