Abstract
To improve the quality of prediction of radioactive contamination, geostatistical methods, and in particular multivariate geostatistical models, are increasingly being used. These methods, however, are optimal only in the case in which the data may be assumed Gaussian and do not properly cope with data measurements that are discrete, nonnegative or show some degree of skewness. To deal with these situations, here we consider a hierarchical model in which non-Gaussian variables of different kind are handled simultaneously. We show that when observations are assumed to be conditionally distributed as Poisson and Gamma, variograms and cross-variograms have convenient simple forms, and estimation of the parameters of the model can be carried out by Monte Carlo EM. This work was inspired by radioactive contamination data from the Maddalena Archipelago (Sardinia, Italy).
Original language | English |
---|---|
Pages (from-to) | 195-213 |
Number of pages | 19 |
Journal | AStA Advances in Statistical Analysis |
Volume | 97 |
Issue number | 2 |
DOIs | |
Publication status | Published - 2013 |
Keywords
- Generalized linear mixed model
- Linear model of coregionalization
- Markov chain Monte Carlo
- Monte Carlo EM
- Spatial factor model
ASJC Scopus subject areas
- Analysis
- Statistics and Probability
- Modelling and Simulation
- Social Sciences (miscellaneous)
- Economics and Econometrics
- Applied Mathematics