Uncertainties in modelling
|
|
- Oswin Armstrong
- 6 years ago
- Views:
Transcription
1 Uncertainties in modelling Luca Vezzaro Modelling and Control of Environmental Systems Padova, 14 th January 2015
2 Introduction: Where are uncertainties? Why worrying about that? Theoretical classification of uncertainty How can be described? Uncertainty Assessment How do we consider uncertainty in modelling? Small examples How can we use uncertainty? 2
3 Introduction: Where are uncertainties? Why worrying about that? Theoretical classification of uncertainty How can be described? Uncertainty Assessment How do we consider uncertainty in modelling? Small examples How can we use uncertainty? 3
4 12 October 2014 The warning was not sent, because the model did not predict the flood model that had never been wrong in the past Wrong (?) model prediction
5 13-14 October 2012 High resolution hours in advance!! From 14 October
6 How weather prediction models work From BBC Two The Code Initial state Ensamble 1 day 3 days 1 week 6
7 Uncertainty affect all mathematical models INPUT SUBMODEL OUTPUT Weather forecast Complexity INPUT SUBMODEL OUTPUT Force equation INPUT MODEL OUTPUT F=m. a 7
8 Let s take to the carnival What time should we leave Padua? Train 8
9 Padova-Venice transport model X t = x /v What about stops on the way? t = x /v + t stop * n stops Fixed Fixed Can vary 9
10 Padova-Venice transport model t = x /v + t stop * n stops X = 40 km n stop = 6 v =? t stop =? What are the values for the model parameters? t obs,1 = 50 min t obs,2 = 55 min 10
11 Let s guess the parameters values X = 37 km n stop = 6 T obs,1 = 50 min T obs,2 = 55 min v [km/h] t stop [min] t [min]
12 On which parameter should we spend time on? Which one are the really important ones? Is there something more important to look at? 12
13 How models are applied System under study MEASUREMENTS Experimenting Optimized System Reality WICKED PROBLEM Virtual Reality PARAMETERS Model of the System Simulating NUMERICAL ERRORS Solution for the System 13
14 Why is important to assess uncertainty? Uncertainty can be used as a wild card Can be used to argue for precaution (or for not acting at all) Pro / Against Global Warming Or for justifying actions ( like invading a country) "As we know, there are known knowns. There are things we know we know. We also know there are known unknowns. That is to say we know there are some things we do not know. But there are also unknown unknowns, the ones we don't know we don't know." Donald Rumsfeld, former U.S. Secretary of Defense 14
15 Can we trust models? Yes, but we need to address uncertainty!! All models are wrong, but some are useful (Box&Draper, 1987) 15
16 Introduction: Where are uncertainties? Why worrying about that? Theoretical classification of uncertainty How can be described? Uncertainty Assessment How do we consider uncertainty in modelling? Small examples How can we use uncertainty? 16
17 What is uncertainty? What is often referred to as "uncertainty" actually hides important technical distinctions (EEA, 2001); Any departure from the unachievable ideal of complete determinism A three dimensional concept (e.g. Walker et al., 2003) Location Nature Level 17
18 Location of Uncertainty (1 st dimension) Generic Locations: Context (e.g. boundaries, framing) Model Structure Input data Calibration data Parameter Uncertainty Model Output (conclusion) Location Level Nature 18
19 Location of uncertainty Context vs. Model Uncertainty Context Uncertainty Wrong definition of the system: the wrong question was answered! 19
20 Let s take to the carnival Train Train The carnival is in Verona 20
21 Wetland modelling You develop a model that is very good at predicting oxygen and ammonia during rain events but you are interested in nutrients loads on an yearly basis 21
22 Location of uncertainty Context vs. Model Uncertainty Input data Model Structure is the description of the internal relationships dominating the system Model structure uncertainty there are competing models describing the system There is not a single way to model your system 22
23 Let s take to the carnival Bus Regional Train High speed train There are many ways to get to Venice 23
24 Uncertainty in environmental models Location The location of uncertainty can be: Model structure Relation between the output of the model, the input u and state variable x e.g. Removal for stormwater pollutants dl dt K a K rem C w P m L dl dt K a K rem C w V L Proportional to intensity or runoff volume? 24
25 Wetland modelling Conceptual models (tanks in series)? hydrodynamic models (Saint- Venant)? Stochastic models???? 25
26 Uncertainty in environmental models Location Souce: Hølten- Lützhøf et al 2008 Model parameters Parameter values can vary with wide ranges e.g. Sorption capacity for benzene (soil/water coefficient k d ) 81 L/kg 85 L/kg 370 L/kg 320 L/kg 280 L/kg 100 L/kg 92 L/kg 59 L/kg 38 L/kg 54 L/kg 44 L/kg 38 L/kg 31 L/kg 79 L/kg 50 L/kg 4 L/kg 31 L/kg 614 L/kg 337 L/kg 95 L/kg 89 L/kg 83 L/kg 85 L/kg 54 L/kg 26 L/kg 13 L/kg 22 L/kg 26
27 Location of uncertainty Data and parameter uncertainty Measurements of global sea temperature Source: Zecca and Chiari, 2009 No model is able to simulate this peak During WWII measurements were biased 27
28 Location of uncertainty Data and parameter uncertainty Rainfall-runoff model Delay in the model s response Simulated but not measured Measured but not simulated 28
29 Location of uncertainty Data and parameter uncertainty Modelled catchment Rain Gauge Input uncertainty Source: Goole Earth,
30 Level of Uncertainty (2 nd dimension) Statistical Uncertainty Scenario Uncertainty Recognised Ignorance Total Ignorance Location Level Nature 30
31 Statistical Uncertainty There exist solid grounds for the assignment of a discrete probability to each of a well-defined set of outcomes. We have a well known functional relationship We have an adequate combination of: (i) number of parameters and (ii) amount and character of data σ=0 μ=
32 Scenario Uncertainty We can describe a set of outcomes to be expected, but we cannot associate probabilities very well. Assumptions; various plausible scenarios; unverified what if? questions; Ambiguous results ???
33 Statistisc and scenario uncertainty Variable, X Scenario 2 Scenario 1 Confidence level Confidence interval Scenario 3 Estimation Planning horizon Time Scenarios don t need to be related to time (e.g. other variables) People usually don t talk about uncertainties in the scenarios. 33
34 Ignorance Statistical Uncertainty Scenario Uncertainty Recognised Ignorance Total Ignorance We do not know the essential functional relationships. There exist neither grounds for the assignment of probabilities, nor even the basis for the definition of a complete set of outcomes. More information may become known later through research, but little is known for the time being. Recognized ignorance: Total ignorance: We know that we don t know! We don t know what we don t know 34
35 Level of Uncertainty Statistical Uncertainty Scenario Uncertainty Recognised Ignorance Total Ignorance Uncertainty Uncertainty Uncertainty known known unknown Probability Probability Probability known unknown unknown That s were we are working so far Location Level Nature 35
36 The Nature of Uncertainty (3 rd dimension) Stochastic Uncertainty: variability and inherent indeterminism in the underlying natural phenomena. You cannot avoid it Epistemic Uncertainty: variability and inherent indeterminism in the underlying natural phenomena. Can be reduced, but.. Ambiguity: no matter how much energy you spend, you will find contrasting (but equally good) results Location Level Nature 36
37 Three dimensions of uncertainty Location of uncertainty Context, model structure, input data, parameters Level of uncertainty Statistical Scenario Recognized Total Uncertainty uncertainty ignorance ignorance Nature of uncertainty Can it be reduced or not? (epistemic or aleatoric) 37
38 How to classify uncertainty Location Level Nature From Warmink et al.,
39 Uncertainty matrix Try to identify the uncertainties in a model Context Level of uncertainty Nature Location of uncertainty Not Statistical Scenario Ignorance Reducible reducible Model structure Input data Parameters 39
40 Introduction: Where are uncertainties? Why worrying about that? Theoretical classification of uncertainty How can be described? Uncertainty Assessment How do we consider uncertainty in modelling? Small examples How can we use uncertainty? 40
41 Uncertainty in environmental models General formulation of a dynamic model m k : observed value time space input parameters But we are living in the real world (which is not perfect): O x, tox, t mk,, u, u, x, tm,, u, u, x, t r observations error parameters error input error model error random error 41
42 The big question which model factors are affecting my output? Rain Inputs System attributes y= f(x(t),θ) Model parameters To find what is relevant To quantify how uncertain Sensitivity analysis Uncertainty analysis 42
43 Addressing the model uncertainties Without observations to compare with Sensitivity Analysis (SA) Identification of sources of uncertainty With observation to compare with Uncertainty Analysis (UA) Quantification of uncertainty 43
44 Padova-Venice transport model v [km/h] t stop [min] t [min] S S v = 0.72 i y y i i i i v [km/h] t stop [min] t [min] S tstop = 0.11 v [km/h] t stop [min] t [min] S tstop =
45 Sensitivity analysis Once-At-Time method Traditionally, One-At-Time (OAT) methods are used: I change the parameter by Δθ I run the model and I obtain Δy I calculate the sensitivity index Not all the parameter space is explored!!! Δθ 2 S y Explored parameter space What about interactions? (θ 1,θ 2 ) Δθ 1 Actual parameter space 47
46 Global Sensitivity Analysis (GSA) Some examples Variance decomposition Identifies the direct and indirect influence σ 2 Output variance σ 2 Variance due to variation of parameter θ i σ 2 Variance due to variation of input Very good overview of the model behaviour Heavy computational requirements Can be applied also to models with many submodels? Variance due to model structural uncertainty 48
47 Sobol sensitivity indices First order index σ 2 Variance due to variation of parameter θ i σ 2 Output variance Total Sensitivity index 1- σ 2 σ 2 Variance due to variation of ALL parameters EXCLUDED θ i Output variance θ 1 θ 2 θ σ 2 θ 1 θ 2 θ σ 2 θ 1 θ 2 θ σ 2 49
48 Global Sensitivity Analysis (GSA) if thousand simulations are too many Morris method OAT method applied in several points of the parameter space Trajectory 2 Δθ 3 Δθ 1 Δθ 2 Source: Pujol, 2009 Many first order sensitivity indices Average μ Variance σ 2 θ 3 Information about: Trajectory 1 Sensitive factors Δθ 3 Δθ 1 Factors with interactions or non-linear effect θ 2 Δθ 2 Limited computational requirements θ 1 for k parameters: (k+1) * size of parameter sample (a couple of dozens) 50
49 Morris method Results analysis Information about: Sensitive factors Factors with interactions or non-linear effect 51
50 Regression approaches Run the model in different point of the parameter space Assume a linear response Estimate regression parametes α, β, γ The regression parameters are proportional to sensitivity indices 52
51 Example of GSA A stormwater quality model 9 model factors 4 parameters for the hydraulic submodel 3 parameters for the quality submodel 2 input error factors Influence on two outputs Mass discharged Concentration in the outlet Sample size = parameter sets 53
52 GSA results 1st order sensitivity indices (what you get from OAT) Deposition rate Washoff rate All the other factors are not relevant With OAT I would stop here 54
53 GSA results look at the interactions Two correlated parameters The hydraulic submodel affects the quality estimation 55
54 Sensitivity analysis step-by-step 1 Define prior distribution for each model parameter 2 Generate N parameter sets N 3 Run model for N parameter sets model out i 4 Evaluate performance of the N simulations out i RMSE 5 Analyze results 56
55 Identification of sources of uncertainty NO Measurements Global Sensitivity Analysis Which factors should I focus on? 757
56 Uncertainty Analysis How to estimate the term ε Various approaches are available: Bayesian and pseudo-bayesian approaches Grey Box modelling Etc. r u m u k O t x u t x u m t x t x O,,,,,,,,,,,, +/- XX%
57 Grey Box method The error is modelled as a stochastic term New parameters, defining the error structure are added to the model The uncertainty of the model is given by the simulation with the parameters which gives the best predictions ˆ 59
58 Grey Box method Parameter estimated through calibration Uncertainty bounds Point prediction 60
59 Statistical methods The Bayes theorem Conditional probability of B given A Prior probability of A Posterior probability Example from Wikipedia: Prior probability of B Suppose there is a school having 60% boys and 40% girls as students P(A)=0,4. The girl students wear trousers or skirts in equal numbers P(B A)=0,5; the boys all wear trousers P(B A )=1 An observer sees a (random) student from a distance; all they can see is that this student is wearing trousers. What is the probability this student is a girl? P(A B)=0,25 61
60 Statistical methods The Bayes theorem Translated in modelling terms Prior probability of parameters θ P ( y) L( y ) P( ) / K Posterior probability of parameters θ given the observation y Likelihood function It s a measure of the fitness of the model output compared to the results (e.g. based on square error, etc.) 62
61 Bayesian approaches P ( y) L( y ) P( ) / K L 2 L 1 L 3 1 Estimate likelihood 2 Estimate posterior distribution 3 Use posterior distribution to obtain uncertainty bounds P(θ/y) θ 63
62 Statistical methods The Generalized Likelihood Uncertainty Estimation (GLUE) Soft version of the Bayesian approach (pseudo-bayesian) Weight the models according to the goodness of their prediction Based on the equifinality principle: Observations can be simulated by many different parameter sets with similar goodness 64
63 Pseudo-Bayesian approaches (GLUE) L 3 L 2 L 1 L 3 1 Estimate informal likelihood P(θ/y) 2 Find behavioural parameter sets 3 Use behavioural sets to obtain uncertainty bounds θ 65
64 Informal likelihoods How do I evaluate my model performance? N i T T i Z Y L M 2 1, N obs i T T i Z Y L M 2 2 1, N t Z t M Z Y L M N obs j j T j i T T i 1 2 exp, N obs j j T j i obs i t Z t M N All based on the Square Error: New ways of considering observation errors 66
65 Statistical methods Bayes vs. GLUE The two approaches are equally complex Bayes is more mathematically elegant and safe GLUE is more relaxed (less assumptions) but results are subjective Mathematicians are still fighting, so we still have to see who will win E.g. four methods similar uncertainty bounds From Dotto et al. (2012) 67
66 Comparison of the approaches How do I choose my method? An example of comparison of different approaches Flow module TSS Module Coverage of GLUE SCEM-UA AMALGAM BAYES observations Computational requirements Coverage of Ability to identify ill-posed models (correlated parameters, poor efficiency and low coverage) Ability to identify sensitive observations Computational requirements parameters Availability Software package Code for use in MATLAB and OCTAVE Software package Software package Required programming skills Min Min Max Max Limitations Subjectivity of acceptance threshold Subjectivity of acceptance threshold Subjectivity of acceptance threshold Knowledge about error distribution From Dotto et al. (2012) 68
67 Simple receipt for UA 1. 1 Choose the parameters that are uncertain 2. 2 Choose a distribution for these parameters (e.g. uniform, triangular, gaussian, etc.) 3. 3 Generate a sample of parameters 4. Choose a (un)formal likelihood measure evaluating the 4 parameters 5. 5 Run the model with the generate sample Bayes GLUE Estimate the parameter posterior distribution Select the parameter sets with better performances 69
68 How the results of UA look like? Confidence bounds Bayes: confidence bounds (e.g. 95% probability that your results lie within the bands) GLUE: prediction bounds (e.g. the best parameters give prediction within the interval) 70
69 How the results of UA look like? Information about parameter intervals Flat distribution = non sensitive Correlation between parameters Peak distribution = sensitive Parameter intervals (you can use afterwards for predictions) 71
70 Introduction: Where are uncertainties? Why worrying about that? Theoretical classification of uncertainty How can be described? Uncertainty assessment How do we consider uncertainty in modelling? Small examples How can we use uncertainty? 72
71 Uncertainty in decision making Everyday pic-nic planning in Denmark Next 48 hours days 73
72 Uncertainty in decision making Everyday pic-nic planning in Denmark Next 48 hours Highly uncertain rain High probability rain days 74
73 Model for stormwater treatment Identification of most important parameters Multi-CSTR model (can simulate different units) Over 12 parameters (related and not-related to PP) What are the major sources of uncertainty? Analysis performed on: 2 Systems (different removal) Pond 6 Substance (different environmental fate) Heavy metals (Cu,Zn) sorbing Pyrene Biofilter Glyphosate Benzene IPBC biodegradable volatile no dominant processes 75
74 Model for stormwater treatment Identification of most important parameters Pond Most important factors are related to TSS (settling/resuspension) Biofilter Majority of important factors are related to hydraulic submodel (evapotranspiration) Do not need to spend a lot of money in measuring micropollutants! 27
75 Discharged MP loads Relative reduction Copper Zinc Fluoranthene 77
76 Compliance with Water Quality criteria EMC ELV is exceeded between XX and YY times per year Return period 78
77 Uncertainty in control of urban drainage networks Didactical example Detention basins Treatment plant
78 Example Real Time Control West Town Objective: Maximize storage East Town
79 Example Model Predictive Control West Town Objective: Maximize future storage East Town Model forecast (without uncertainty)
80 Example Risk-based Model Predictive Control? Rainfall evolution is uncertain
81 Example Risk-based Model Predictive Control Risk of overflow West Town East Town Target Target If we do not consider uncertainty If we consider uncertainty Objective: Minimize CSO risk
82 Conclusion Models are an essential tool for managing environmental systems But they are all wrong!!! (but some are useful) Acknowledging the models limits, failure and uncertainties is essential Managing model uncertainties is a complex issue, which a lot of people struggle to understand When you present your model result don t give a number, give an interval! 84
83 Thank you for listening!!! For any question:
SYDE 372 Introduction to Pattern Recognition. Probability Measures for Classification: Part I
SYDE 372 Introduction to Pattern Recognition Probability Measures for Classification: Part I Alexander Wong Department of Systems Design Engineering University of Waterloo Outline 1 2 3 4 Why use probability
More informationThe Generalized Likelihood Uncertainty Estimation methodology
CHAPTER 4 The Generalized Likelihood Uncertainty Estimation methodology Calibration and uncertainty estimation based upon a statistical framework is aimed at finding an optimal set of models, parameters
More informationTHE DEVELOPMENT OF RAIN-BASED URBAN FLOOD FORECASTING METHOD FOR RIVER MANAGEMENT PRACTICE USING X-MP RADAR OBSERVATION
Research Paper Advances in River Engineering, JSCE, Vol.19, 2013,June THE DEVELOPMENT OF RAIN-BASED URBAN FLOOD FORECASTING METHOD FOR RIVER MANAGEMENT PRACTICE USING X-MP RADAR OBSERVATION Seongsim YOON
More informationFlash Flood Guidance System On-going Enhancements
Flash Flood Guidance System On-going Enhancements Hydrologic Research Center, USA Technical Developer SAOFFG Steering Committee Meeting 1 10-12 July 2017 Jakarta, INDONESIA Theresa M. Modrick Hansen, PhD
More informationUncertainty propagation in a sequential model for flood forecasting
Predictions in Ungauged Basins: Promise and Progress (Proceedings of symposium S7 held during the Seventh IAHS Scientific Assembly at Foz do Iguaçu, Brazil, April 2005). IAHS Publ. 303, 2006. 177 Uncertainty
More informationNaïve Bayes classification
Naïve Bayes classification 1 Probability theory Random variable: a variable whose possible values are numerical outcomes of a random phenomenon. Examples: A person s height, the outcome of a coin toss
More informationCPSC 340: Machine Learning and Data Mining. MLE and MAP Fall 2017
CPSC 340: Machine Learning and Data Mining MLE and MAP Fall 2017 Assignment 3: Admin 1 late day to hand in tonight, 2 late days for Wednesday. Assignment 4: Due Friday of next week. Last Time: Multi-Class
More informationDevelopment of Stochastic Artificial Neural Networks for Hydrological Prediction
Development of Stochastic Artificial Neural Networks for Hydrological Prediction G. B. Kingston, M. F. Lambert and H. R. Maier Centre for Applied Modelling in Water Engineering, School of Civil and Environmental
More informationState-space calibration of radar rainfall and stochastic flow forecasting for use in real-time Control of urban drainage systems
State-space calibration of radar rainfall and stochastic flow forecasting for use in real-time Control of urban drainage systems 9th Int. Conf. On Urban Drainage Modeling Roland Löwe *, Peter Steen Mikkelsen,
More informationHarvard University. Rigorous Research in Engineering Education
Statistical Inference Kari Lock Harvard University Department of Statistics Rigorous Research in Engineering Education 12/3/09 Statistical Inference You have a sample and want to use the data collected
More informationNaïve Bayes classification. p ij 11/15/16. Probability theory. Probability theory. Probability theory. X P (X = x i )=1 i. Marginal Probability
Probability theory Naïve Bayes classification Random variable: a variable whose possible values are numerical outcomes of a random phenomenon. s: A person s height, the outcome of a coin toss Distinguish
More informationBayesian Learning (II)
Universität Potsdam Institut für Informatik Lehrstuhl Maschinelles Lernen Bayesian Learning (II) Niels Landwehr Overview Probabilities, expected values, variance Basic concepts of Bayesian learning MAP
More informationWhaddya know? Bayesian and Frequentist approaches to inverse problems
Whaddya know? Bayesian and Frequentist approaches to inverse problems Philip B. Stark Department of Statistics University of California, Berkeley Inverse Problems: Practical Applications and Advanced Analysis
More informationStatistical Methods in Particle Physics
Statistical Methods in Particle Physics Lecture 11 January 7, 2013 Silvia Masciocchi, GSI Darmstadt s.masciocchi@gsi.de Winter Semester 2012 / 13 Outline How to communicate the statistical uncertainty
More informationStat 535 C - Statistical Computing & Monte Carlo Methods. Arnaud Doucet.
Stat 535 C - Statistical Computing & Monte Carlo Methods Arnaud Doucet Email: arnaud@cs.ubc.ca 1 CS students: don t forget to re-register in CS-535D. Even if you just audit this course, please do register.
More informationLecture 3. STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher
Lecture 3 STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher Previous lectures What is machine learning? Objectives of machine learning Supervised and
More informationMachine Learning. Classification. Bayes Classifier. Representing data: Choosing hypothesis class. Learning: h:x a Y. Eric Xing
Machine Learning 10-701/15 701/15-781, 781, Spring 2008 Naïve Bayes Classifier Eric Xing Lecture 3, January 23, 2006 Reading: Chap. 4 CB and handouts Classification Representing data: Choosing hypothesis
More informationRadar-raingauge data combination techniques: A revision and analysis of their suitability for urban hydrology
9th International Conference on Urban Drainage Modelling Belgrade 2012 Radar-raingauge data combination techniques: A revision and analysis of their suitability for urban hydrology L. Wang, S. Ochoa, N.
More informationAchilles: Now I know how powerful computers are going to become!
A Sigmoid Dialogue By Anders Sandberg Achilles: Now I know how powerful computers are going to become! Tortoise: How? Achilles: I did curve fitting to Moore s law. I know you are going to object that technological
More informationUniversität Potsdam Institut für Informatik Lehrstuhl Maschinelles Lernen. Bayesian Learning. Tobias Scheffer, Niels Landwehr
Universität Potsdam Institut für Informatik Lehrstuhl Maschinelles Lernen Bayesian Learning Tobias Scheffer, Niels Landwehr Remember: Normal Distribution Distribution over x. Density function with parameters
More informationModern Methods of Data Analysis - SS 2009
Modern Methods of Data Analysis Lecture X (2.6.09) Contents: Frequentist vs. Bayesian approach Confidence Level Re: Bayes' Theorem (1) Conditional ( bedingte ) probability Examples: Rare Disease (1) Probability
More informationParameter estimation Conditional risk
Parameter estimation Conditional risk Formalizing the problem Specify random variables we care about e.g., Commute Time e.g., Heights of buildings in a city We might then pick a particular distribution
More informationModelling Under Risk and Uncertainty
Modelling Under Risk and Uncertainty An Introduction to Statistical, Phenomenological and Computational Methods Etienne de Rocquigny Ecole Centrale Paris, Universite Paris-Saclay, France WILEY A John Wiley
More informationMarkov Chain Monte Carlo
Department of Statistics The University of Auckland https://www.stat.auckland.ac.nz/~brewer/ Emphasis I will try to emphasise the underlying ideas of the methods. I will not be teaching specific software
More informationShould all Machine Learning be Bayesian? Should all Bayesian models be non-parametric?
Should all Machine Learning be Bayesian? Should all Bayesian models be non-parametric? Zoubin Ghahramani Department of Engineering University of Cambridge, UK zoubin@eng.cam.ac.uk http://learning.eng.cam.ac.uk/zoubin/
More informationCurrent best practice of uncertainty forecast for wind energy
Current best practice of uncertainty forecast for wind energy Dr. Matthias Lange Stochastic Methods for Management and Valuation of Energy Storage in the Future German Energy System 17 March 2016 Overview
More informationDay 5: Generative models, structured classification
Day 5: Generative models, structured classification Introduction to Machine Learning Summer School June 18, 2018 - June 29, 2018, Chicago Instructor: Suriya Gunasekar, TTI Chicago 22 June 2018 Linear regression
More informationIntroduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak
Introduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak 1 Introduction. Random variables During the course we are interested in reasoning about considered phenomenon. In other words,
More informationConfidence intervals CE 311S
CE 311S PREVIEW OF STATISTICS The first part of the class was about probability. P(H) = 0.5 P(T) = 0.5 HTTHHTTTTHHTHTHH If we know how a random process works, what will we see in the field? Preview of
More information2. I will provide two examples to contrast the two schools. Hopefully, in the end you will fall in love with the Bayesian approach.
Bayesian Statistics 1. In the new era of big data, machine learning and artificial intelligence, it is important for students to know the vocabulary of Bayesian statistics, which has been competing with
More informationOverview. Probabilistic Interpretation of Linear Regression Maximum Likelihood Estimation Bayesian Estimation MAP Estimation
Overview Probabilistic Interpretation of Linear Regression Maximum Likelihood Estimation Bayesian Estimation MAP Estimation Probabilistic Interpretation: Linear Regression Assume output y is generated
More informationRegional Flash Flood Guidance and Early Warning System
WMO Training for Trainers Workshop on Integrated approach to flash flood and flood risk management 24-28 October 2010 Kathmandu, Nepal Regional Flash Flood Guidance and Early Warning System Dr. W. E. Grabs
More informationFitting a Straight Line to Data
Fitting a Straight Line to Data Thanks for your patience. Finally we ll take a shot at real data! The data set in question is baryonic Tully-Fisher data from http://astroweb.cwru.edu/sparc/btfr Lelli2016a.mrt,
More information(1) Introduction to Bayesian statistics
Spring, 2018 A motivating example Student 1 will write down a number and then flip a coin If the flip is heads, they will honestly tell student 2 if the number is even or odd If the flip is tails, they
More informationBAYESIAN DECISION THEORY
Last updated: September 17, 2012 BAYESIAN DECISION THEORY Problems 2 The following problems from the textbook are relevant: 2.1 2.9, 2.11, 2.17 For this week, please at least solve Problem 2.3. We will
More informationUncertainty in energy system models
Uncertainty in energy system models Amy Wilson Durham University May 2015 Table of Contents 1 Model uncertainty 2 3 Example - generation investment 4 Conclusion Model uncertainty Contents 1 Model uncertainty
More informationUncertainties in the Projections of Future Heavy Rainfall: Climate Non-stationarity and Urban Flood Risk in Greater Chicago
Uncertainties in the Projections of Future Heavy Rainfall: Climate Non-stationarity and Urban Flood Risk in Greater Chicago Momcilo Markus and Gregory Byard Rainfall frequency sources TP-40, ISWS Bulletin
More informationClassification and Regression Trees
Classification and Regression Trees Ryan P Adams So far, we have primarily examined linear classifiers and regressors, and considered several different ways to train them When we ve found the linearity
More informationAnalysis of Regression and Bayesian Predictive Uncertainty Measures
Analysis of and Predictive Uncertainty Measures Dan Lu, Mary C. Hill, Ming Ye Florida State University, dl7f@fsu.edu, mye@fsu.edu, Tallahassee, FL, USA U.S. Geological Survey, mchill@usgs.gov, Boulder,
More informationCPSC 340: Machine Learning and Data Mining
CPSC 340: Machine Learning and Data Mining MLE and MAP Original version of these slides by Mark Schmidt, with modifications by Mike Gelbart. 1 Admin Assignment 4: Due tonight. Assignment 5: Will be released
More informationFeature-specific verification of ensemble forecasts
Feature-specific verification of ensemble forecasts www.cawcr.gov.au Beth Ebert CAWCR Weather & Environmental Prediction Group Uncertainty information in forecasting For high impact events, forecasters
More informationA Data-driven Approach for Remaining Useful Life Prediction of Critical Components
GT S3 : Sûreté, Surveillance, Supervision Meeting GdR Modélisation, Analyse et Conduite des Systèmes Dynamiques (MACS) January 28 th, 2014 A Data-driven Approach for Remaining Useful Life Prediction of
More informationMeasurement And Uncertainty
Measurement And Uncertainty Based on Guidelines for Evaluating and Expressing the Uncertainty of NIST Measurement Results, NIST Technical Note 1297, 1994 Edition PHYS 407 1 Measurement approximates or
More informationAt the start of the talk will be a trivia question. Be prepared to write your answer.
Operational hydrometeorological forecasting activities of the Australian Bureau of Meteorology Thomas Pagano At the start of the talk will be a trivia question. Be prepared to write your answer. http://scottbridle.com/
More informationTesting models as hypotheses: how can we be scientific in the face of epistemic errors?
UUEM Workshop, LSE, January 2014 Testing models as hypotheses: how can we be scientific in the face of epistemic errors? Keith Beven Lancaster University and Uppsala University My background I am a hydrologist
More informationChapter 5 : Design calculations
Chapter 5 : Design calculations 5. Purpose of design calculations For design calculations, it is important to assess the maximum flow, maximum water level and maximum velocity which occur in every node
More informationAvailable online at ScienceDirect. Procedia Engineering 119 (2015 ) 13 18
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 119 (2015 ) 13 18 13th Computer Control for Water Industry Conference, CCWI 2015 Real-time burst detection in water distribution
More informationBayesian Interpretations of Regularization
Bayesian Interpretations of Regularization Charlie Frogner 9.50 Class 15 April 1, 009 The Plan Regularized least squares maps {(x i, y i )} n i=1 to a function that minimizes the regularized loss: f S
More informationKan vi få et bedre miljø med smartere kloakker? Lektor Luca Vezzaro Forskning Døgn Slagelse, d. 26. april 2018
Kan vi få et bedre miljø med smartere kloakker? Lektor Luca Vezzaro Forskning Døgn Slagelse, d. 26. april 2018 Lidt om mig Født i Padova, tæt på Venedig Uddannet som miljøingeniør Kom til Danmark som udvekslingsstudent
More informationIntroduction to Imprecise Probability and Imprecise Statistical Methods
Introduction to Imprecise Probability and Imprecise Statistical Methods Frank Coolen UTOPIAE Training School, Strathclyde University 22 November 2017 (UTOPIAE) Intro to IP 1 / 20 Probability Classical,
More informationCSC321 Lecture 18: Learning Probabilistic Models
CSC321 Lecture 18: Learning Probabilistic Models Roger Grosse Roger Grosse CSC321 Lecture 18: Learning Probabilistic Models 1 / 25 Overview So far in this course: mainly supervised learning Language modeling
More informationIntroduction to Machine Learning CMU-10701
Introduction to Machine Learning CMU-10701 2. MLE, MAP, Bayes classification Barnabás Póczos & Aarti Singh 2014 Spring Administration http://www.cs.cmu.edu/~aarti/class/10701_spring14/index.html Blackboard
More informationParametric Techniques Lecture 3
Parametric Techniques Lecture 3 Jason Corso SUNY at Buffalo 22 January 2009 J. Corso (SUNY at Buffalo) Parametric Techniques Lecture 3 22 January 2009 1 / 39 Introduction In Lecture 2, we learned how to
More informationLecture Notes 1: Decisions and Data. In these notes, I describe some basic ideas in decision theory. theory is constructed from
Topics in Data Analysis Steven N. Durlauf University of Wisconsin Lecture Notes : Decisions and Data In these notes, I describe some basic ideas in decision theory. theory is constructed from The Data:
More informationDS-GA 1002 Lecture notes 11 Fall Bayesian statistics
DS-GA 100 Lecture notes 11 Fall 016 Bayesian statistics In the frequentist paradigm we model the data as realizations from a distribution that depends on deterministic parameters. In contrast, in Bayesian
More informationUncertainty and Bias UIUC, 403 Advanced Physics Laboratory, Fall 2014
Uncertainty and Bias UIUC, 403 Advanced Physics Laboratory, Fall 2014 Liang Yang* There are three kinds of lies: lies, damned lies and statistics. Benjamin Disraeli If your experiment needs statistics,
More informationCOMP61011 : Machine Learning. Probabilis*c Models + Bayes Theorem
COMP61011 : Machine Learning Probabilis*c Models + Bayes Theorem Probabilis*c Models - one of the most active areas of ML research in last 15 years - foundation of numerous new technologies - enables decision-making
More informationUncertainty Quantification for Machine Learning and Statistical Models
Uncertainty Quantification for Machine Learning and Statistical Models David J. Stracuzzi Joint work with: Max Chen, Michael Darling, Stephen Dauphin, Matt Peterson, and Chris Young Sandia National Laboratories
More informationE. Santovetti lesson 4 Maximum likelihood Interval estimation
E. Santovetti lesson 4 Maximum likelihood Interval estimation 1 Extended Maximum Likelihood Sometimes the number of total events measurements of the experiment n is not fixed, but, for example, is a Poisson
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Course Contents Introduction to Random Variables (RVs) Probability Distributions
More informationManipulating Radicals
Lesson 40 Mathematics Assessment Project Formative Assessment Lesson Materials Manipulating Radicals MARS Shell Center University of Nottingham & UC Berkeley Alpha Version Please Note: These materials
More informationBayesian Classifiers and Probability Estimation. Vassilis Athitsos CSE 4308/5360: Artificial Intelligence I University of Texas at Arlington
Bayesian Classifiers and Probability Estimation Vassilis Athitsos CSE 4308/5360: Artificial Intelligence I University of Texas at Arlington 1 Data Space Suppose that we have a classification problem The
More informationBayes Formula. MATH 107: Finite Mathematics University of Louisville. March 26, 2014
Bayes Formula MATH 07: Finite Mathematics University of Louisville March 26, 204 Test Accuracy Conditional reversal 2 / 5 A motivating question A rare disease occurs in out of every 0,000 people. A test
More informationEnsembles, Uncertainty and Climate Projections. Chris Brierley (Room 117)
Ensembles, Uncertainty and Climate Projections Chris Brierley (Room 117) Definitions Ensemble: group of model simulations Uncertainty: doubt and ambiguity about future conditions Climate Projection: modeled
More informationFORECASTING: A REVIEW OF STATUS AND CHALLENGES. Eric Grimit and Kristin Larson 3TIER, Inc. Pacific Northwest Weather Workshop March 5-6, 2010
SHORT-TERM TERM WIND POWER FORECASTING: A REVIEW OF STATUS AND CHALLENGES Eric Grimit and Kristin Larson 3TIER, Inc. Pacific Northwest Weather Workshop March 5-6, 2010 Integrating Renewable Energy» Variable
More informationInnovative Ways to Monitor Land Displacement
ARTICLE Innovative Ways to Monitor Land Displacement When people think about landslides, they usually imagine large mud streams which cause considerable loss of life. Whereas such large-scale disasters
More informationIntroduction to Machine Learning
Introduction to Machine Learning CS4375 --- Fall 2018 Bayesian a Learning Reading: Sections 13.1-13.6, 20.1-20.2, R&N Sections 6.1-6.3, 6.7, 6.9, Mitchell 1 Uncertainty Most real-world problems deal with
More informationAdvanced /Surface Hydrology Dr. Jagadish Torlapati Fall 2017 MODULE 2 - ROUTING METHODS
Routing MODULE - ROUTING METHODS Routing is the process of find the distribution of flow rate and depth in space and time along a river or storm sewer. Routing is also called Flow routing or flood routing.
More informationMATH 10 INTRODUCTORY STATISTICS
MATH 10 INTRODUCTORY STATISTICS Ramesh Yapalparvi It is Time for Homework! ( ω `) First homework + data will be posted on the website, under the homework tab. And also sent out via email. 30% weekly homework.
More information9/12/17. Types of learning. Modeling data. Supervised learning: Classification. Supervised learning: Regression. Unsupervised learning: Clustering
Types of learning Modeling data Supervised: we know input and targets Goal is to learn a model that, given input data, accurately predicts target data Unsupervised: we know the input only and want to make
More informationConfidence Intervals. First ICFA Instrumentation School/Workshop. Harrison B. Prosper Florida State University
Confidence Intervals First ICFA Instrumentation School/Workshop At Morelia,, Mexico, November 18-29, 2002 Harrison B. Prosper Florida State University Outline Lecture 1 Introduction Confidence Intervals
More informationIntroduction to Machine Learning
Uncertainty Introduction to Machine Learning CS4375 --- Fall 2018 a Bayesian Learning Reading: Sections 13.1-13.6, 20.1-20.2, R&N Sections 6.1-6.3, 6.7, 6.9, Mitchell Most real-world problems deal with
More informationStatistics for extreme & sparse data
Statistics for extreme & sparse data University of Bath December 6, 2018 Plan 1 2 3 4 5 6 The Problem Climate Change = Bad! 4 key problems Volcanic eruptions/catastrophic event prediction. Windstorms
More informationDer SPP 1167-PQP und die stochastische Wettervorhersage
Der SPP 1167-PQP und die stochastische Wettervorhersage Andreas Hense 9. November 2007 Overview The priority program SPP1167: mission and structure The stochastic weather forecasting Introduction, Probabilities
More informationHow to integrate wetland processes in river basin modeling? A West African case study
How to integrate wetland processes in river basin modeling? A West African case study stefan.liersch@pik-potsdam.de fred.hattermann@pik-potsdam.de June 2011 Outline Why is an inundation module required?
More informationComputational Perception. Bayesian Inference
Computational Perception 15-485/785 January 24, 2008 Bayesian Inference The process of probabilistic inference 1. define model of problem 2. derive posterior distributions and estimators 3. estimate parameters
More information12 Statistical Justifications; the Bias-Variance Decomposition
Statistical Justifications; the Bias-Variance Decomposition 65 12 Statistical Justifications; the Bias-Variance Decomposition STATISTICAL JUSTIFICATIONS FOR REGRESSION [So far, I ve talked about regression
More informationMachine Learning. Gaussian Mixture Models. Zhiyao Duan & Bryan Pardo, Machine Learning: EECS 349 Fall
Machine Learning Gaussian Mixture Models Zhiyao Duan & Bryan Pardo, Machine Learning: EECS 349 Fall 2012 1 The Generative Model POV We think of the data as being generated from some process. We assume
More informationModeling Great Britain s Flood Defenses. Flood Defense in Great Britain. By Dr. Yizhong Qu
Modeling Great Britain s Flood Defenses AIRCurrents Editor s note: AIR launched its Inland Flood Model for Great Britain in December 2008. The hazard module captures the physical processes of rainfall-runoff
More informationApplication of Real-Time Rainfall Information System to CSO control. 2 October 2011 Naruhito Funatsu METAWATER Co., Ltd.
Application of Real-Time Rainfall Information System to CSO control 2 October 2011 Naruhito Funatsu METAWATER Co., Ltd. Presentation Points Objectives To verify the applicability of the real-time rainfall
More informationStatistical Methods for SVM
Statistical Methods for SVM Support Vector Machines Here we approach the two-class classification problem in a direct way: We try and find a plane that separates the classes in feature space. If we cannot,
More informationMAE 493G, CpE 493M, Mobile Robotics. 6. Basic Probability
MAE 493G, CpE 493M, Mobile Robotics 6. Basic Probability Instructor: Yu Gu, Fall 2013 Uncertainties in Robotics Robot environments are inherently unpredictable; Sensors and data acquisition systems are
More informationFlood Forecasting. Fredrik Wetterhall European Centre for Medium-Range Weather Forecasts
Flood Forecasting Fredrik Wetterhall (fredrik.wetterhall@ecmwf.int) European Centre for Medium-Range Weather Forecasts Slide 1 Flooding a global challenge Number of floods Slide 2 Flooding a global challenge
More informationClimate Adaptation Challenges for Boston s Water and Sewer Systems
National Association of Flood & Stormwater Management Agencies Climate Adaptation Challenges for Boston s Water and Sewer Systems John P Sullivan P.E. October 15,2014 Boston 1630 Boston 1630-2012 Boston
More informationBEST ESTIMATE PLUS UNCERTAINTY SAFETY STUDIES AT THE CONCEPTUAL DESIGN PHASE OF THE ASTRID DEMONSTRATOR
BEST ESTIMATE PLUS UNCERTAINTY SAFETY STUDIES AT THE CONCEPTUAL DESIGN PHASE OF THE ASTRID DEMONSTRATOR M. Marquès CEA, DEN, DER F-13108, Saint-Paul-lez-Durance, France Advanced simulation in support to
More informationWorkshop on MCCOE Radar Meteorology /Climatology in Indonesia. Segel Ginting Wanny K. Adidarma
Workshop on MCCOE Radar Meteorology /Climatology in Indonesia BPPT, 28 Februari 2013 JAKARTA FLOOD EARLY WARNING SYSTEM (J-FEWS) Segel Ginting Wanny K. Adidarma JCP (Joint Cooperation Program) Indonesia
More informationComponents for Accurate Forecasting & Continuous Forecast Improvement
Components for Accurate Forecasting & Continuous Forecast Improvement An ISIS Solutions White Paper November 2009 Page 1 Achieving forecast accuracy for business applications one year in advance requires
More informationStrong Lens Modeling (II): Statistical Methods
Strong Lens Modeling (II): Statistical Methods Chuck Keeton Rutgers, the State University of New Jersey Probability theory multiple random variables, a and b joint distribution p(a, b) conditional distribution
More informationParameter Estimation. William H. Jefferys University of Texas at Austin Parameter Estimation 7/26/05 1
Parameter Estimation William H. Jefferys University of Texas at Austin bill@bayesrules.net Parameter Estimation 7/26/05 1 Elements of Inference Inference problems contain two indispensable elements: Data
More informationFeature selection. Micha Elsner. January 29, 2014
Feature selection Micha Elsner January 29, 2014 2 Using megam as max-ent learner Hal Daume III from UMD wrote a max-ent learner Pretty typical of many classifiers out there... Step one: create a text file
More informationQuantile POD for Hit-Miss Data
Quantile POD for Hit-Miss Data Yew-Meng Koh a and William Q. Meeker a a Center for Nondestructive Evaluation, Department of Statistics, Iowa State niversity, Ames, Iowa 50010 Abstract. Probability of detection
More informationIntroduction to AI Learning Bayesian networks. Vibhav Gogate
Introduction to AI Learning Bayesian networks Vibhav Gogate Inductive Learning in a nutshell Given: Data Examples of a function (X, F(X)) Predict function F(X) for new examples X Discrete F(X): Classification
More informationHydrological forecasting and decision making in Australia
Hydrological forecasting and decision making in Australia Justin Robinson, Jeff Perkins and Bruce Quig Bureau of Meteorology, Australia The Bureau's Hydrological Forecasting Services Seasonal Forecasts
More informationHow advances in atmospheric modelling are used for improved flood forecasting. Dr Michaela Bray Cardiff University
How advances in atmospheric modelling are used for improved flood forecasting Dr Michaela Bray Cardiff University Overview of current short term rainfall forecasting Advancements and on going research
More informationChapter-1 Introduction
Modeling of rainfall variability and drought assessment in Sabarmati basin, Gujarat, India Chapter-1 Introduction 1.1 General Many researchers had studied variability of rainfall at spatial as well as
More informationRegression with correlation for the Sales Data
Regression with correlation for the Sales Data Scatter with Loess Curve Time Series Plot Sales 30 35 40 45 Sales 30 35 40 45 0 10 20 30 40 50 Week 0 10 20 30 40 50 Week Sales Data What is our goal with
More informationA Support Vector Regression Model for Forecasting Rainfall
A Support Vector Regression for Forecasting Nasimul Hasan 1, Nayan Chandra Nath 1, Risul Islam Rasel 2 Department of Computer Science and Engineering, International Islamic University Chittagong, Bangladesh
More informationChoosing among models
Eco 515 Fall 2014 Chris Sims Choosing among models September 18, 2014 c 2014 by Christopher A. Sims. This document is licensed under the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported
More informationImproving Operational Use of Scanning Rain Radar Estimates with Vertically Pointing Radar
Improving Operational Use of Scanning Rain Radar Estimates with Vertically Pointing Radar Luke Sutherland-Stacey Weather Radar NZ www.weatherradar.co.nz Tom Joseph Mott MacDonald Geoff Austin Weather Radar
More informationDEVELOPMENT OF A FORECAST EARLY WARNING SYSTEM ethekwini Municipality, Durban, RSA. Clint Chrystal, Natasha Ramdass, Mlondi Hlongwae
DEVELOPMENT OF A FORECAST EARLY WARNING SYSTEM ethekwini Municipality, Durban, RSA Clint Chrystal, Natasha Ramdass, Mlondi Hlongwae LOCATION DETAILS AND BOUNDARIES ethekwini Municipal Area = 2297 km 2
More information