Related papers: BeyondPlanck II. CMB map-making through Gibbs samp…
We describe a Bayesian framework for estimating the time-domain noise covariance of CMB observations, typically parametrized in terms of a 1/f frequency profile. This framework is based on the Gibbs sampling algorithm, which allows for…
In this paper we describe how MAP inference can be used to sample efficiently from Gibbs distributions. Specifically, we provide means for drawing either approximate or unbiased samples from Gibbs' distributions by introducing low…
We present a Bayesian calibration algorithm for CMB observations as implemented within the global end-to-end BeyondPlanck (BP) framework, and apply this to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck…
Map-making presents a significant computational challenge to the next generation of kilopixel CMB polarisation experiments. Years worth of time ordered data (TOD) from thousands of detectors will need to be compressed into maps of the T, Q…
We present a new map-making method for CMB measurements. The method is based on the destriping technique, but it also utilizes information about the noise spectrum. The low-frequency component of the instrument noise stream is modelled as a…
We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BeyondPlanck, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis.…
We study different variants of the Gibbs sampler algorithm from the perspective of their applicability to the estimation of power spectra of the cosmic microwave background (CMB) anisotropies. These include approaches studied earlier in the…
We propose a solution to the CMB component separation problem based on standard parameter estimation techniques. We assume a parametric spectral model for each signal component, and fit the corresponding parameters pixel by pixel in a…
CMB data analysis is in general done through two main steps : map-making of the time data streams and power spectrum extraction from the maps. The latter basically consists in the separation between the variance of the CMB and that of the…
Line-intensity mapping (LIM) is an emerging cosmological technique that traces large-scale structure through the integrated spectral-line emission of unresolved sources. Reconstructing unbiased sky maps requires careful joint treatment of…
A central task in many applications is reasoning about processes that change over continuous time. Continuous-Time Bayesian Networks is a general compact representation language for multi-component continuous-time processes. However, exact…
This paper derives two new optimization-driven Monte Carlo algorithms inspired from variable splitting and data augmentation. In particular, the formulation of one of the proposed approaches is closely related to the alternating direction…
This paper addresses the mapping problem. Using a conjugate prior form, we derive the exact theoretical batch multi-object posterior density of the map given a set of measurements. The landmarks in the map are modeled as extended objects,…
We present in this article two different ways to make CMB maps in practice, from large timelines. One is to make a simple destriping, fitting the data and using the scan intercepts to remove the low frequency noise (stripes). The second,…
Monte Carlo methods are essential tools for Bayesian inference. Gibbs sampling is a well-known Markov chain Monte Carlo (MCMC) algorithm, extensively used in signal processing, machine learning, and statistics, employed to draw samples from…
The Hamiltonian Monte Carlo (HMC) algorithm is often lauded for its ability to effectively sample from high-dimensional distributions. In this paper we challenge the presumed domination of HMC for the Bayesian analysis of GLMs. By utilizing…
We present a parallel implementation of a map-making algorithm for CMB anisotropy experiments which is both fast and efficient. We show for the first time a Maximum Likelihood, minimum variance map obtained by processing the entire data…
Future cosmic microwave background (CMB) polarisation experiments aim to measure an unprecedentedly small signal - the primordial gravity wave component of the polarisation field B-mode. To achieve this, they will analyse huge datasets,…
A fundamental task in machine learning and related fields is to perform inference on Bayesian networks. Since exact inference takes exponential time in general, a variety of approximate methods are used. Gibbs sampling is one of the most…
We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are 1) conditional…