Hello, I am creating error bounds based on simulated data across multiple data durations. I was wondering if there is a way to update an object class from "regfit" from lmomRFA package? The reason is for consistency across durations. Example below: library(lmom); library(lmomRFA) data <- c(0.42, 0.13, 0.59, 0.12, 0.78, 0.17, 0.3, 0.41, 0.28, 0.79) # random data reg <- regsamlmu(data) # calc l-moments *org_gev <- regfit(reg,"gev")* # original gev fit # UPDATE *"org_gev*" values (below) and save as "update_gev" # xi = 0.65 # alpha = 0. 51 # k = -0.023 Thank you for the help, Doug -- ----------------------------------------- Douglas M. Hultstrand, MS Senior Hydrometeorologist Applied Weather Associates Monument, Colorado mobile: 720-771-5840 www.appliedweatherassociates.com dhultstrand at appliedweatherassociates.com ----------------------------------------- [[alternative HTML version deleted]]