Page History

...

The plotting of ordered data on extreme probability paper is done according to a general plotting position function: P = (m - a) / (N + 1 - 2a). Constant 'a' is an input variable and is default set to 0.3. Many different plotting functions are used, some of them can be reproduced by changing the constant 'a'.

For the normal distribution the Blom plotting function can best be used, for the Gumbel distribution Gringorton gives the best results. A plotting function which can be used for all distribution functions is Chegadayev. More information on plotting positions can be found in many hydrological handbooks (i.e. Applied Hydrology: pages. 394-396).
The frequency distributions are briefly described in then next Section. For a detailed description of the frequency distributions reference is made to other text books. The selection procedure of data is dealt with in Section 10.2.3.

Frequency

...

Distributions

Normal

...

Distribution Function:

with:
m_{x μx} = mean of X,
s σ_x = standard deviation of X.
The estimates for m_x and sσ_x are obtained from the equations explained in the Basic Statistics function.

Log-normal

...

Distribution Function:

Image Added!05 - Time Series Analysis^image026.gif!
with:
Y = ln (X - X₀)
m μ_x = mean of X,
s_x = standard deviation of X.
X₀ = location parameter, conditioned by (X₀ < X).

The location parameter X₀ may be determined (3-parameter distribution) or can be provided by the user (2-parameter distribution). The model parameters can be estimated by: 1.

1. method of moments, or

...

1. modified maximum likelihood method.

Box-Cox

...

Transformation to

...

Normality:

Arbitrarily distributed variables may be transformed to normality by the following transformation:

for l ¹0
for l= 0

where:
X₀ = location parameter, (X₀ < X)
l λ = transformation parameter,

The location parameter X₀ is provided by the user while the other parameters are estimated by the maximum likelihood method.

Standard

...

Iincomplete Gamma Function:

This function forms the basis of a number of distribution functions. It has the following form: !05 - Time Series Analysis^image032.gif!

Image Added

with:
!05 - Time Series Analysis^image034.gif! Image Added

where: 
g γ = shape parameter
Z = argument

Pearson Type III or Gamma

...

Distribution Function:

P_PIII (X) = F_G (Z)

with:
Z = (X - X₀)/ bβ
X₀ = location parameter, (X₀ < X)
b β = scale parameter.

The location parameters X₀ may be provided (2-parameter distribution) or estimated (3-parameter distribution). The parameters X₀, b and gcan be estimated by:

1. method of moments, and
2. modified maximum likelihood method

Raleigh

...

Distribution Function:

P_R (X) = F_G (Z)

with:
Z = ((X - X₀) / bβ )² and g and γ = 1
X₀ = location parameter, (X₀ < X)
b β = scale parameter.

The parameters X₀ and bare estimated from the moments.

Exponential

...

Distribution Function:

P_E (X) = F_G (Z)

with:
Z = (X - X₀) / b and g  β and γ = 1
X₀ = location parameter, (X₀ < X)
b β = scale parameter.

The parameters X₀ and bare estimated from the moments.

General Pearson

...

Distribution Function:

P_GP (X) = 0.5 + (k/|k|)(F_G (Z) - 0.5)

with:
Z = ((X - X₀ )/bβ)^k
b β = scale parameter.
K = = type parameter (integer):
1 : Pearson Type III (Exponential distribution for g 1for γ = 1),
2 : Raleigh (g γ = 1) and Maxwell (gγ = 1.5) distributions,
-1 : Pearson Type V distribution.
The type parameter k is provided by the user. The parameters X₀, (X₀ < X) or g may be provided by the user (2-parameter distribution) or can be estimated (3-parameter distribution). The parameters X₀, band gare  β and γ are estimated by a mixed moment-maximum likelihood method.

Log-Pearson Type III

...

Distribution Function:

P_LP (X) = F_G (Z)

with:
Z = (ln(X - X₀ )-Y{~}0 ) / bβ,
X₀ ,Y₀ = location parameters,
b β = scale parameter.
g β = shape parameter.

The location parameter X₀, (X₀ < X) is provided by the user, while the parameters Y₀, band gare estimated by: 1.

1. mixed moment-maximum likelihood method,

...

1. modified maximum likelihood method on Y = ln (X - X₀)

Extreme Type I or Gumbel distribution:

...

with:
X₀ = location parameter,
b β = scale parameter.

The parameters X₀ and bcan be estimated by: 1.

1. method of moments,

...

1. modified maximum likelihood method, (with or without censoring).

The parameters can also be determined by application of the second method for part of the data set leaving the lowest N₁ and highest N₂ values out of the analysis provided:
N - (N₁ + N₂) ³ 5.

...

with:
X₀ = location parameter, (X₀ < X) !!
b β = scale parameter,
g γ = shape parameter. (k < 0)
The location parameter may be given (2-parameter distribution) or estimated (3-parameter distribution). The parameters can be estimated by the modified maximum likelihood method.

...

with:
X₀ = location parameter, (X₀ < X) !!
b β = scale parameter,
g γ = shape parameter. (k > 0)

For the estimation of parameters the same applies as for the Extreme Type II distribution.

...

with:
X₀ = location parameter, (X₀ < X) !!
b β = scale parameter,
g γ = shape parameter. (k > 0)

For the estimation of parameters the same applies as for the Extreme Type II distribution.

...

• actual values
• annual minimum values, and
• annual maximum values.

In case of actual values a threshold selection menu is displayed from which no threshold, minimum, maximum, both, peaks over threshold (POT) or peaks under threshold can be selected. If peaks over threshold is selected and no use will be made of the POT-method, a message will appear on screen, the number of values will be made equal to the number of years and the highest values will be selected. If peaks are selected also a value for the 'horizon' can be entered (a period in time interval units of the original series, e.g. days, which is used to skip lower peaks within that period before and after a peak). Default value = 1 (no horizon). A maximum of 1000 values/peaks of a maximum of 50000 input values will enter analysis.
The annual minimum/maximum values may be the minimum/maximum of full years or of a part of years, like seasonally or monthly extremes.

Computation Period

The computation period can be set to 'Full years' or 'Part of years'. When 'Full years' is selected all data values of the complete year will be taken into the analysis, make sure the start date and end date of the processing period are form the first of January to the first of January. When 'Part of years' is selected a start date and end date for the sub-period must be entered. When data for only the month of March must be selected, the start date of the sub-period must be set to "01-03" and the end date of the sub-period to "01-04".

...

For use of the Peaks over Threshold (pot)-method <actual values> have to be selected. In the potoption the number of years is input. In an earlier stage (see above) a threshold level X~0~ and a horizon are entered. The data series used in the method are all peaks between successive up-crossings and down-crossings taking into account the given horizon. For a more extensive explanation of the POT method, see the Basic Statistics function.

note : the amount of data is generally not equal to the number years. Hence the interpretation of probability of non-exceedance versus return period differs from the usual one.
References:

• Chow, V.T., Maidment, D.R., and Mays, L.W., Applied Hydrology, McGraw-Hill Book Co., 1988.