Researcher profile

C. P. Robert

C. P. Robert contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
6works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2014arXiv

Constraining the thick disc formation scenario of the Milky Way

We study the shape of the thick disc using photometric data at high and intermediate latitudes from SDSS and 2MASS surveys. We use the population synthesis approach using an ABC-MCMC method to characterize the thick disc shape, scale height, scale length, local density and flare, and we investigate the extend of the thick disc formation period by simulating several formation episodes. We find that the vertical variation in density is not exponential, but much closer to an hyperbolic secant squared. Assuming a single formation epoch, the thick disc is better fitted with a sech2 scale height of 470 pc and a scale length of 2.3 kpc. However if one simulates two successive formation episodes, mimicking an extended formation period, the older episode has a higher scale height and a larger scale length than the younger episode, indicating a contraction during the collapse phase. The scale height decreases from 800 pc to 340 pc, and the scale length from 3.2 kpc to 2 kpc. The star formation increases from the old episode to the young one. During the fitting process, the halo parameters are also determined. The constraint on the halo shows that a transition between a inner and outer halo, if exists, cannot be at a distance of less than about 30 kpc, which is the limit of our investigation using turnoff halo stars. Finally, we show that extrapolating the thick disc towards the bulge region explains well the stellar populations observed there, that there is no longer need to invoke a classical bulge. To explain these results, the most probable scenario for the thick disc is that it formed while the Galaxy was gravitationally collapsing from well mixed gas-rich giant clumps sustained by high turbulence which awhile prevent a thin disc to form, as proposed by Bournaud et al. (2009). This scenario explains the observations in the thick disc region as well as in the bulge region. (abridged)

preprint2013arXiv

Objective Bayesian hypothesis testing in binomial regression models with integral prior distributions

In this work we apply the methodology of integral priors to handle Bayesian model selection in binomial regression models with a general link function. These models are very often used to investigate associations and risks in epidemiological studies where one goal is to exhibit whether or not an exposure is a risk factor for developing a certain disease; the purpose of the current paper is to test the effect of specific exposure factors. We formulate the problem as a Bayesian model selection case and solve it using objective Bayes factors. To construct the reference prior distributions on the regression coefficients of the binomial regression models, we rely on the methodology of integral priors that is nearly automatic as it only requires the specification of estimation reference priors and it does not depend on tuning parameters or on hyperparameters within these priors.

preprint2013arXiv

Relevant statistics for Bayesian model choice

The choice of the summary statistics used in Bayesian inference and in particular in ABC algorithms has bearings on the validation of the resulting inference. Those statistics are nonetheless customarily used in ABC algorithms without consistency checks. We derive necessary and sufficient conditions on summary statistics for the corresponding Bayes factor to be convergent, namely to asymptotically select the true model. Those conditions, which amount to the expectations of the summary statistics to asymptotically differ under both models, are quite natural and can be exploited in ABC settings to infer whether or not a choice of summary statistics is appropriate, via a Monte Carlo validation.

preprint2013arXiv

Two discussions of the paper "Bayesian measures of model complexity and fit" by D. Spiegelhalter et al., Read before The Royal Statistical Society at a meeting organized by the Research Section on Wednesday, March 13th, 2002

These are the written discussions of the paper "Bayesian measures of model complexity and fit" by D. Spiegelhalter et al. (2002), following the discussions given at the Annual Meeting of the Royal Statistical Society in Newcastle-upon-Tyne on September 3rd, 2013.

preprint2012arXiv

Bayesian computation via empirical likelihood

Approximate Bayesian computation (ABC) has become an essential tool for the analysis of complex stochastic models when the likelihood function is numerically unavailable. However, the well-established statistical method of empirical likelihood provides another route to such settings that bypasses simulations from the model and the choices of the ABC parameters (summary statistics, distance, tolerance), while being convergent in the number of observations. Furthermore, bypassing model simulations may lead to significant time savings in complex models, for instance those found in population genetics. The BCel algorithm we develop in this paper also provides an evaluation of its own performance through an associated effective sample size. The method is illustrated using several examples, including estimation of standard distributions, time series, and population genetics models.

preprint2007arXiv

Convergence of adaptive mixtures of importance sampling schemes

In the design of efficient simulation algorithms, one is often beset with a poor choice of proposal distributions. Although the performance of a given simulation kernel can clarify a posteriori how adequate this kernel is for the problem at hand, a permanent on-line modification of kernels causes concerns about the validity of the resulting algorithm. While the issue is most often intractable for MCMC algorithms, the equivalent version for importance sampling algorithms can be validated quite precisely. We derive sufficient convergence conditions for adaptive mixtures of population Monte Carlo algorithms and show that Rao--Blackwellized versions asymptotically achieve an optimum in terms of a Kullback divergence criterion, while more rudimentary versions do not benefit from repeated updating.