# Get an annual maximum sample, estimate the parameters, and estimate the 50-year RP
am_27090 <- GetAM(27090)
pars <- as.numeric(GumbelPars(am_27090$Flow))
GumbelEst(pars[1], pars[2], RP = 50)
# Estimate the RP for a 600 m^3/s discharge
GumbelEst(pars[1], pars[2], q = 600)
Run the code above in your browser using DataLab