Local–Linear Two-Stage Estimation of Local Autoregressive Geographically and Temporally Weighted Regression Model

A geographically and temporally weighted regression (GTWR) model is an effective tool for dealing with spatial heterogeneity and temporal non-stationarity simultaneously. As an important characteristic of spatiotemporal data, spatiotemporal autocorrelation should be considered when constructing spatiotemporally varying coefficient models. The proposed local autoregressive geographically and temporally weighted regression (GTWRLAR) model can simultaneously handle spatiotemporal autocorrelations among response variables and the spatiotemporal heterogeneity of regression relationships. The two-stage weighted least squares (2SLS) estimation can effectively reduce computational complexity. However, the weighted least squares estimation is essentially a Nadaraya–Watson kernel-smoothing approach for nonparametric regression models, and it suffers from a boundary effect. For spatiotemporally varying coefficient models, the three-dimensional spatiotemporal coefficients (longitude, latitude, and time) inherently exhibit larger boundaries than one-dimensional intervals. Therefore, the boundary effect of the 2SLS estimation of GTWRLAR will be more serious. A local–linear geographically and temporally weighted 2SLS (GTWRLAR-L) estimation is proposed to correct the boundary effect in both the spatial and temporal dimensions of GTWRLAR and simultaneously improve parameter estimation accuracy. The simulation experiment shows that the GTWRLAR-L method reduces the root mean square error (RMSE) of parameter estimates compared to the standard GTWRLAR approach. Empirical analyses of carbon emissions in China’s Yellow River Basin (2017–2021) show that GTWRLAR-L enhances the adjusted R2 from 0.888 to 0.893.