On comparing some algorithms for finding the optimal bandwidth in Geographically Weighted Regression

Abstract Geographically Weighted Regression (GWR) models are sensitive to the choice of the bandwidth, and this choice is traditionally made through the golden section search algorithm. This algorithm is applied in a function, known as Cross-Validation (CV), which quantifies the efficiency of the model, therefore looking for the optimal parameter that results in the best model. In this paper, the behavior of the CV function was studied, and it was verified that when it is not strictly convex, the golden section search algorithm converges to local minimums instead of the global one. Three algorithms have been used to find the optimal bandwidth: the lightning search algorithm, the harmony search algorithm and an adaptation of the golden section search algorithm. In addition, comparisons were made between them to check the suitability of each one in GWR models. It was found that the golden section search algorithm is not the most adequate in this situation because, in more than one simulation, it resulted in a value too far from the optimal bandwidth. It was also verified that the models with the bandwidth far from the optimal value showed differences in the significance of the parameter estimates compared to the models with the optimal bandwidth.

[1]  M. D. Noskov,et al.  Modeling the development of the stepped leader of a lightning discharge , 1999 .

[2]  A S Fotheringham,et al.  Geographically weighted Poisson regression for disease association mapping , 2005, Statistics in medicine.

[3]  A. Fotheringham,et al.  The Multiple Testing Issue in Geographically Weighted Regression , 2016 .

[4]  Zong Woo Geem,et al.  A New Heuristic Optimization Algorithm: Harmony Search , 2001, Simul..

[5]  David G. Luenberger,et al.  Linear and nonlinear programming , 1984 .

[6]  William H. Press,et al.  Numerical recipes , 1990 .

[7]  Alan Ricardo da Silva,et al.  Geographically Weighted Negative Binomial Regression—incorporating overdispersion , 2014, Stat. Comput..

[8]  Chris Brunsdon,et al.  Geographically Weighted Regression: The Analysis of Spatially Varying Relationships , 2002 .

[9]  A. Berkopec Fast particles as initiators of stepped leaders in CG and IC lightnings , 2012 .

[10]  Gaige Wang,et al.  A Novel Hybrid Bat Algorithm with Harmony Search for Global Numerical Optimization , 2013, J. Appl. Math..

[11]  L. Zhang,et al.  Comparison of bandwidth selection in application of geographically weighted regression : a case study , 2008 .

[12]  Steven Farber,et al.  A systematic investigation of cross-validation in GWR model estimation: empirical analysis and Monte Carlo simulations , 2007, J. Geogr. Syst..

[13]  S. Fotheringham,et al.  Geographically Weighted Regression , 1998 .

[14]  Azah Mohamed,et al.  A binary variant of lightning search algorithm: BLSA , 2017, Soft Comput..

[15]  S. Fotheringham,et al.  Geographically weighted regression : modelling spatial non-stationarity , 1998 .

[16]  Hussain Shareef,et al.  Lightning search algorithm , 2015, Appl. Soft Comput..