Related papers: The Information Content of Projected Galaxy Fields
The normalisation of the matter power spectrum, sigma_8, is an essential ingredient to predict the phenomenology of the low redshift universe. It has been measured using several methods, such as X-ray cluster counts, weak lensing and the…
One of the major science goals over the coming decade is to test fundamental physics with probes of the cosmic large-scale structure out to high redshift. Here we present a fully Bayesian approach to infer the three-dimensional cosmic…
Markov chain Monte Carlo (MCMC) sampling of densities restricted to linearly constrained domains is an important task arising in Bayesian treatment of inverse problems in the natural sciences. While efficient algorithms for uniform polytope…
We present the first cosmological measurement derived from a galaxy density field subject to a `clipping' transformation. By enforcing an upper bound on the galaxy number density field in the Galaxy and Mass Assembly survey (GAMA),…
Extracting the non-Gaussian information of the cosmic large-scale structure (LSS) is vital in unlocking the full potential of the rich datasets from the upcoming stage-IV galaxy surveys. Galaxy skew spectra serve as efficient…
Plenty of crucial information about our Universe is encoded in the cosmic large-scale structure (LSS). However, the extractions of these information are usually hindered by the nonlinearities of the LSS, which can be largely alleviated by…
We present a Hamiltonian Monte Carlo algorithm to sample from multivariate Gaussian distributions in which the target space is constrained by linear and quadratic inequalities or products thereof. The Hamiltonian equations of motion can be…
High-dimensional data sets are expected from the next generation of large-scale surveys. These data sets will carry a wealth of information about the early stages of galaxy formation and cosmic reionization. Extracting the maximum amount of…
In this paper, we discuss an extension of the Split Hamiltonian Monte Carlo (Split HMC) method for Gaussian process model (GPM). This method is based on splitting the Hamiltonian in a way that allows much of the movement around the state…
A new method for the reconstruction of the projected mass distribution of clusters of galaxies from the image distortion of background galaxies is discussed. This method is essentially equivalent to the one we developed previously, i.e.,…
The generalized inverse Gaussian-Poisson (GIGP) distribution proposed by Sichel in the 1970s has proved to be a flexible fitting tool for diverse frequency data, collectively described using the item production model. In this paper, we…
We present a new approach to calculating the statistical distributions for magnification, shear, and rotation of images of cosmological sources due to gravitational lensing. In this approach one specifies an underlying Robertson-Walker…
In this paper, we describe a procedure for modelling strong lensing galaxy clusters with parametric methods, and to rank models quantitatively using the Bayesian evidence. We use a publicly available Markov chain Monte-Carlo (MCMC) sampler…
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…
Weak gravitational lensing induces flux dependent fluctuations in the observed galaxy number density distribution. This cosmic magnification (magnification bias) effect in principle enables lensing reconstruction alternative to cosmic shear…
We introduce a variant of the Hybrid Monte Carlo (HMC) algorithm to address large-deviation statistics in stochastic hydrodynamics. Based on the path-integral approach to stochastic (partial) differential equations, our HMC algorithm…
The redistribution of baryonic matter in massive halos through processes like active galactic nuclei feedback and star formation leads to a suppression of the matter power spectrum on small scales. This redistribution can be measured…
Interest rises to exploit the full shape information of the galaxy power spectrum, as well as pushing analyses to smaller non-linear scales. Here I use the halo model to quantify the information content in the tomographic angular power…
Compressed sensing is a theory which guarantees the exact recovery of sparse signals from a small number of linear projections. The sampling schemes suggested by current compressed sensing theories are often of little practical relevance…
We consider the problem of recovering $n$ i.i.d samples from a zero mean multivariate Gaussian distribution with an unknown covariance matrix, from their modulo wrapped measurements, i.e., measurement where each coordinate is reduced modulo…