On bandwidth choice for spatial data density estimation. (21st April 2020)

• Record Type: Journal Article
• Authors: Jiang, Zhenyu; Ling, Nengxiang; Lu, Zudi; Tj⊘stheim, Dag; Zhang, Qiang
• Is Part Of: Journal of the Royal Statistical Society. Volume 82:Number 3(2020)
• Page Start: 817
• Page End: 840
• Publication Date: 2020-04-21
• Subjects: Cross‐validation -- Kernel density estimation -- Optimal bandwidth -- Spatial lattice data -- Spatially adaptive bandwidth choice; Statistics -- Periodicals; Great Britain -- Statistics -- Periodicals; 519.2
• DOI: 10.1111/rssb.12367
• Languages: English
• ISSNs: 1369-7412

Abstract

Summary: Bandwidth choice is crucial in spatial kernel estimation in exploring non‐Gaussian complex spatial data. The paper investigates the choice of adaptive and non‐adaptive bandwidths for density estimation given data on a spatial lattice. An adaptive bandwidth depends on local data and hence adaptively conforms with local features of the spatial data. We propose a spatial cross‐validation (SCV) choice of a global bandwidth. This is done first with a pilot density involved in the expression for the adaptive bandwidth. The optimality of the procedure is established, and it is shown that a non‐adaptive bandwidth choice comes out as a special case. Although the cross‐validation idea has been popular for choosing a non‐adaptive bandwidth in data‐driven smoothing of independent and time series data, its theory and application have not been much investigated for spatial data. For the adaptive case, there is little theory even for independent data. Conditions that ensure asymptotic optimality of the SCV‐selected bandwidth are derived, actually, also extending time series and independent data optimality results. Further, for the adaptive bandwidth with an estimated pilot density, oracle properties of the resultant density estimator are obtained asymptotically as if the true pilot were known. Numerical simulations show that finite sample performance of the SCV adaptive bandwidth choice works quite well. It outperforms the existing R routines such as the 'rule of thumb' and the so‐called 'second‐generation' Sheather–Jones bandwidths for moderate and big data sets. An empirical application to a set of spatial soil data is further implemented with non‐Gaussian features significantly identified.