Generalised Pareto distribution estimates from parameters
Estimated quantiles as function of return period (RP) and vice versa, from user input parameters
GenParetoEst(loc, scale, shape, q = NULL, RP = 100, ppy = 1)
loc |
location parameter |
scale |
scale parameter |
shape |
shape parameter |
q |
quantile. magnitude of the variable under consideration |
RP |
return period |
ppy |
peaks per year. Default is one |
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 average number of peaks per year argument (ppy) is for the function to convert from the peaks over threshold (POT) scale to the annual scale. For example, if there are 3 peaks per year, the probability associated with the 100-yr return period estimate would be 0.01/3 (i.e. an RP of 300 on the POT scale) rather than 0.01.
quantile as a function of RP or vice versa
Anthony Hammond
#Get a POT sample, estimate the parameters, and estimate 50-year RP ThamesPOT <- POTextract(ThamesPQ[,c(1,3)], thresh = 0.90) GenParetoPars(ThamesPOT$peak) GenParetoEst(174.2862, 127.4085, 0.1805716, ppy = 1.867, RP = 50) #Estimate the RP for a 600m3/s discharge GenParetoEst(174.2862, 127.4085, 0.1805716, ppy = 1.867, q = 600)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.