returnpotpareto {INQC}

Threshold percentile for the Pareto outliers

Description

This function returns a value of a threshold percentile for the Pareto outliers

Usage

returnpotpareto(pato, ret, w = 1.65)

Arguments

pato

list with results of modelling/fitting the generalized Pareto distribution

ret

pseudo-return period (in yr)

w

average sampling frequency (in 1/yr), a parameter to equate to return period to a temporal interval (recall the approach is not block maxima but peak over threshold. Typical value of w to equate the return period to years is 1.65 (See Wilks, 2011. Statistical Analysis for the Atmospheric Sciences)

Value

for a given Pareto distribution, returns the value (the quantile) representing a requested return period

Examples

#Extract the ECA&D precipitation data file from the example data folder
path2inptfl<-system.file("extdata", "RR_SOUID132730.txt", package = "INQC")
#Read the data file
y<-readecad(input=path2inptfl,missing= -9999)[,4]
#Fit the Generalized Pareto distribution
pato<-potpareto(y)
#Define the quantile corresponding to the requested return period 25 years (ret=25)
returnpotpareto(pato,25)
#Define the quantile assuming the existence of 2 precipitation peaks/extreme values 
#every year (on average)
returnpotpareto(pato,25,w=2)