An evaluation of likelihood-based bandwidth selectors for spatial and spatiotemporal kernel estimates

ABSTRACT Spatial point pattern data sets are commonplace in a variety of different research disciplines. The use of kernel methods to smooth such data is a flexible way to explore spatial trends and make inference about underlying processes without, or perhaps prior to, the design and fitting of more intricate semiparametric or parametric models to quantify specific effects. The long-standing issue of ‘optimal’ data-driven bandwidth selection is complicated in these settings by issues such as high heterogeneity in observed patterns and the need to consider edge correction factors. We scrutinize bandwidth selectors built on leave-one-out cross-validation approximation to likelihood functions. A key outcome relates to previously unconsidered adaptive smoothing regimens for spatiotemporal density and multitype conditional probability surface estimation, whereby we propose a novel simultaneous pilot-global selection strategy. Motivated by applications in epidemiology, the results of both simulated and real-world analyses suggest this strategy to be largely preferable to classical fixed-bandwidth estimation for such data.

[1]  Andrew B. Lawson,et al.  An evaluation of non-parametric relative risk estimators for disease maps , 2004, Comput. Stat. Data Anal..

[2]  G. Terrell The Maximal Smoothing Principle in Density Estimation , 1990 .

[3]  Matthew P. Wand,et al.  Kernel Smoothing , 1995 .

[4]  Emanuele Giorgi,et al.  Spatial point patterns:methodology and applications with R , 2017 .

[5]  P. Diggle A Kernel Method for Smoothing Point Process Data , 1985 .

[6]  J. S. Marron,et al.  Variable window width kernel estimates of probability densities , 1988 .

[7]  Martin L. Hazelton,et al.  On the utility of asymptotic bandwidth selectors for spatially adaptive kernel density estimation , 2018, Statistics & Probability Letters.

[8]  Shuowen Hu,et al.  Bayesian adaptive bandwidth kernel density estimation of irregular multivariate distributions , 2012, Comput. Stat. Data Anal..

[9]  Martin L Hazelton,et al.  Generalizing the spatial relative risk function. , 2014, Spatial and spatio-temporal epidemiology.

[10]  Martin L. Hazelton,et al.  Symmetric adaptive smoothing regimens for estimation of the spatial relative risk function , 2016, Comput. Stat. Data Anal..

[11]  Bernard W. Silverman,et al.  Density Estimation for Statistics and Data Analysis , 1987 .

[12]  Peter J. Diggle,et al.  Spatial and Space-Time Point Pattern Analysis , 2015 .

[13]  A B Lawson,et al.  Spatial competing risk models in disease mapping. , 2000, Statistics in medicine.

[14]  Gerard Rushton,et al.  Modeling the probability distribution of positional errors incurred by residential address geocoding , 2007 .

[15]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[16]  P. Diggle,et al.  Non-parametric estimation of spatial variation in relative risk. , 1995, Statistics in medicine.

[17]  Martin L. Hazelton,et al.  Assessing minimum contrast parameter estimation for spatial and spatiotemporal log-Gaussian Cox processes , 2013 .

[18]  Peter J. Diggle,et al.  Point process methodology for on‐line spatio‐temporal disease surveillance , 2005 .

[19]  Tilman M. Davies,et al.  Fast computation of spatially adaptive kernel estimates , 2017, Statistics and Computing.

[20]  M. C. Jones,et al.  Simple boundary correction for kernel density estimation , 1993 .

[21]  Modelling dichotomously marked muscle fibre configurations. , 2013, Statistics in medicine.

[22]  A. Cuevas,et al.  A comparative study of several smoothing methods in density estimation , 1994 .

[23]  Dai Feng,et al.  Computing and Displaying Isosurfaces in R , 2008 .

[24]  A B Lawson,et al.  Applications of extraction mapping in environmental epidemiology. , 1993, Statistics in medicine.

[25]  J. Møller,et al.  Log Gaussian Cox Processes , 1998 .

[26]  Peter J. Diggle,et al.  Nonparametric estimation of spatial segregation in a multivariate point process: bovine tuberculosis in Cornwall, UK , 2005 .

[27]  Rob J. Hyndman,et al.  A Bayesian approach to bandwidth selection for multivariate kernel density estimation , 2006, Comput. Stat. Data Anal..

[28]  D. Wheeler A comparison of spatial clustering and cluster detection techniques for childhood leukemia incidence in Ohio, 1996 – 2003 , 2007, International journal of health geographics.

[29]  P. J. Green,et al.  Density Estimation for Statistics and Data Analysis , 1987 .

[30]  Ian Abramson On Bandwidth Variation in Kernel Estimates-A Square Root Law , 1982 .

[31]  Kernel estimation of risk surfaces without the need for edge correction , 2008, Statistics in medicine.

[32]  Martin L. Hazelton,et al.  Boundary kernels for adaptive density estimators on regions with irregular boundaries , 2010, J. Multivar. Anal..

[33]  A. Gelfand,et al.  ANALYSIS OF MINNESOTA COLON AND RECTUM CANCER POINT PATTERNS WITH SPATIAL AND NONSPATIAL COVARIATE INFORMATION. , 2009, The annals of applied statistics.

[34]  Tilman M. Davies,et al.  Inference Based on Kernel Estimates of the Relative Risk Function in Geographical Epidemiology , 2009, Biometrical journal. Biometrische Zeitschrift.

[35]  Frederik P. Agterberg,et al.  Interactive spatial data analysis , 1996 .

[36]  Tilman M. Davies,et al.  Tutorial on kernel estimation of continuous spatial and spatiotemporal relative risk , 2017, Statistics in medicine.

[37]  D. W. Scott,et al.  Cross-Validation of Multivariate Densities , 1994 .

[38]  P. Diggle,et al.  Kernel estimation of relative risk , 1995 .

[39]  J. Bithell Estimation of relative risk functions. , 1991, Statistics in medicine.

[40]  P. Diggle,et al.  Spatial variation in risk of disease: a nonparametric binary regression approach , 2002 .

[41]  Peter J. Diggle,et al.  Edge-correction for spatial kernel smoothing methods? When is it necessary? , 2004 .

[42]  J. Bithell An application of density estimation to geographical epidemiology. , 1990, Statistics in medicine.

[43]  Andrew B. Lawson,et al.  Armadale: A Case‐Study in Environmental Epidemiology , 1994 .

[44]  Tilman M. Davies,et al.  Identification of high-risk regions for schistosomiasis in the Guichi region of China: an adaptive kernel density estimation-based approach , 2013, Parasitology.

[45]  A. Gelfand,et al.  Spatial Point Patterns , 2010 .

[46]  M. C. Jones,et al.  A Brief Survey of Bandwidth Selection for Density Estimation , 1996 .

[47]  Tilman M. Davies,et al.  Adaptive kernel estimation of spatial relative risk , 2010, Statistics in medicine.