14:00 - 15:00
Jackie Wong Siaw Tze
Lectures, talks and seminars
Mathematical Sciences, Department of
Osama Mahmoud email@example.com
These Departmental Seminars are for everyone in Maths. We encourage anyone interested in the subject in general, or in the particular subject of the seminar, to come along. It's a great opportunity to meet people in the Maths Department and join in with our community.
Computation of normalising constants is a fundamental mathematical problem in various disciplines, particularly in Bayesian model selection problems. A sampling-based technique known as bridge sampling (Meng and Wong, 1996) has been found to produce accurate estimates of normalising constants and is shown to possess good asymptotic properties.
For small to moderate sample sizes (as in situations with limited computational resources), we demonstrate that the (optimal) bridge sampler produces biased estimates. Specifically, when one density (we denote as p2) is constructed to be close to the target density (we denote as p1) using method of moments, our simulation-based results indicate that the correlation-induced bias through the moment-matching procedure is on-negligible.
More crucially, the bias amplifies as the dimensionality of the problem increases. Thus, a series of theoretical as well as empirical investigations is carried out to identify the nature and origin of the bias.
We then examine the effect of sample size allocation on the accuracy of bridge sampling estimates and discovered that one possibility of reducing both the bias and standard error with a small increase in computational effort is by drawing extra samples from the moment-matched density p2 (which we assume easy to sample from), provided that the evaluation of p1 is not too expensive. We proceed to show how the simple adaptive approach we termed “splitting” manages to alleviate the correlation-induced bias at the expense of a higher standard error, irrespective of the dimensionality involved. We also slightly modified the strategy suggested by Wang et al. (2019) to address the issue of the increase in standard error due to splitting, which is later generalised to further improve the efficiency.
We conclude the paper by offering our insights of the application of a combination of these adaptive methods to improve the accuracy of bridge sampling estimates in Bayesian applications (where posterior samples are typically expensive to generate) based on the preceding investigations, with an application to a practical example.
Jackie Wong Sia Tze, University of Essex
If not a member of the Dept. Mathematical Science at the University of Essex, you can register your interest in attending the seminar and request the Zoom’s meeting password by emailing Dr Osama Mahmoud (firstname.lastname@example.org)