Assesses the model's calibration quality with the help of the pairwise complete modelled as well as observed loads and the following metrics:

  • NSE: Nash-Sutcliffe Efficiency

  • mNSE: Modified Nash-Sutcliffe Efficiency (j = 1)

  • KGE: Modified Kling-Gupta Efficiency

  • RMSE: Root Mean Square Error

  • PBIAS: Percent Bias

  • RSR: Ratio of the RMSE to the standard deviation of the observations

  • RCV: Ratio of the coefficients of variation

  • GMRAE: Geometric Mean Relative Absolute Error

  • MdRAE: Median Relative Absolute Error

In addition, a scatter plot with the observed river loads on the x- and the modelled river loads on the y-axis is displayed and provides a visual impression of the model performance. Other elements of this plot are an identity line (solid) and plus/minus 30% deviation lines (dashed).

# S4 method for RPhosFate
calibrationQuality(x, substance, col)



An S4 RPhosFate river catchment object.


A character string specifying the substance to calculate.


A character string specifying the calibration data column with the respective substance river loads.


A named numeric vector containing the assessed metrics along with the in-channel retention ratio (one minus sum of xxt at catchment outlet(s) divided by sum of xxt_inp).


Nash, J.E., Sutcliffe, J.V., 1970. River flow forecasting through conceptual models part I – a discussion of principles. Journal of Hydrology 10, 282–290.

Legates, D.R., McCabe Jr., G.J., 1999. Evaluating the use of “goodness-of-fit” measures in hydrologic and hydroclimatic model validation. Water Resources Research 35, 233–241.

Kling, H., Fuchs, M., Paulin, M., 2012. Runoff conditions in the upper Danube basin under an ensemble of climate change scenarios. Journal of Hydrology 424–425, 264–277.

Moriasi, D.N., Arnold, J.G., Van Liew, M.W., Bingner, R.L., Harmel, R.D., Veith, T.L., 2007. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Transactions of the ASABE 50, 885–900.


# \donttest{
# temporary demonstration project copy
cv_dir <- demoProject()
#> Warning: A folder called "demoProject" already exists and is left as is.
# load temporary demonstration project
x <- RPhosFate(
  cv_dir = cv_dir,
  ls_ini = TRUE
# presupposed method calls
x <- firstRun(x, "SS")
x <- snapGauges(x)

calibrationQuality(x, "SS", "SS_load")# }
#> NSE:   0.7191828
#> mNSE:  0.5345233
#> KGE:   0.7323813
#> RMSE:  5.656556
#> PBIAS: -26.6
#> RSR:   0.4326794
#> RCV:   0.9761362
#> GMRAE: 0.3842057
#> MdRAE: 0.8072484
#> In-channel retention ratio: -2.220446e-16

#>                     NSE                    mNSE                     KGE 
#>            7.191828e-01            5.345233e-01            7.323813e-01 
#>                    RMSE                   PBIAS                     RSR 
#>            5.656556e+00           -2.660000e+01            4.326794e-01 
#>                     RCV                   GMRAE                   MdRAE 
#>            9.761362e-01            3.842057e-01            8.072484e-01 
#> inChannelRetentionRatio 
#>           -2.220446e-16