Scalable Monte Carlo Methods
Arguably the main deterrent to more wide-spread use of Bayesian methods is their reliance on Monte Carlo methods, such as MCMC, which scale poorly to big-data settings, and are often unsuitable for implementation in a parallel computing environment. Approximate approaches, such as variational methods, do offer scalable alternatives. While these can perform well in terms of approximating the body of the posterior and making point predictions, they often give unreliable approximations to the tails. However, in many health science applications it is the tails of the posterior that are crucial in determining the best decisions.
The challenge with these methods stems from the extra difficulty in simulating a continuous-time process, and is around developing general implementations for important classes of statistical model. This will be an initial focus, particularly using recent insights we have developed in terms of using such algorithms in the place of reversible-jump MCMC for problems of model-choice within association studies. In some situations the challenges of simulating the continuous-time dynamics may lead to excessive computational overhead, and thus we will investigate theoretically and empirically the properties of algorithms that use approximate simulation ideas.