Strategic Monte Carlo Methods for State and Parameter Estimation in High Dimensional Nonlinear Problems
share
In statistical data assimilation one seeks the largest maximum of the conditional probability distribution P(X,p|Y) of model states, X, and parameters,p, conditioned on observations Y through minimizing the `action', A(X) = - P(X,p|Y). This determines the dominant contribution to the expected values of functions of X but does not give information about the structure of P(X,p|Y) away from the maximum. We introduce a Monte Carlo sampling method, called Strategic Monte Carlo (SMC) sampling, for estimating P(X, p|Y) in the neighborhood of its largest maximum to remedy this limitation. SMC begins with a systematic variational annealing (VA) procedure for finding the smallest minimum of A(X). SMC generates accurate estimates for the mean, standard deviation and other higher moments of P(X,p|Y). Additionally, the random search allows for an understanding of any multimodal structure that may underly the dynamics of the problem. SMC generates a gaussian probability control term based on the paths determined by VA to minimize a cost function A(X,p). This probability is sampled during the Monte Carlo search of the cost function to constrain the search to high probability regions of the surface thus substantially reducing the time necessary to sufficiently explore the space.
READ FULL TEXT