Simulation based sequential Monte Carlo methods for discretely observed Markov processes

Parameter estimation for discretely observed Markov processes is a challenging problem. However, simulation of Markov processes is straightforward using the Gillespie algorithm. We exploit this ease of simulation to develop an effective sequential Monte Carlo (SMC) algorithm for obtaining samples from the posterior distribution of the parameters. In particular, we introduce two key innovations, coupled simulations, which allow us to study multiple parameter values on the basis of a single simulation, and a simple, yet effective, importance sampling scheme for steering simulations towards the observed data. These innovations substantially improve the efficiency of the SMC algorithm with minimal effect on the speed of the simulation process. The SMC algorithm is successfully applied to two examples, a Lotka-Volterra model and a Repressilator model.Comment: 27 pages, 5 figure

Optimal scaling for partially updating MCMC algorithms

In this paper we shall consider optimal scaling problems for high-dimensional Metropolis--Hastings algorithms where updates can be chosen to be lower dimensional than the target density itself. We find that the optimal scaling rule for the Metropolis algorithm, which tunes the overall algorithm acceptance rate to be 0.234, holds for the so-called Metropolis-within-Gibbs algorithm as well. Furthermore, the optimal efficiency obtainable is independent of the dimensionality of the update rule. This has important implications for the MCMC practitioner since high-dimensional updates are generally computationally more demanding, so that lower-dimensional updates are therefore to be preferred. Similar results with rather different conclusions are given for so-called Langevin updates. In this case, it is found that high-dimensional updates are frequently most efficient, even taking into account computing costs.Comment: Published at http://dx.doi.org/10.1214/105051605000000791 in the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org

Multitype randomized Reed--Frost epidemics and epidemics upon random graphs

We consider a multitype epidemic model which is a natural extension of the randomized Reed--Frost epidemic model. The main result is the derivation of an asymptotic Gaussian limit theorem for the final size of the epidemic. The method of proof is simpler, and more direct, than is used for similar results elsewhere in the epidemics literature. In particular, the results are specialized to epidemics upon extensions of the Bernoulli random graph.Comment: Published at http://dx.doi.org/10.1214/105051606000000123 in the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org

The basic reproduction number, R0R_0, in structured populations

In this paper, we provide a straightforward approach to defining and deriving the key epidemiological quantity, the basic reproduction number, $R_0$, for Markovian epidemics in structured populations. The methodology derived is applicable to, and demonstrated on, both $SIR$ and $SIS$ epidemics and allows for population as well as epidemic dynamics. The approach taken is to consider the epidemic process as a multitype process by identifying and classifying the different types of infectious units along with the infections from, and the transitions between, infectious units. For the household model, we show that our expression for $R_0$ agrees with earlier work despite the alternative nature of the construction of the mean reproductive matrix, and hence, the basic reproduction number.Comment: 26 page

Optimal scaling of the independence sampler : theory and practice

The independence sampler is one of the most commonly used MCMC algorithms usually as a component of a Metropolis-within-Gibbs algorithm. The common focus for the independence sampler is on the choice of proposal distribution to obtain an as high as possible acceptance rate. In this paper we have a somewhat different focus concentrating on the use of the independence sampler for updating augmented data in a Bayesian framework where a natural proposal distribution for the independence sampler exists. Thus we concentrate on the proportion of the augmented data to update to optimise the independence sampler. Generic guidelines for optimising the independence sampler are obtained for independent and identically distributed product densities mirroring findings for the random walk Metropolis algorithm. The generic guidelines are shown to be informative beyond the narrow confines of idealised product densities in two epidemic examples

Prediction of forces and moments for hypersonic flight vehicle control effectors

This research project includes three distinct phases. For completeness, all three phases of the work are briefly described in this report. The goal was to develop methods of predicting flight control forces and moments for hypersonic vehicles which could be used in a preliminary design environment. The first phase included a preliminary assessment of subsonic/supersonic panel methods and hypersonic local flow inclination methods for such predictions. While these findings clearly indicated the usefulness of such methods for conceptual design activities, deficiencies exist in some areas. Thus, a second phase of research was conducted in which a better understanding was sought for the reasons behind the successes and failures of the methods considered, particularly for the cases at hypersonic Mach numbers. This second phase involved using computational fluid dynamics methods to examine the flow fields in detail. Through these detailed predictions, the deficiencies in the simple surface inclination methods were determined. In the third phase of this work, an improvement to the surface inclination methods was developed. This used a novel method for including viscous effects by modifying the geometry to include the viscous/shock layer