Bayesian data augmentation reintroduces missing LISA data segments as auxiliary variables during posterior sampling to enable consistent parameter estimation for galactic binaries despite gaps.
Bayesian and Maximum Likelihood Estimation for Gaussian Processes on an Incomplete Lattice
1 Pith paper cite this work. Polarity classification is still indexing.
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.
fields
gr-qc 1years
2019 1verdicts
UNVERDICTED 1representative citing papers
citing papers explorer
-
Gravitational-wave parameter estimation with gaps in LISA: a Bayesian data augmentation method
Bayesian data augmentation reintroduces missing LISA data segments as auxiliary variables during posterior sampling to enable consistent parameter estimation for galactic binaries despite gaps.