Finite area smoothing with generalized distance splines

Environmental and Ecological Statistics - Tập 21 - Trang 715-731 - 2014
David L. Miller1, Simon N. Wood2
1Centre for Research into Ecological and Environmental Modelling, University of St Andrews, The Observatory, Scotland
2Department of Mathematical Sciences, University of Bath, Claverton Down, Bath, UK

Tóm tắt

Most conventional spatial smoothers smooth with respect to the Euclidean distance between observations, even though this distance may not be a meaningful measure of spatial proximity, especially when boundary features are present. When domains have complicated boundaries leakage (the inappropriate linking of parts of the domain which are separated by physical barriers) can occur. To overcome this problem, we develop a method of smoothing with respect to generalized distances, such as within domain distances. We obtain the generalized distances between our points and then use multidimensional scaling to find a configuration of our observations in a Euclidean space of 2 or more dimensions, such that the Euclidian distances between points in that space closely approximate the generalized distances between the points. Smoothing is performed over this new point configuration, using a conventional smoother. To mitigate the problems associated with smoothing in high dimensions we use a generalization of thin plate spline smoothers proposed by Duchon (Constructive theory of functions of several variables, pp 85–100, 1977). This general method for smoothing with respect to generalized distances improves on the performance of previous within-domain distance spatial smoothers, and often provides a more natural model than the soap film approach of Wood et al. (J R Stat Soc Ser B Stat Methodol 70(5):931–955, 2008). The smoothers are of the linear basis with quadratic penalty type easily incorporated into a range of statistical models.

Tài liệu tham khảo

Bernstein M, De Silva V, Langford J, Tenenbaum J (2000) Graph approximations to geodesics on embedded manifolds. Technical report, Department of Psychology, Stanford University. ftp://ftp-sop.inria.fr/prisme/boissonnat/ImageManifolds/isomap.pdf

Curriero F (2005) On the use of non-euclidean isotropy in geostatistics. Technical report 94, Johns Hopkins University, Department of Biostatistics. http://www.bepress.com/cgi/viewcontent.cgi?article=1094&context=jhubiostat

Driscoll TA, Trefethen L (2002) Schwartz-Christoffel transform. Cambridge University Press, Cambridge

Duchon J (1977) Splines minimizing rotation-invariant semi-norms in Sobolev spaces. Constructive theory of functions of several variables, pp 85–100

Hastie TJ, Tibshirani RJ (1990) Generalized additive models. Monographs on statistics and applied probability. Taylor & Francis, New York

Jensen OP, Christman MC, Miller TJ (2006) Landscape-based geostatistics: a case study of the distribution of blue crab in Chesapeake Bay. Environmetrics 17(6):605–621. doi:10.1002/env.767

Løland A, Høst G (2003) Spatial covariance modelling in a complex coastal domain by multidimensional scaling. Environmetrics 14(3):307–321. doi:10.1002/env.588

Oh MS, Raftery AE (2001) Bayesian multidimensional scaling and. J Am Stat Assoc 96(455):1031

Ruppert D, Wand M, Carroll RJ (2003) Semiparametric regression. Cambridge series on statistical and probabilistic mathematics. Cambridge University Press, Cambridge

Scott-Hayward LAS, MacKenzie ML, Donovan CR, Walker CG, Ashe E (2013) Complex region spatial smoother (CReSS). J Comput Graph. Stat. doi: 10.1080/10618600.2012.762920

Vretblad A (2003) Fourier analysis and its applications. Graduate texts in mathematics. Springer, Berlin

Williams R, Hedley SL, Branch TA, Bravington MV, Zerbini AN, Findlay KP (2011) Chilean blue whales as a case study to illustrate methods to estimate abundance and evaluate conservation status of rare species. Conserv Biol 25(3):526–535. doi:10.1111/j.1523-1739.2011.01656.x