pith. sign in

arxiv: 1402.4281 · v1 · pith:JEML4KC4new · submitted 2014-02-18 · 📊 stat.CO · stat.ME

Bayesian and Maximum Likelihood Estimation for Gaussian Processes on an Incomplete Lattice

classification 📊 stat.CO stat.ME
keywords inferencelikelihoodapproachbayesiandatamaximumestimationgaussian
0
0 comments X
read the original abstract

This paper proposes a new approach for Bayesian and maximum likelihood parameter estimation for stationary Gaussian processes observed on a large lattice with missing values. We propose an MCMC approach for Bayesian inference, and a Monte Carlo EM algorithm for maximum likelihood inference. Our approach uses data augmentation and circulant embedding of the covariance matrix, and provides exact inference for the parameters and the missing data. Using simulated data and an application to satellite sea surface temperatures in the Pacific Ocean, we show that our method provides accurate inference on lattices of sizes up to 512 x 512, and outperforms two popular methods: composite likelihood and spectral approximations.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Gravitational-wave parameter estimation with gaps in LISA: a Bayesian data augmentation method

    gr-qc 2019-07 unverdicted novelty 6.0

    Bayesian data augmentation reintroduces missing LISA data segments as auxiliary variables during posterior sampling to enable consistent parameter estimation for galactic binaries despite gaps.