I am using the "lmomRFA" package for frequency analysis. I have the best fitted distribution for the region when doing quantile estimation for return periods.
My problem is that i need the standard deviation for the quantile estimates for the different return periods, i use this code but returns "Error in pelxxx("gev", lmom) : pelgev: missing values in L-moment vector"
Code: regsimq(qfunc = quagev, para = c(0.8254, 0.2492, -0.1118), cor = 0, index = NULL, cod, nrep = 1000, fit = "gev", f = c(0.01, 0.8, 0.9, 0.98, 0.99), boundprob = c(0.05, 0.95), save = TRUE)
Does anyone know whats wrong with the code?