Lecture 08
September 23, 2024
Text: VSRIKRISH to 22333
Glib Answer: A lack of certainty!
More Seriously: Uncertainty refers to an inability to exactly describe current or future values or states.
Two (broad) types of uncertainties:
Designing and managing environmental systems is often about minimizing or managing risk:
The Society for Risk Analysis definition:
“risk” involves uncertainty about the effects/implications of an activity with respect to something that humans value (such as health, well-being, wealth, property or the environment), often focusing on negative, undesirable consequences.
Important: “Risk” is not just another words for probability, but:
Multiple components which contribute to risk:
Risk management is often a key consideration in systems analysis. For example, consider regulatory standards.
We often represent uncertainties with probabilities. What is a probability?
We often don’t want to just know if a particular event \(A\) has a certain probability, but also how other events (call them \(B\)) might depend on that outcome.
In other words:
We want the conditional probability of \(B\) given \(A\), denoted \(\mathbb{P}(B|A)\).
We can write conditional probabilities in terms of unconditional probabilities:
\[\mathbb{P}(B|A) = \frac{\mathbb{P}(BA)}{\mathbb{P}(A)}.\]
Conditional probabilities can be inverted according to Bayes’ Theorem:
\[\mathbb{P}(A|B) = \frac{\mathbb{P}(B|A) \times \mathbb{P}(A)}{\mathbb{B}}.\]
The probability of possible values of an unknown quantity are often represented as a probability distribution.
Probability distributions associate a probability to every event under consideration (the event space) and have to follow certain rules (for example, total probability = 1).
The specification of distributions can strongly influence the analysis.
A distribution implicitly answers questions like:
The tails of distributions represent the probability of high-impact outcomes.
Key consideration: Small changes to these (low) probabilities can greatly influence risk.
Monte Carlo simulation: Propagating random samples through a model to estimate a value (usually an expectation or a quantile).
Monte Carlo is a broad method, which can be used to:
Monte Carlo estimation involves framing the quantity of interest as a summary statistic (such as an expected value).
Finding \(\pi\) by sampling random values from the unit square and computing the fraction in the unit circle. This is an example of Monte Carlo integration.
\[\frac{\text{Area of Circle}}{\text{Area of Square}} = \frac{\pi}{4}\]
What is the probability of rolling 4 dice for a total of 19?
Can simulate dice rolls and find the frequency of 19s among the samples.
This type of estimation can be repeated with any simulation model that has a stochastic component.
For example, consider our dissolved oxygen model. Suppose that we have a probability distribution for the inflow DO.
How could we compute the probability of DO falling below the regulatory standard somewhere downstream?
Wednesday: Monte Carlo Lab (clone before class, maybe instantiate environment too)
Monday: Why Does Monte Carlo Work?
HW3: Released today, due 10/03 at 9pm.