Beyond Gaussian processes: Flexible Bayesian modeling and inference for geostatistical processes
This work proposes a novel family of geostatistical models to account for features that cannot be properly accommodated by traditional Gaussian processes. The family is specified hierarchically and combines the infinite dimensional dynamics of Gaussian processes to that of any multivariate continuous distribution. This combination is stochastically defined through a latent Poisson process and the new family is called the Poisson-Gaussian Mixture Process - POGAMP. Whilst the attempt of defining a geostatistical process by assigning some arbitrary continuous distributions to be the finite-dimension distributions usually leads to non-valid processes, the POGAMP can have its finite-dimensional distributions to be arbitrarily close to any continuous distribution and still be a valid process. Formal results to establish its existence and other important properties, such as absolute continuity w.r.t. a Gaussian process measure, are provided. Also, a MCMC algorithm is carefully devised to perform Bayesian inference when the POGAMP is discretely observed in some space domain. Simulations are performed to empirically investigate the modelling properties of the POGAMP and the efficiency of the MCMC algorithm. Finally, some real datasets are analysed to illustrate the applicability of the proposed methodology.
READ FULL TEXT