Abstract

We propose a hierarchical model coupled to geostatistics to deal with a non-gaussian data distribution and take explicitly into account complex spatial structures (i.e. trends, patchiness and random fluctuations). A common characteristic of animal count data is a distribution that is both zero-inflated and heavy tailed. In such cases, empirical variograms are no more robust and most structural analyses result in poor and noisy estimated spatial variogram structures. Thus kriged maps feature a broad variance of prediction. Moreover, due to the heterogeneity of wildlife population habitats, a nonstationary model is often required. To avoid these difficulties, we propose a hierarchical model that assumes that the count data follow a Poisson distribution given a theoretical sighting density which is a latent variable to be estimate. This density is modelled as the product of a positive long range trend by a positive stationary random field, characterized by a unit mean and a variogram function. A first estimate of the drift is used to obtain an estimate of the variogram of residuals including a correction term for variance coming from the Poisson distribution and weights due to the non-constant spatial mean. Then a kriging procedure similar to a modified universal kriging is implemented to directly map the latent density from raw count data. An application on fin whale data illustrates the effectiveness of the method in mapping animal density in a context that is presumably non-stationary. E. Bellier and P. Monestiez Biostatistique et Processus Spatiaux, INRA, Domaine Saint-Paul, Site Agroparc, 84914 Avignon cedex 9, France E. Bellier ( ) Norwegian Institute for Nature Research NINA, NO-7485 Trondheim, NORWAY e-mail: edwige.bellier@nina.no C. Guinet Centre d’Etudes Biologiques de Chize, CNRS, 79360 Villiers-en-Bois, France P.M. Atkinson and C.D. Lloyd (eds.), geoENV VII – Geostatistics for Environmental Applications, Quantitative Geology and Geostatistics 16, DOI 10.1007/978-90-481-2322-3 1, c Springer Science+Business Media B.V. 2010 1

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call