Bayesian Updating, Model Class Selection and Robust Stochastic Predictions of Structural Response

A fundamental issue when predicting structural response by using mathematical models is how to treat both modeling and excitation uncertainty. A general framework for this is presented which uses probability as a multi-valued conditional logic for quantitative plausible reasoning in the presence of uncertainty due to incomplete information. The fundamental probability models that represent the structure’s uncertain behavior are specified by the choice of a stochastic system model class: a set of input-output probability models for the structure and a prior probability distribution over this set that quantifies the relative plausibility of each model. A model class can be constructed from a parameterized deterministic structural model by stochastic embedding utilizing Jaynes’ Principle of Maximum Information Entropy. Robust predictive analyses use the entire model class with the probabilistic predictions of each model being weighted by its prior probability, or if structural response data is available, by its posterior probability from Bayes’ Theorem for the model class. Additional robustness to modeling uncertainty comes from combining the robust predictions of each model class in a set of competing candidates weighted by the prior or posterior probability of the model class, the latter being computed from Bayes’ Theorem. This higherlevel application of Bayes’ Theorem automatically applies a quantitative Ockham razor that penalizes the data-fit of more complex model classes that extract more information from the data. Robust predictive analyses involve integrals over highdimensional spaces that usually must be evaluated numerically. Published applications have used Laplace's method of asymptotic approximation or Markov Chain Monte Carlo algorithms

Approximate Bayesian Computation by Subset Simulation

A new Approximate Bayesian Computation (ABC) algorithm for Bayesian updating of model parameters is proposed in this paper, which combines the ABC principles with the technique of Subset Simulation for efficient rare-event simulation, first developed in S.K. Au and J.L. Beck [1]. It has been named ABC- SubSim. The idea is to choose the nested decreasing sequence of regions in Subset Simulation as the regions that correspond to increasingly closer approximations of the actual data vector in observation space. The efficiency of the algorithm is demonstrated in two examples that illustrate some of the challenges faced in real-world applications of ABC. We show that the proposed algorithm outperforms other recent sequential ABC algorithms in terms of computational efficiency while achieving the same, or better, measure of ac- curacy in the posterior distribution. We also show that ABC-SubSim readily provides an estimate of the evidence (marginal likelihood) for posterior model class assessment, as a by-product

Learning about a Moving Target in Resource Management: Optimal Bayesian Disease Control

Resource managers must often make difficult choices in the face of imperfectly observed and dynamically changing systems (e.g., livestock, fisheries, water, and invasive species). A rich set of techniques exists for identifying optimal choices when that uncertainty is assumed to be understood and irreducible. Standard optimization approaches, however, cannot address situations in which reducible uncertainty applies to either system behavior or environmental states. The adaptive management literature overcomes this limitation with tools for optimal learning, but has been limited to highly simplified models with state and action spaces that are discrete and small. We overcome this problem by using a recently developed extension of the Partially Observable Markov Decision Process (POMDP) framework to allow for learning about a continuous state. We illustrate this methodology by exploring optimal control of bovine tuberculosis in New Zealand cattle. Disease testing—the control variable—serves to identify herds for treatment and provides information on prevalence, which is both imperfectly observed and subject to change due to controllable and uncontrollable factors. We find substantial efficiency losses from both ignoring learning (standard stochastic optimization) and from simplifying system dynamics (to facilitate a typical, simple learning model), though the latter effect dominates in our setting. We also find that under an adaptive management approach, simplifying dynamics can lead to a belief trap in which information gathering ceases, beliefs become increasingly inaccurate, and losses abound

Potentials and Limits of Bayesian Networks to Deal with Uncertainty in the Assessment of Climate Change Adaptation Policies

Bayesian networks (BNs) have been increasingly applied to support management and decision-making processes under conditions of environmental variability and uncertainty, providing logical and holistic reasoning in complex systems since they succinctly and effectively translate causal assertions between variables into patterns of probabilistic dependence. Through a theoretical assessment of the features and the statistical rationale of BNs, and a review of specific applications to ecological modelling, natural resource management, and climate change policy issues, the present paper analyses the effectiveness of the BN model as a synthesis framework, which would allow the user to manage the uncertainty characterising the definition and implementation of climate change adaptation policies. The review will let emerge the potentials of the model to characterise, incorporate and communicate the uncertainty, with the aim to provide an efficient support to an informed and transparent decision making process. The possible drawbacks arising from the implementation of BNs are also analysed, providing potential solutions to overcome them.Adaptation to Climate Change, Bayesian Network, Uncertainty

Simulation methods for reliability-based design optimization and model updating of civil engineering structures and systems

This thesis presents a collection of original contributions pertaining to the subjects of reliability-based design optimization (RBDO) and model updating of civil engineering structures and systems. In this regard, probability theory concepts and tools are instrumental in the formulation of the herein reported developments. Firstly, two approaches are devised for the RBDO of structural dynamical systems under stochastic excitation. Namely, a stochastic search technique is proposed for constrained and unconstrained RBDO problems involving continuous, discrete and mixed discrete-continuous design spaces, whereas an efficient sensitivity assessment framework for linear stochastic structures is implemented to identify optimal designs and evaluate their sensitivities. Moreover, two classes of model updating problems are considered. In this context, the Bayesian interpretation of probability theory plays a key role in the proposed solution schemes. Specifically, contaminant source detection in water distribution networks is addressed by resorting to a sampling-based Bayesian model class selection framework. Furthermore, an effective strategy for Bayesian model updating with structural reliability methods is presented to treat identification problems involving structural dynamical systems, measured response data, and high-dimensional parameter spaces. The approaches proposed in this thesis integrate stochastic simulation techniques as an essential part of their formulation, which allows obtaining non-trivial information about the systems of interest as a byproduct of the solution processes. Overall, the findings presented in this thesis suggest that the reported methods can be potentially adopted as supportive tools for a number of practical decision-making processes in civil engineering.Diese Arbeit stellt eine Sammlung von Beiträgen vor, die sich mit der Reliability-based-Design-Optimization (RBDO) und dem Model updating von Strukturen und Systemen im Bauwesen befassen. In diesem Zusammenhang sind wahrscheinlichkeitstheoretische Konzepte für die Formulierung der hier vorgestellten Entwicklungen von entscheidender Bedeutung. Zunächst werden zwei Ansätze für eine RBDO von strukturdynamischen Systemen unter stochastischer Anregung entwickelt. Es wird eine stochastische Suchtechnik für beschränkte und unbeschränkte RBDO-Probleme vorgeschlagen. Diese beziehen kontinuierliche, diskrete und gemischt diskret-kontinuierliche Designräume ein. Gleichzeitig wird ein effizientes Framework zur Bewertung der Sensitivität lineare stochastische Strukturen implementiert, um optimale Designs zu identifizieren und ihre Sensitivitäten zu bewerten. Darüber hinaus werden zwei Klassen von Problem aus dem Model updating betrachtet. Der Fokus wird hierbei auf die Erkennung von Kontaminationsquellen in Wasserverteilungsnetzen mithilfe eines auf Stichproben basierenden Bayesian-Model-Class-selection-Framework gelegt. Ferner wird eine effektive Strategie zur Bearbeitung von Problemen des Bayesian-Model-updating, die strukturdynamischen Systeme, gemessene Systemantwortdaten und hochdimensionale Parameterräume umfassen, vorgestellt. Die beschriebenen Ansätze verwenden stochastische Simulationstechniken als wesentlicher Bestandteil ihrer Formulierung, wodurch nicht-triviale Informationen über betrachtete Systeme als Nebenprodukt der Lösungsprozesse gewonnen werden können. Insgesamt deuten die vorgestellten Ergebnisse dieser Arbeit darauf hin, dass die beschriebenen Methoden potenziell als unterstützende Elemente in praktischen Entscheidungsproblemen im Zusammenhang mit Strukturen und Systemen im Bauwesen eingesetzt werden können

Dual state-parameter estimation of hydrological models using ensemble Kalman filter

Hydrologic models are twofold: models for understanding physical processes and models for prediction. This study addresses the latter, which modelers use to predict, for example, streamflow at some future time given knowledge of the current state of the system and model parameters. In this respect, good estimates of the parameters and state variables are needed to enable the model to generate accurate forecasts. In this paper, a dual state-parameter estimation approach is presented based on the Ensemble Kalman Filter (EnKF) for sequential estimation of both parameters and state variables of a hydrologic model. A systematic approach for identification of the perturbation factors used for ensemble generation and for selection of ensemble size is discussed. The dual EnKF methodology introduces a number of novel features: (1) both model states and parameters can be estimated simultaneously; (2) the algorithm is recursive and therefore does not require storage of all past information, as is the case in the batch calibration procedures; and (3) the various sources of uncertainties can be properly addressed, including input, output, and parameter uncertainties. The applicability and usefulness of the dual EnKF approach for ensemble streamflow forecasting is demonstrated using a conceptual rainfall-runoff model. © 2004 Elsevier Ltd. All rights reserved

Bayesian sequential experimental design for binary response data with application to electromyographic experiments

We develop a sequential Monte Carlo approach for Bayesian analysis of the experimental design for binary response data. Our work is motivated by surface electromyographic (SEMG) experiments, which can be used to provide information about the functionality of subjects' motor units. These experiments involve a series of stimuli being applied to a motor unit, with whether or not the motor unit res for each stimulus being recorded. The aim is to learn about how the probability of ring depends on the applied stimulus (the so-called stimulus response curve); One such excitability parameter is an estimate of the stimulus level for which the motor unit has a 50% chance of ring. Within such an experiment we are able to choose the next stimulus level based on the past observations. We show how sequential Monte Carlo can be used to analyse such data in an online manner. We then use the current estimate of the posterior distribution in order to choose the next stimulus level. The aim is to select a stimulus level that mimimises the expected loss. We will apply this loss function to the estimates of target quantiles from the stimulus-response curve. Through simulation we show that this approach is more ecient than existing sequential design methods for choosing the stimulus values. If applied in practice, it could more than halve the length of SEMG experiments

The safety case and the lessons learned for the reliability and maintainability case

This paper examine the safety case and the lessons learned for the reliability and maintainability case