Related papers: Fast optimal CMB power spectrum estimation with Ha…
Hamiltonian Monte Carlo (HMC) is a Markov chain algorithm for sampling from a high-dimensional distribution with density $e^{-f(x)}$, given access to the gradient of $f$. A particular case of interest is that of a $d$-dimensional Gaussian…
Recently, the Hamilton Monte Carlo (HMC) has become widespread as one of the more reliable approaches to efficient sample generation processes. However, HMC is difficult to sample in a multimodal posterior distribution because the HMC chain…
We improve the algorithm of Kosowsky, Milosavljevic, and Jimenez (2002) for computing power spectra of the cosmic microwave background. The present algorithm computes not only the temperature power spectrum but also the E-mode polarization…
Efficient sampling from high-dimensional distributions is a challenging issue which is encountered in many large data recovery problems involving Markov chain Monte Carlo schemes. In this context, sampling using Hamiltonian dynamics is one…
Hamiltonian Monte Carlo (HMC) is an efficient Bayesian sampling method that can make distant proposals in the parameter space by simulating a Hamiltonian dynamical system. Despite its popularity in machine learning and data science, HMC is…
A maximum-likelihood method is presented for estimating the power spectrum of anisotropies in the cosmic microwave background (CMB) from interferometer observations. The method calculates flat band-power estimates in separate bins in…
Estimation of the angular power spectrum is one of the important steps in Cosmic Microwave Background (CMB) data analysis. Here, we present a nonparametric estimate of the temperature angular power spectrum for the Planck 2013 CMB data. The…
We use a Quadratic Maximum Likelihood (QML) method to estimate the angular power spectrum of the cross-correlation between cosmic microwave background and large scale structure maps as well as their individual auto-spectra. We describe our…
Szapudi et al (2001) introduced the method of estimating angular power spectrum of the CMB sky via heuristically weighted correlation functions. Part of the new technique is that all (co)variances are evaluated by massive Monte Carlo…
We describe an algorithm for the extraction of the angular power spectrum of an intensity field, such as the cosmic microwave background (CMB), from interferometer data. This new method, based on the gridding of interferometer visibilities…
We describe a fast and accurate method for estimation of the cosmic microwave background (CMB) anisotropy angular power spectrum --- Monte Carlo Apodised Spherical Transform EstimatoR. Originally devised for use in the interpretation of the…
In Hamiltonian Monte Carlo sampling, the shape of the potential and the choice of the momentum distribution jointly give rise to the Hamiltonian dynamics of the sampler. An efficient sampler propagates quickly in all regions of the…
We present two novel methods for the estimation of the angular power spectrum of cosmic microwave background (CMB) anisotropies. We assume an absolute CMB experiment with arbitrary asymmetric beams and arbitrary sky coverage. The methods…
Pulsar timing arrays (PTAs) detect low-frequency gravitational waves (GWs) by looking for correlated deviations in pulse arrival times. Current Bayesian searches use Markov Chain Monte Carlo (MCMC) methods, which struggle to sample the…
This paper generalises the hybrid power spectrum estimator developed in Efstathiou (2004a) to the estimation of polarization power spectra of the cosmic microwave background radiation. The hybrid power spectrum estimator is unbiased and we…
Power spectrum estimation and evaluation of associated errors in the presence of incomplete sky coverage; non-homogeneous, correlated instrumental noise; and foreground emission is a problem of central importance for the extraction of…
It is widely believed that maximum likelihood estimators must be used to provide optimal estimates of power spectra. Since such estimators require require of order N_d^3 operations they are computationally prohibitive for N_d greater than a…
We develop the XFaster Cosmic Microwave Background (CMB) temperature and polarization anisotropy power spectrum and likelihood technique for the Planck CMB satellite mission. We give an overview of this estimator and its current…
Hamiltonian Monte Carlo (HMC) is a state of the art method for sampling from distributions with differentiable densities, but can converge slowly when applied to challenging multimodal problems. Running HMC with a time varying Hamiltonian,…
We test the hemispherical power asymmetry of the WMAP 7-year low-resolution temperature and polarization maps. We consider two natural estimators for such an asymmetry and exploit our implementation of an optimal angular power spectrum…