Atelier réunissant astrophysiciens et statisticiens
8-9 déc. 2011 Grenoble (France)
Estimation of cosmological parameters using adaptive importance sampling
Darren Wraith, Martin Kilbinger, Karim Benabed, Olivier Cappé, Jean-François Cardoso, Gersende Fort, Simon Prunet, Christian Robert  1  
1 : CEntre de REcherches en MAthématiques de la DEcision  (CEREMADE)  -  Site web
CNRS : UMR7534, Université Paris Dauphine - Paris IX
Place du Maréchal de Lattre de Tassigny 75775 - Paris Cedex 16 -  France

We present a Bayesian sampling algorithm called adaptive importance sampling or Population Monte Carlo (PMC), whose computational workload is easily parallelizable and thus has the potential to considerably reduce the wall-clock time required for sampling, along with providing other benefits. To assess the performance of the approach for cosmological problems, we use simulated and
actual data consisting of CMB anisotropies, supernovae of type Ia, and weak cosmological lensing, and provide a comparison of results to those obtained using state-of-the-art Markov Chain Monte Carlo (MCMC). For both types of data sets, we find comparable parameter estimates for PMC and MCMC, with the advantage of a significantly lower computational time for PMC. In the case of
WMAP5 data, for example, the wall-clock time reduces from several days for MCMC to a few hours using PMC on a cluster of processors. Other benefits of the PMC approach, along with potential difficulties in using the approach, are analysed and discussed.

Personnes connectées : 1