EncProb | R Documentation |
Calculates the probability of experiencing at least n events with a given return period (RP), over a given number of years
EncProb(n, yrs, RP, dist = "Poisson")
n |
number of events |
yrs |
number of years |
RP |
return period of the events |
dist |
choice of probability distribution. Either "Poisson" or "Binomial" |
The choice of binomial or Poisson distributions for calculating encounter probablities is akin to annual maximum (AM) versus peaks over threshold (POT) approaches to extreme value analysis. AM and binomial assume only one "event" can occur in the blocked time period. Whereas Poisson and POT don't make this assumption. In the case of most catchments in the UK, it is rare to have less than two independent "events" per year; in which case the Poisson and POT choices are more suitable. In large catchments, with seasonally distinctive baseflow, there may only be one independent peak in the year. However, the results from both methods converge with increasing magnitude, yielding insignificant difference beyond a 20-year return period.
A probability
Anthony Hammond
#Calculate the probability of exceeding at least one 50-yr RP event
#over a 10 year period, using the Poisson distribution.
EncProb(n = 1, yrs = 10, RP = 50)
#Calculate the probability of exceeding at least two 100-yr RP events
#over a 100 year period, using the binomial distribution.
EncProb(n = 2, yrs = 100, RP = 100, dist = "Binomial")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.