Some functions useful to perform a Peak Over Threshold
        analysis in univariate and bivariate cases, see Beirlant et al. (2004)
        <doi:10.1002/0470012382>. A user guide is available in the vignette.

Provides functions to compute the potential model as defined by
    Stewart (1941) <doi:10.1126/science.93.2404.89>. Several options are available
    to customize the model, such as the possibility to fine-tune the distance
    friction functions or to use custom distance matrices. Some computations are
    parallelized to improve their efficiency.

Prepares groundwater-level observations from measured hydraulic
    head or from depth-to-water and land-surface elevation, interpolates
    potentiometric surfaces using thin-plate splines, inverse-distance
    weighting, ordinary Kriging, universal Kriging, or user-supplied methods,
    and creates raster, contour, diagnostic, support, and hydraulic-gradient
    products for review and export. Functions retain method conditions and fit
    diagnostics, validate explicit prediction tasks, inspect model-conditional
    uncertainty and monitoring-network sensitivity, identify limited prediction
    support, and check whether scaled hydraulic-gradient arrows remain within
    finite raster support and end at lower modeled head. Raster processing uses methods from
    'terra' (Hijmans 2025) <doi:10.32614/CRAN.package.terra>, thin-plate
    splines use methods from 'fields' (Nychka et al. 2021)
    <doi:10.5065/D6W957CT>, and geostatistical interpolation uses methods from
    'gstat' (Pebesma 2004) <doi:10.1016/j.cageo.2004.03.012>.

Calculate POTH for treatment hierarchies from frequentist and Bayesian network meta-analysis. POTH quantifies the certainty in a treatment hierarchy. Subset POTH, POTH residuals, and best k treatments POTH can also be calculated to improve interpretation of treatment hierarchies.

Store and retrieve data from options() using syntax derived from
    the 'here' package. 'potions' makes it straightforward to update and 
    retrieve options, either in the workspace or during package development, 
    without overwriting global options.