Adaptive algorithms for Stratified Sampling Monte Carlo

Add to your list(s) Download to your calendar using vCal

• Alexandra Carpentier (Statistical Laboratory)
• Wednesday 13 February 2013, 14:00-15:30
• MR9, CMS.

If you have a question about this talk, please contact Robin Evans.

Note the earlier time and different room to allow attendance at the RSS Ordinary Meeting

We consider the problem of estimating the integral of a function f over a domain. Although no analytic expression for f is available, it is possible to obtain n samples from f, chosen anywhere in the domain. A popular method for computing the integral of the function is to stratify the space in strata and sample points in the strata.

We propose an algorithm for returning a stratified estimate of the integral. We prove that this algorithm adapts online the number of samples in each stratum to the amount of variation of the function in the stratum. In particular, this enables to allocate more samples where the function varies more, and be almost as efficient as an “oracle” strategy that has access to the variations of the functions in each stratum. More precisions on this aspect is in paper (Carpentier and Munos, 2011). We also provide some results on (i) how to choose the number of strata in an efficient way and (ii) how to adapt the strata themselves to the specific shape of the function. We express those results with finite-time bounds on a proxy of the variance of the estimate (returned by the algorithms we present).

This talk is part of the Statistics Reading Group series.

This talk is included in these lists:

• All CMS events
• All Talks (aka the CURE list)
• CMS Events
• DPMMS Lists
• DPMMS info aggregator
• MR9, CMS
• School of Physical Sciences
• Statistical Laboratory info aggregator
• Statistics Group
• Statistics Reading Group

Note that ex-directory lists are not shown.