************************************************** * * * lmomRFA package news * * * ************************************************** Version 3.6 2023 Aug 29 o Fortran source files rearranged, to make it clear which routines are copied from package "lmom", now specifically "lmom" version 3.0. o [invisible to users] Cleanup of package "lmomRFA"'s own Fortran code: modifications to avoid warnings from CRAN's compiler (GCC 12.2.0) about "Fortran 2018 deleted features"; specific names of intrinsic functions changed to generic. o Documentation of "regdata" now includes definition of L-CV. Version 3.5 2023 Feb 19 o New vignette "RegSim". o In regsimq(), regquantbounds(), and sitequantbounds(), there are changes in calculation methods for accuracy measures (RMSEs and error bounds). For the relative RMSE of the regional growth curve (computed by regsimq()), the absolute RMSE of the regional growth curve (computed by regquantbounds()), and for RMSEs and error bounds of quantiles at sites not used in regional estimation (computed by sitequantbounds() with 'index' and 'seindex' specified), the computational method is different from that used in previous versions of the package. The new approach is fully defined in vignette "RegSim". It uses more information from the simulations and should be more accurate, though results are typically little different from prior versions. The calculation method for RMSEs and error bounds of quantiles at sites used in regional estimation (computed by sitequantbounds() with 'index' not specified) is unchanged. Note also that in regsimq(), version 3.4, the relative RMSE of the regional growth curve was computed incorrectly. The relative RMSE was underestimated, though typically only by a small amount. The RMSE returned by regquantbounds() is a scaled version of the relative RMSE computed by regsimq() and was similarly affected. A fix of this issue was superseded by the change in calculation method noted above. o regsimq()'s return value gains elements "nrep", "true.asgc", and "sim.rgc", but no longer has an element "sim.rgcratio". o In sitequantbounds() when 'index' and 'seindex' are both specified, the distribution of the relative estimation error of the index flood is now assumed to be gamma rather than lognormal. This is felt (by the package maintainer) to be more realistic for most practical applications of regional frequency analysis. o [invisible to users] Changes to internal Fortran code: * Added routine SAMLM (from package "lmom", version 2.9). * Removed routines SAMLMU, SORT. * Minor changes to routines REGTST, LMRGPA. o Fixed: print.regsimq() of the result of a call to regsimq() with 'save=FALSE' would not print the number of simulations. o Fixed: regsimq() would give an error if 'boundprob' had only 1 element. Version 3.4 2022 May 29 o In the value returned by regfit(), component 'qfunc' (a function) now has a smaller environment and uses less memory. o Documentation of regsamlmu() notes that it gives different results from the corresponding program in the LMOMENTS Fortran package, version 3.04. o Fixed: regsimq() was computing the relative RMSE for the regional growth curve incorrectly. The error was small, typically less than 1%, and always in the positive direction. Version 3.3 2019 Dec 5 o Fixed: evplot.rfd(), when adding a site quantile function to an existing plot, would plot an incorrectly scaled version of the quantile function (it would plot the unscaled quantiles of the regional frequency distribution rather than applying the appropriate site-specific scale factor ("index flood")). o [invisible to users] Internal code no longer assumes that class() of a matrix has length 1 (would have caused error in R version 4.0.0). Version 3.2 2019 Mar 10 o Fixed: in cluinf(), incorrect test of whether argument 'merge' could be the result of a call to cluagg(). Version 3.1 2017 Aug 1 o "Native routines" are now "registered". o The simulation procedure used in regtst() involves fitting a kappa distribution to the regional data. The kappa fit may fail, in which case a generalized logistic distribution is used instead. This has always been the case, but is now documented. o Fixed: regtst.s() would give an error for a region containing a single site. Version 3.0-1 2015 Feb 1 o Maintainer's email address changed. No changes to code. Version 3.0 2014 Jan 13 o Increase in major version number, occasioned by some significant recent changes: * Fix to regsimq() noted below: earlier versions of the package would give incorrect results, though errors were typically less than 10%. * Change to computation of relative RMSE in version 2.5. Results returned by regsimq() in earlier versions of the package should be regarded as unreliable. o Now explicitly imports standard package "stats", and functions from standard package "graphics". o In regfit() and regsimq(), the search path used to find the functions specified in the arguments regfit(dist=) and regsimq(fit=) is now explicitly defined. (The search path was previously undocumented, and was not consistent across all prior versions of the package.) o regfit() and regsimq() give more informative error messages when the functions specified in their respective arguments 'dist' and 'fit' cannot be found. o regquantbounds() and sitequantbounds() should now never return negative error bounds. Instead, NA or Inf is returned for error bounds corresponding to negative or positive estimated quantiles, respectively. o regsimq() is faster, by about 25% in the author's tests. o Fixed: regquantbounds() and sitequantbounds() would give negative values for RMSEs corresponding to negative estimated quantiles. o Fixed: regsimq() computed incorrect values in its return value component "relbounds.by site". o Fixed: sitequantbounds(), when 'sitenames' is missing and 'index' has names, is documented to use the names of 'index' as the site names but did not in fact do so. Version 2.5 2013 Jun 17 o Now depends on package lmom, version >= 2.0. o New function regtst.s(). o In regsimq(), the computation of relative RMSE of estimates of quantiles of the regional growth curve has been modified, and is now more consistent with the error bounds for different probabilities. The interpretation of the relative RMSE has not changed: it is still the relative RMSE of the estimated regional growth curve regarded as an estimator of a randomly chosen site growth curve. o In 'regsimq(...,save=TRUE)', the matrix that is component 'sim.rgcratio' of the return value now has row names (derived from the value of the function's argument 'f'). o Component "D" of the return value from regtst() now has names (the site names). o LICENSE file no longer included in package, to avoid a NOTE when building the package for CRAN under R version 3.0.0. The license can be found at http://www.eclipse.org/legal/cpl-v10.html. Version 2.4 2012 Jan 23 o The package uses lazy loading of data, so its datasets (Appalach, Cascades, Maxwind) can be accessed without first calling data(). o regtst() now responds more quickly to a user interrupt (e.g. pressing Esc from the R Windows GUI). o Fixed: regtst() would give an error when the discordancy measure could not be calculated because the covariance matrix of the sample L-moment ratios (matrix A of Hosking and Wallis, 1997, eq. (3.2)) was singular. Version 2.3 2010 Aug 21 o regsimq() no longer gives an error when asked to compute the estimation accuracy for an infinite quantile; instead it returns NaN values in the relevant places (i.e., elements corresponding to infinite quantiles in the "rel.RMSE" and "rel.bound..." variables of the data frames in components "relbounds.rgc" and "relbounds.by.site"). o Fixed: when a data frame with only 6 columns was supplied as argument 'regdata' of regtst(), and summary() was called on the resulting object of class "regtst", an error would occur ("Error in quawak ... missing values in parameter vector"). Version 2.2 2009 Nov 16 o evplot.rfd() gains an argument 'add'. o In the return values of regquantbounds() and sitequantbounds(), column names of "bound" items no longer contain an embedded space: their format is "bound.0.05" etc. o The data frame supplied as argument 'regdata' of regtst() now need have only 6 rather than 7 columns; if it does, some output quantities will not be calculated. (Originally claimed as a feature of Version 2.0, but not actually implemented therein.) Version 2.1 2009 Jul 10 o evplot.rfd() is now exported from the package namespace. o Fixed: as.regdata() would give an error when its argument 'x' was a data frame with default column names ("V1, "V2", etc.). Version 2.0 2009 Apr 30 o This version of the package aims to enable straightforward computation of all the procedures in chapters 3-6 of Hosking and Wallis (1997). To that end there are several new functions: * regavlmom() computes regional average L-moments. It is similar to reglmr() but its interface is more consistent with other functions in the package. * regfit() fits a regional frequency distribution to a regional data set. * regsimh() runs simulations to estimate the distribution of heterogeneity and goodness-of-fit measures. * regsimq() runs simulations to estimate the variability of quantile estimates from a regional frequency distribution. * regqfunc() and siteqfunc() return the regional growth curve and the quantile functions for individual sites, respectively, from a regional frequency distribution fitted by regfit(). * regquant() and sitequant() directly compute quantiles of the regional growth curve and of distributions for individual sites, respectively, from a regional frequency distribution fitted by regfit(). * regquantbounds() and sitequantbounds() compute error bounds for the regional growth curve and for quantiles at individual sites, respectively, from a regional frequency distribution fitted by regfit(). * print.rfd(), print.regsimh(), and print.regsimq() are print methods for the objects returned by regfit(), regsimh(), and regsimq(), respectively. o Class "regdata" is introduced as a standard structure for storing L-moments and other key data for the sites in a region. o New function as.regdata() coerces suitably structured data frames to class "regdata". o New function evplot.rfd() is a method for generic function evplot() in package 'lmom'. It facilitates graphical display of the results returned by regfit(), regquantbounds(), and sitequantbounds(). o regsamlmu(), with suitable values of its arguments, now returns an object of class "regdata". o Argument 'regdat' of regtst() is renamed 'regdata', and is expected to be of class "regdata" (after coercion if necessary). o reglmr() is deprecated; regavlmom() is the recommended replacement. Version 1.0 2009 Mar 1 o First public release. Functions included: cluagg cluinf clukm print.regtst print.summary.regtst reglmr regsamlmu regtst summary.regtst Data sets included: Appalach Cascades Maxwind