Mt Rainier

Mt Rainier
Mt Rainier
Showing posts with label stochastic modeling. Show all posts
Showing posts with label stochastic modeling. Show all posts

Saturday, May 21, 2011

Calculating Odds and Probabilities



1986 Liberty Half Dollar Tossed on Grass (Heads) (Image on Photoshelter)


A team will have an opportunity to call “heads or tails” for a coin toss at the beginning of a football game. If the call is successful, heads in the case of the coin pictured above, the team will have the opportunity to chose to take the ball (go on offense) or defend a particular goal. The choices are non-trivial as a team’s defense may be superior to their offense, or the wind may determine the choice of the goal defended. These choices represent the team’s expectations of their success, given their strengths, weaknesses and prevailing conditions.

The coin toss, for a fair coin, will give a probability (p) of .5, or 50% that the call (either heads or tails) will be successful. Coins flipped successively will tend, over a period of time to “revert to the mean”, or tend towards their expected value of 50% heads and 50% tails. The odds ratio for the coin toss is 1:1. Odds ratio is calcuated as the ratio of the probability of the event happening (p) to the probability that it won’t happen (1-p). Thus, out of 2 possible outcomes (heads or tails), one outcome (out of 2) is successful, one outcome (out of 2) is unsuccessful, for a ratio of 1:1.

While the coin toss is random, there is a balance of risks involved in the result. For example, a team with a superior offense can choose to go on offense, sputter on three downs, punt, and turn the ball over to the other team in good field position, despite their expectation of offensive success.


Two Sixes Dice Toss (infrared digital effect) (image on Photoshelter)

People use probabilities and odd ratios in their daily lives, consciously and unconsciously computing expectations and making decisions based on those expectations.

Some may choose to gamble on the roll of the dice. In this case, the probability of rolling two sixes is (1/36) while the odds are 1 to 35 (1 good outcome, 35 losing outcomes). When gambling in an establishment, the house will include a margin in for themselves, to cover expenses and profit. Given fair dice, the underlying odds (and the probabilities) will tend towards their mathematical expected value (1 to 35 and 1/36), thus giving over time a margin to the house.

Mathematical models can calculate the odds that there will be a city rivalry World Series this year, the odds of a particular candidate winning the 2012 election, and many other events.

Of particular interest are the mathematical models that deal with environmental issues, such as weather and climate change, given the emergence of extreme weather events.

Albert Einstein was quoted as saying “I, at any rate, am convinced that He does not throw dice.” (paraphrased as “God does not play dice with the universe”). However Stephen Hawking is quoted as saying “So Einstein was wrong when he said, "God does not play dice." Consideration of black holes suggests, not only that God does play dice, but that he sometimes confuses us by throwing them where they can't be seen”.

The National Oceanographic and Atmospheric Administration (NOAA) has sophisticated climatological models which forecast the weather, and they are continually improving these models in order to provide the maximum amount of warning possible to the populace. Through their efforts, NOAA is seeking to minimize the amount of risk, and maximize the information provided, seeking the “load the dice” in the publics favor.

We rely on weather forecasts in our daily plans. Should we call off the company picnic with a 70% chance of thunderstorms? How about that outdoor wedding? Do you pack the rain slicker or umbrella or leave it behind? Governments and businesses rely on weather forecasts for their daily operations and long term plans. Seasonal hurricane forecasts are looked at with interest by many sectors in planning ahead into the coming hurricane season.



Two United States Coast Guard Cutters monitoring the Westport, Washington Bar entrance (image on Photoshelter)

Nowhere is the issue of balance of risks and weather forecasting more apparent than in the maritime arena. When is the storm predicted to arrive? Is there a confidence interval around that arrival time, as well as its expected intensity? Sailors will balance the risk of venturing out and making headway to their destination versus the risk that the weather will take an unexpected turn for the worse, making for difficult, or hazardous headway, or worse.

The United States Coast Guard (USCG) is a branch of the U.S. Armed forces and is charged with maritime safety, security and stewardship. In this image, the two Coast Guard Cutters are monitoring the closed Westport Bar / navigational entrance from the Pacific Ocean to Westport, Washington. The USCG monitors ongoing conditions and weather forecasts to decide whether or not to close the bar, thus managing safety for mariners. The United States Coast Guard, in its decision-making, relies on the concept of balance of risk, as it makes decisions that affect the public. The public, in turn, makes decisions using balance of risk as it relies on NOAA guidance, the actions of the Coast Guard, and other factors.

It is human nature to calculate probabilities, calculate odds, consciously or unconsciously assess balance of risks. We do so continuously in our daily lives and rely on others to help us do so. We need to expand that awareness to include the risks we face due to environmental changes, consider the balance of risks inherent in the evolving situation, and find a way to meet the considerable challenges head on.

Tuesday, April 26, 2011

Chernobyl 25th Anniversary



An aerial view of Chernobyl Nuclear Power Plant in April, 1986, with the red glow towards the center showing the heat from Unit #4. Source: epa.gov

The Chernobyl Disaster was a nuclear accident which occurred on April 26, 1986. It occurred in the Ukraine Republic (formerly part of the former USSR (Soviet Union)). Now the Chernobyl site is part of the country Ukraine. Today marks the 25th anniversary of the Chernobyl disaster.

It is particularly compelling to consider the impacts of Chernobyl today, twenty-fives years later, as we witness the unfolding of another nuclear disaster at Fukushima, Japan, following the 9.0 earthquake and tsunami on March 11, 2011. The two disasters arose from different circumstances and unfolded differently, however they share in common the impact of a low probability-high risk event.

They are set apart in both time and space, one occurring in vastness of the then-Soviet Union; the other set on the more densely populated island of Japan. The Chernobyl Disaster and the Fukushima disasters were both were graded as a 7 on the International Nuclear Event Scale. (Fukushima was raised from a 5 to a 7 on April 12, 2011, one month and a day following the earthquake and tsunami on March 11th).

How do we fathom such events as we seek to understand the risks of pursuing nuclear energy? How do we internalize these low probability-high risk events so that we carefully assess risks yet do not fall prey to unwarranted fears and suspicions? Do we hide our heads in the sands of improbability and ignore the potential of a very small yet very dangerous risk? Do we pour millions and billions of dollars to hedge against a risk that may never, in our lifetime occur? Do we even care about the impact of our decisions on those ancestors who may follow us many generations down the turnpike?

It is critical how we answer these questions, because our fate and the fate of our planet may hang in the balance. With the burgeoning population on this planet, the growing scarcity of resources, and the challenges presented by rising carbon dioxide levels, and other indications of planetary strain, we must find a way to make informed decisions that appropriately incorporate the low probability/high risk event in our search for and use of energy resources.

In my last blog posting, Emperor Penguin Energy-Risk Model - Part 2 , I discussed mathematical modeling, random variables and evolved a stochastic model of emperor penguin energy-risk behavior. I discussed some of the variables that may be considered in the emerging emperor penguin population, including mortality, morbidity and accident. I introduced the concept of a low probability event into the model (an eruption of an Antarctic volcano), and discussed the values of the stochastic process in informing results. The post was intended to discuss energy seeking behavior in a different (emperor penguin) population as an energy seeking risk example using stochastic modeling.

In my last post I stated “The objective of the stochastic processes is to help us inform our decision making process, to help us understand the impact of variables under a wide range of assumptions, conditions, and scenarios. Thus a stochastic process should inform us about the expectations of the model under a wide variety of conditions, including the impacts of low probability / high risk events.”

However, a stochastic process cannot inform without reasonable assumptions. Assumptions must be developed to allow the stochastic process to produce a credible range of results that will indeed be informative for the intended usages. There are many variables to consider, and assumptions to be made in analyzing risk. For a variety of reasons it may be difficult to obtain a robust set of assumptions that everyone agrees with for all potential uses.

Thursday, April 21, 2011

Emperor Penguin Energy-Risk Model - Part 2



Emperor Penguin Diving onto Ice Shelf from Sea, Stancomb Wells Ice Edge, Weddell Sea, Antarctica (Image on Alamy.com)

In my last blog post “An Emperor Penguin Energy-Risk Model” on April 14, 2011, l discussed the predator-prey relationship between the leopard seal and emperor penguin in Antarctica. The leopard seal waits at the edge of the ice shelf and opportunistically picks off emperor penguins entering or leaving the sea. For the emperor penguin, feeding at sea is a decision between the need to feed to live and the risk of dying in the mouth of a leopard seal..

In the blog post I state: “From studying the emperor penguin and the leopard seal we know the emperor penguins will continue to feed, but so will the leopard seal. Some emperor penguins, despite their various risk protection strategies, will get eaten. It is important to note that in a probabalistic sense, we know that some penguins will be eaten by the leopard seal, but we don’t know which specific penguins will “bite the dust”.” This is true casually looking at a row of emperor penguins lined up to go into the sea in search of food.

However, upon closer analysis and study, over a period of time, it might be possible to determine which emperor penguins have a bit of catch in their step, have been injured in a narrow escape from a leopard seal, or have slowed down. These emperor penguins might come a belly-flop short of landing on the ice, and end up as prey in the mouth of a leopard seal. However, it is also possible, that a healthy, fit, member of the emperor penguin colony might suffer a particularly ill-fated episode of bad luck. This penguin might be in the wrong place at the wrong time when the leopard seal is rising out of the water with its mouth wide open ready for business. In fact, you could have the emperor penguin equivalent of the 4.0-40 yard dash champion, and end up as leopard seal “dinner”, with some bad luck and timing.

Looking at emperor penguin energy-seeking behavior and risk, it becomes apparent that probabilities have a great deal to do with the outcome but are not deterministic. You may attach a relatively higher probability of being eaten to the more fragile members of the emperor penguin population and a relatively lower probability of being eaten to those fitter members. The larger the colony size, and the more emperor penguins entering the sea at the same time, the lower the risk, the probability of being eaten, for any particular emperor penguin as there are more penguins entering the sea. (“there’s safety in numbers”).

You can run scenarios with differing proportions of fragile and fit emperor penguins, with higher and lower probabilities of being eaten (mortality rates), varying degrees of illness (morbidity rates) or accident, including leopard seal attack. In such scenarios, the leopard seal would most likely pick off different emperor penguins each time the scenario is run, however there would be objective tendencies to pick off more members of the more fragile group versus those of the fitter group.

In performing mathematical modeling of the fate of the emperor penguins by running scenarios with objective data and assumptions, we may set up a stochastic process which helps us to understand the behavior of the system as it evolves under a variety of scenarios.

Mathematical models involve expressing real world problems in mathematical language. This entails defining variables and establishing a formulaic process which will express the model as evolves. Variables are elements in the model which may change during the model. Because they may change, the model needs to calculate how they change over the course of the model and how they interact with other model variables, and are affected by the constants assumed by the model. Constants may arise from established data or may be assumptions plugged in to the model.

Stochastic processes incorporate non-deterministic, random elements into a mathematical model. The result may vary with time and with each model run. In comparison, a deterministic model will always produce the same result given the same assumptions and initial state. Thus, a stochastic process is run using random processes, employing a variety of assumptions and probability distributions informing objective tendencies for various model events to occur..

The random process in stochastic modeling will randomly choose which penguins are attacked, survive, suffer morbidity or injury from accident, and die over a period of time. Each run will be unique, as specific, members of the colony are differently impacted by the random process each time. By running many such models, one can get a picture of the survival data for the colony as a whole under a wide range of assumptions. Depending on the characteristics of the data, model and variables, results may be similar on an overall group basis, while differing by individual members impacted over time.

Under a normal range of assumptions and outcomes, this model may well predict overall group behavior over a period time. However modeling becomes much challenging when very low probability events enter into the model or rear their head in actual life.

For example, a eruption of an Antarctic volcano may be infrequent, however it could certainly impact emperor penguins. If the model assumed a volcanic eruption with a low probability, a robust number of stochastic model runs may randomly select such an event resulting in a BBQ penguin supper for the leopard seals.

The objective of the stochastic processes is to help us inform our decision making process, to help us understand the impact of variables under a wide range of assumptions, conditions, and scenarios. Thus a stochastic process should inform us about the expectations of the model under a wide variety of conditions, including the impacts of low probability / high risk events.