GEVAM: Generalised extreme value distribution - estimates directly from sample
Description
Estimated quantiles as function of return period (RP) and vice versa, directly from the data
Usage
GEVAM(x, RP = 100, q = NULL)
Value
quantile as a function of RP or vice versa.
Arguments
x
numeric vector (block maxima sample)
RP
return period (default = 100)
q
quantile (magnitude of variable)
Author
Anthony Hammond
Details
If the argument q is used, it overrides RP and provides RP as a function of q (magnitude of variable) as opposed to q as a function of RP. The parameters are estimated by the method of L-moments, as detailed in 'Hosking J. Wallis J. 1997 Regional Frequency Analysis: An Approach Based on L-moments. Cambridge University Press, New York'.
#Get an annual maximum sample and estimate the 50-year RPAM.27090 <- GetAM(27090)
GEVAM(AM.27090$Flow, RP = 50)
#Estimate the RP for a 600m3/s dischargeGEVAM(AM.27090$Flow, q = 600)