Realtime Sequential Inference of Static Parameters with Expensive Likelihood Calculations

Abstract

Summary A methodology is developed for making inference about parameters of a possible covert chemical or biological atmospheric release from sensor readings. The key difficulty in performing this inference is that the results must be obtained in a very short timescale (5 min) to make use of the inference for protection. The methodology that is developed uses some of the components in a sequential Monte Carlo algorithm. However, this inference problem is different from many other sequential Monte Carlo problems, in that there are no state evolution equations, the forward model is highly non-linear and the likelihoods are non-Gaussian. The algorithm that is developed can use stored output from complex physics models for more rapid update of the posterior from new data without having to rerun the models. The use of differential evolution Markov chain sampling allows new samples to diverge rapidly from degenerate sample sets. Results for inferences made of atmospheric releases (both real and simulated) of material are presented, demonstrating that the sampling scheme performs adequately despite constraints of a short time span for calculations.