A Probabilistic Characterization of Simulation Model Uncertainties

Similar documents
Outlier-tolerant parameter estimation

Review - Probabilistic Classification

Analyzing Frequencies

CHAPTER 7d. DIFFERENTIATION AND INTEGRATION

Grand Canonical Ensemble

External Equivalent. EE 521 Analysis of Power Systems. Chen-Ching Liu, Boeing Distinguished Professor Washington State University

ST 524 NCSU - Fall 2008 One way Analysis of variance Variances not homogeneous

EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Lecture 12

Representation and Reasoning with Uncertain Temporal Relations

The Hyperelastic material is examined in this section.

10/7/14. Mixture Models. Comp 135 Introduction to Machine Learning and Data Mining. Maximum likelihood estimation. Mixture of Normals in 1D

A Note on Estimability in Linear Models

A NON-LINEAR MODEL FOR STUDYING THE MOTION OF A HUMAN BODY. Piteşti, , Romania 2 Department of Automotive, University of Piteşti

PHA 5128 Answer CASE STUDY 3 Question #1: Model

COMPLEX NUMBER PAIRWISE COMPARISON AND COMPLEX NUMBER AHP

UNIT 8 TWO-WAY ANOVA WITH m OBSERVATIONS PER CELL

Soft k-means Clustering. Comp 135 Machine Learning Computer Science Tufts University. Mixture Models. Mixture of Normals in 1D

te Finance (4th Edition), July 2017.

??? Dynamic Causal Modelling for M/EEG. Electroencephalography (EEG) Dynamic Causal Modelling. M/EEG analysis at sensor level. time.

Physics of Very High Frequency (VHF) Capacitively Coupled Plasma Discharges

Econ107 Applied Econometrics Topic 10: Dummy Dependent Variable (Studenmund, Chapter 13)

4D SIMPLICIAL QUANTUM GRAVITY

Economics 600: August, 2007 Dynamic Part: Problem Set 5. Problems on Differential Equations and Continuous Time Optimization

Outline. Types of Experimental Designs. Terminology. EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Lecture 12

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

Fakultät III Univ.-Prof. Dr. Jan Franke-Viebach

Stability and convergence of the Crank-Nicolson scheme for a class of variable-coefficient tempered fractional diffusion equations

Epistemic Foundations of Game Theory. Lecture 1

Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari

Unbalanced Panel Data Models

A general N-dimensional vector consists of N values. They can be arranged as a column or a row and can be real or complex.

2008 AP Calculus BC Multiple Choice Exam

Fault prediction using the combination of regularized OS-ELM and strong tracking SCKF

Physics 256: Lecture 2. Physics

3.4 Properties of the Stress Tensor

Chapter 6 Student Lecture Notes 6-1

Logistic Regression I. HRP 261 2/10/ am

Introduction to logistic regression

Today s logistic regression topics. Lecture 15: Effect modification, and confounding in logistic regression. Variables. Example

NATURAL VIBRATION ANALYSIS OF CONTINUOUS HORIZONTALLY CURVED GIRDER BRIDGES USING DQEM *

Small Area Estimation for Crop Acreage in Remote Sensing Assisted Crop Survey

Application of Local Influence Diagnostics to the Linear Logistic Regression Models

SCITECH Volume 5, Issue 1 RESEARCH ORGANISATION November 17, 2015

ENGR 323 BHW 15 Van Bonn 1/7

Jones vector & matrices

Security of Reconciliation in Continuous Variable Quantum Key Distribution *

Phys 2310 Fri. Nov. 7, 2014 Today s Topics. Begin Chapter 15: The Superposition of Waves Reading for Next Time

Maxwellian Collisions

Stat 543 Exam 2 Spring 2016

Math 102. Rumbos Spring Solutions to Assignment #8. Solution: The matrix, A, corresponding to the system in (1) is

Preference and Demand Examples

Lecture 14. Relic neutrinos Temperature at neutrino decoupling and today Effective degeneracy factor Neutrino mass limits Saha equation

8-node quadrilateral element. Numerical integration

4.8 Huffman Codes. Wordle. Encoding Text. Encoding Text. Prefix Codes. Encoding Text

What is LP? LP is an optimization technique that allocates limited resources among competing activities in the best possible manner.

Stat 543 Exam 2 Spring 2016

A Recursive Algorithm Based on Fuzzy Neural Network for Target Number

EXST Regression Techniques Page 1

Fakultät III Wirtschaftswissenschaften Univ.-Prof. Dr. Jan Franke-Viebach

Estimating the Variance in a Simulation Study of Balanced Two Stage Predictors of Realized Random Cluster Means Ed Stanek

Stress-Based Finite Element Methods for Dynamics Analysis of Euler-Bernoulli Beams with Various Boundary Conditions

Comparisons of the Variance of Predictors with PPS sampling (update of c04ed26.doc) Ed Stanek

Calculus II (MAC )

Exercises for lectures 7 Steady state, tracking and disturbance rejection

Linear Algebra. Definition The inverse of an n by n matrix A is an n by n matrix B where, Properties of Matrix Inverse. Minors and cofactors

Problem Set 4 Solutions Distributed: February 26, 2016 Due: March 4, 2016

sin VALUE LIMIT IN DOUBLE CHOOZ EXPERIMENT Pi-Jung Chang

2. Grundlegende Verfahren zur Übertragung digitaler Signale (Zusammenfassung) Informationstechnik Universität Ulm

Convergence Theorems for Two Iterative Methods. A stationary iterative method for solving the linear system: (1.1)

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional

(A) the function is an eigenfunction with eigenvalue Physical Chemistry (I) First Quiz

Spectral stochastic finite element analysis of structures with random field parameters under bounded-but-uncertain forces

Optimal Ordering Policy in a Two-Level Supply Chain with Budget Constraint

CHAPTER 33: PARTICLE PHYSICS

CSE-571 Probabilistic Robotics

Differentiation of Exponential Functions

Hostel Occupancy Survey (YHOS) Methodology

Stat 642, Lecture notes for 01/27/ d i = 1 t. n i t nj. n j

Lecture Slides for. ETHEM ALPAYDIN The MIT Press,

Consider a system of 2 simultaneous first order linear equations

On Selection of Best Sensitive Logistic Estimator in the Presence of Collinearity

Where k is either given or determined from the data and c is an arbitrary constant.

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

FREE VIBRATION ANALYSIS OF FUNCTIONALLY GRADED BEAMS

APPLICABILITY OF LINEARIZED DUSTY GAS MODEL FOR MULTICOMPONENT DIFFUSION OF GAS MIXTURES IN POROUS SOLIDS. Jelena Markovi and Radovan Omorjan

You already learned about dummies as independent variables. But. what do you do if the dependent variable is a dummy?

Basic Electrical Engineering for Welding [ ] --- Introduction ---

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Bayesian Estimation of the Logormal Distrbution Mean Using Ranked SET Sampling

Determination of Photon Interaction Parameters of CaO and MgO for Multi-Energetic Photons using γ-ray Attenuation Technique

Phy213: General Physics III 4/10/2008 Chapter 22 Worksheet 1. d = 0.1 m

BAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS. Dariusz Biskup

1997 AP Calculus AB: Section I, Part A

Code Design for the Low SNR Noncoherent MIMO Block Rayleigh Fading Channel

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

10. Limits involving infinity

Radial Cataphoresis in Hg-Ar Fluorescent Lamp Discharges at High Power Density

Fermi-Dirac statistics

From Structural Analysis to FEM. Dhiman Basu

Transcription:

A Proalstc Charactrzaton of Sulaton Modl Uncrtants Vctor Ontvros Mohaad Modarrs Cntr for Rsk and Rlalty Unvrsty of Maryland 1

Introducton Thr s uncrtanty n odl prdctons as wll as uncrtanty n xprnts Th uncrtanty n th xprnt and uncrtanty n th odl prdctons ar consdrd ndpndnt Any coparson twn th rsults of xprnt and odl prdctons ay usd to stat uncrtanty of th ral valu of ntrst Th odl uncrtanty s th only uncrtanty that should consdrd whn th dstruton of ral valu gvn a odl prdcton s statd. Thr dffrnt odl rrors ar prsntd: th addtv rror prcntag rror odl and ultplcatv rror odl Th ultplcatv rror odl shows a ttr agrnt

Addtv Error Modl: Assuptons Th addtv rror of th rsults of xprnt copard to th ral valus of ntrst s a norally dstrutd rando nur wth gvn an and standard dvaton Th addtv rror of odl prdctons copard to th ral valus of ntrst s also a norally dstrutd rando nur Error twn odl and xprnt can now assssd 3

Addtv Error Modl = + E ; E ~ N s Modl Prdcton Rsult of Exprnt = + E ; E ~ N s whr : : Ral Quantty : Rsult of xprnt : Modl prdcton : Man and SD of xprntal addtv rror : Man and SD of odl addtv rror E ~ s + s N 4

Ptfalls Error plottd twn odl prdcton and xprntal asurnts ust dntcally dstrutd whch s not always tru as thy ay ncras at hghr rangs. Th ntroducd addtv rror can ngatv zro or postv Th E whch s th addtv rror twn odl prdcton and xprnt can ngatv postv or zro. Ths lts th chocs of lklhood functon for ths rando varal to noral dstruton. Whn th data s wdly scattrd th noral dstruton assupton rsults n ngatv lowr ounds wth no anngful physcal ntrprtaton 5

Prcntag Error Modl: Assuptons Th prcntag rror of th rsults of xprnt copard to th ral valus of ntrst s a norally dstrutd rando nur wth gvn an and standard dvaton Th prcntag rror of odl prdctons copard to th ral valus of ntrst s a norally dstrutd rando nur Th prcntag rror of odl prdctons copard to th rsults of xprnt s a functon of th two rando varals ntroducd arlr. Th dstruton of ths rando nur wll usd to rprsnt th lklhood of data 6

Prcntag Error Modl Modl Prdcton Rsult of Exprnt whr : : Ral Quantty : Rsult of xprnt :Modl prdcton :Man and SD of xprntal prcntag rror : Man and SD of odl prcntag rror = E = E = E ; ; E E ~ N ~ N 1 E E E = E + 1 : Man and SD prcntag rrorof Exprnt copard to odl Indpndncy of E E } ~ N 1+ 1 + + 1 + 1 + 1 + + 1 + E 3 4 Approxaton! 7

Ptfalls Th ntroducd prcntag rror can ngatv zro or postv. Ths ascally forcs th noral dstruton assupton for prcntag rrors E & E. Th E whch s th prcntag rror twn odl prdcton and xprnt can ngatv postv or zro. Ths lts th chocs of lklhood functon for ths rando varal to noral dstruton. Th xact dstruton of E can not analytcally drvd Whn th data s wdly scattrd th noral dstruton assupton rsults n ngatv lowr ounds wth no anngful physcal ntrprtaton 8

Multplcatv Error: Assuptons Th odl prdcton rsult of xprnt and ral valu of ntrst hav th sa sgn all postv or all ngatv Th rato of ral valu and xprntal rsults s a rando varal wth lognoral dstruton for whch confdnc ounds ar known Th rato of ral valu and odl prdcton s a rando varal wth lognoral dstruton wth paratrs to dtrnd Th rato of odl prdctons and rsults of xprnt s a functon of th two rando varals ntroducd arlr. Th dstruton of ths rando varal s lognoral and wll usd to rprsnt th lklhood of data Havng th aov assuptons th dstruton of ral quantty of ntrst gvn a odl prdcton wll a lognoral dstruton 9

Multplcatv Error Modl Modl Prdcton Rsult of Exprnt = = = t Indpndncy of t ~ whr : : Ral Quantty : Rsult of xprnt : Modl prdcton :Th rror factor for odl prdctons : Man and SD of xprntal rror factor :Th rror factor for xprntal data : Man and SD of odl rror factor LN = = ; ; + ~ LN ~ LN 10

11 Multplcatv Error: Baysan Postror n LN LN f f L whr d d L f L f f π ln ~ ~ prdcton odl gvn as followng : as statd wll of dstruton th as such prdcton odl Gvn a Paratrs Postror Jont Dstruton of : Paratrs Pror Jont Dstruton of : 1 : 3 0 ln 1 1 0 0 + = + = = + =

Exapl IVE Radant Hat lux Baysan Approach 1

Proalty of Excdanc: HGL Tpratur 1 0.9 0.8 0.7 Excdanc Proalty 0.6 0.5 0.4 0.3 0. 0.1 Excdanc fro 18C Excdanc fro 330C 18 330 0 150 00 50 300 350 400 450 500 Prdctd HGL Pak Tpratur ºC 13

Conclusons Th rror statd y coparng xprnts and odl prdctons s not th uncrtanty of th odl. It s rathr a conaton of xprntal and odl uncrtants. Th dstruton of th ral valu of ntrst gvn a odl prdcton dpnds only on th uncrtanty of th odl. Th Baysan frawork allows dffrnt wghts and xprt judgnts to latr consdrd whn dalng wth non-hoognous populaton of xprnts or odl prdctons Postrors fro Baysan analyss can usd as pror to updatd y nw data ponts whn co avalal Multplcatv rror odl provds good rsults 14