Related papers: Non Gaussian extrema counts for CMB maps
A multivariate distribution can be described by a triangular transport map from the target distribution to a simple reference distribution. We propose Bayesian nonparametric inference on the transport map by modeling its components using…
Deriving Bayesian inference for exponential random graph models (ERGMs) is a challenging "doubly intractable" problem as the normalizing constants of the likelihood and posterior density are both intractable. Markov chain Monte Carlo (MCMC)…
Since the inhomogeneous instrument noise can produce extra non-Gaussianity in the CMB anisotropy, its effect should be carefully subtracted in the primordial non-Gaussianity estimation. We calculate the probability distribution function of…
The recently proposed non-Gaussian Mat\'{e}rn random field models, generated through Stochastic Partial differential equations (SPDEs), are extended by considering the class of Generalized Hyperbolic processes as noise forcings. The models…
The non-Gaussianity of inflationary perturbations, as encoded in the bispectrum (or 3-point correlator), has become an important additional way of distinguishing between inflation models, going beyond the linear Gaussian perturbation…
While the probability density function (PDF) of the cosmic microwave background (CMB) convergence field approximately follows a Gaussian distribution, primordial non-Gaussianities and small contributions from structures at low redshifts…
Bayesian inference for models that have an intractable partition function is known as a doubly intractable problem, where standard Monte Carlo methods are not applicable. The past decade has seen the development of auxiliary variable Monte…
We develop a novel computational method for evaluating the extreme excursion probabilities arising from random initialization of nonlinear dynamical systems. The method uses excursion probability theory to formulate a sequence of Bayesian…
We develop a biased Monte Carlo algorithm to measure probabilities of rare events in cluster-cluster aggregation for arbitrary collision kernels. Given a trajectory with a fixed number of collisions, the algorithm modifies both the waiting…
In cosmic web analysis, complementary to traditional cosmological probes, the extrema (e.g. peaks and voids) two-point correlation functions (2PCFs) are of particular interest for the study of both astrophysical phenomena and cosmological…
Our focus is on the design and analysis of efficient Monte Carlo methods for computing tail probabilities for the suprema of Gaussian random fields, along with conditional expectations of functionals of the fields given the existence of…
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…
Peak counts have been shown to be an excellent tool to extract the non-Gaussian part of the weak lensing signal. Recently, we developped a fast stochastic forward model to predict weak-lensing peak counts. Our model is able to reconstruct…
Analytic expressions for the statistics of peaks of random fields with weak non-Gaussianity are provided. Specifically, the abundance and spatial correlation of peaks are represented by formulas which can be evaluated only by virtually…
We present a new method based on the N-point probability distribution (pdf) to study non-Gaussianity in cosmic microwave background (CMB) maps. Likelihood and Bayesian estimation are applied to a local non-linear perturbed model up to third…
The skew-spectrum statistic introduced by Munshi & Heavens (2010) has recently been used in studies of non-Gaussianity from diverse cosmological data sets including the detection of primary and secondary non-Gaussianity of Cosmic Microwave…
This paper deals with Gibbs samplers that include high dimensional conditional Gaussian distributions. It proposes an efficient algorithm that avoids the high dimensional Gaussian sampling and relies on a random excursion along a small set…
Non-Gaussian distributions of cosmic microwave background (CMB) anisotropies have been proposed to reconcile the discrepancies between different experiments at half-degree scales (Coulson et al. 1994). Each experiment probes a different…
We discuss the tests of nonGaussianity in CMB maps using morphological statistics known as Minkowski functionals. As an example we test degree-scale cosmic microwave background (CMB) anisotropy for Gaussianity by studying the \qmask map…
We focus on the problem of estimating and quantifying uncertainties on the excursion set of a function under a limited evaluation budget. We adopt a Bayesian approach where the objective function is assumed to be a realization of a Gaussian…