Changes in version 1.3.1 (2026-02-27)                  

Bug Fixes

  - Replaced all \dontrun{} with unwrapped examples or \donttest{} per
    CRAN policy.
  - Wrapped Trehafod and glme.gev11() examples in \donttest{} (> 5 sec
    on Debian).
  - Fixed invalid URL for UK National River Flow Archive Peak Flow
    Dataset.
  - Single-quoted acronyms and proper names in DESCRIPTION to avoid
    spell-check notes.

                        Changes in version 1.3.0                        

Breaking Changes

  - BREAKING: glme.gev() output names changed for consistency:
      - lme → para.lme
      - glme → para.glme
      - nllh.pref → nllh.glme
      - covinv → covinv.lmom
  - BREAKING: ma.gev() output names changed for consistency:
      - zp.ma → qua.ma
      - zp.bma → qua.bma
      - fin.se.ma → fixw.se.ma
      - adj.se.ma → ranw.se.ma
      - numk_ma and numk_bma → run.numk
      - pick_xi_ma and pick_xi_bma → pick_xi
      - remle1 → para.remle1 (in return and internal use)
      - remle2 → para.remle2 (in return and internal use)

Internal Function Renames

  - Renamed internal functions for consistency (these are not exported
    but may affect code using ::: accessor):
      - init.glme() → init.gevmax()
      - new_pf_norm() → pk.norm.stnary() (with backward-compatible
        alias)
      - gev.rl.delta_new() → gev.rl.delta()
      - lme.boots.new() → lme.boots()
      - cand.xi.new.paper() → cand.xi()
      - weight.com.new() → weight.com()
      - cov.interp.new() → cov.interp()
      - gev.profxi.mdfy.paper() → gev.profxi.mdfy()
      - comp.prof.ci.new() → comp.prof.ci()
      - gev1.CD() → mle.gev.CD()
      - gev.remle() → remle.gev()
      - ginit.max() → init.gevmax()
      - pargev.xifix.ma() → pargev.xifix()
  - Penalty functions in set.prior() now use pk.beta.stnary() from
    glme.gev.R

                        Changes in version 1.2.0                        

Breaking Changes

  - BREAKING: glme.gev11() output para.jkss renamed to para.lme for
    consistency. Users accessing result$para.jkss should update their
    code to use result$para.lme.
  - BREAKING: glme.gev11() and gado.prop_11() output strup.final renamed
    to para.wls. Users accessing result$strup.final should update their
    code to use result$para.wls.
  - BREAKING: glme.gev11() no longer returns strup.sta in its output.

New Features

  - Enhanced glme.gev11() with new parameters:
      - glme.pre = "wls": Pre-estimation method selection ("wls" or
        "gado")
      - choose = "gof": Model selection criterion ("gof" for
        goodness-of-fit, "nllh" for negative log-likelihood)
      - pen.choice = 6: Default penalty hyperparameter choice changed
        from NULL to 6
  - New quagev.NS() function for calculating quantiles from
    non-stationary GEV models
      - Supports GEV11, GEV10, GEV20, and stationary GEV00 models
      - Returns time-varying quantiles as vector or matrix
  - Enhanced ma.gev() with new estimation options:
      - CD = TRUE: Coles-Dixon penalized MLE for shape parameter
        regularization
      - remle = TRUE: Restricted MLE with mean/median constraints
      - Returns mle.CD, qua.CD, remle1, remle2, qua.remle1, qua.remle2
      - Returns quant in output for convenience
      - BMA outputs now include bma.se.between and bma.se.within
  - New diagnostic plotting functions for MAGEV:
      - magev.ksensplot(): K sensitivity analysis to select optimal
        number of submodels
      - magev.qqplot(): 2x2 Q-Q diagnostic plot comparing MLE, LME,
        surrogate, and REMLE
      - magev.rlplot(): Return level plot with 95% confidence intervals
  - Added bangkok dataset: Annual maximum daily rainfall from Bangkok,
    Thailand
  - Added haenam dataset: Annual maximum daily rainfall from Haenam,
    South Korea

Bug Fixes

  - Improved handling of single quantile (numq = 1) in ma.gev()

                        Changes in version 1.1.0                        

  - Added Model Averaging GEV estimation (ma.gev()) for high quantile
    estimation.
      - Combines MLE and L-moment estimates through weighted model
        averaging
      - Multiple weighting schemes: like, gLd, med, cvt and variants
      - Optional Bayesian model averaging (bma=TRUE) with normal/beta
        priors
      - Returns model-averaged quantiles (qua.ma) with standard errors
  - New dependencies: ismev, Rsolnp, zoo.
  - Reference: Shin et al. (2026) SERA, 40(2), 47 (MAGEV method)

                        Changes in version 1.0.0                        

  - Initial CRAN release.
  - Stationary GEV parameter estimation using generalized L-moments
    (glme.gev()).
  - Non-stationary GEV11 model estimation (glme.gev11()) where location
    (mu) and scale (sigma) parameters vary linearly with time.
      - para.glme: Proposed GLME estimates
      - para.lme: L-moment based estimates for non-stationary model
  - Compatibility functions for Shin et al. (2025, J. Korean Stat.
    Soc.):
      - nsgev(): Simple interface for L-moment based non-stationary
        estimation
      - gado.prop_11(): Comprehensive estimation with multiple methods
  - Multiple penalty function options: "beta" (default), "norm", "ms"
    (Martins-Stedinger), "park", "cannon", "cd" (Coles-Dixon), and "no"
    (no penalty).
  - Flexible hyperparameter specification via pen.choice or direct
    parameters (p, c1, c2 for beta; mu, std for normal penalty).
  - Included datasets: streamflow, PhliuAgromet, and Trehafod.
  - References:
      - Shin et al. (2025a) arXiv:2512.20385 (GLME method)
      - Shin et al. (2025b) J. Korean Stat. Soc. 54:947-970
        (Non-stationary L-moment)