hydroDC_Index {hydroDCindex}

Duration Curve Hydrological Model Indexes

Description

Duration Curve Hydrological Model Indexes

Usage

hydroDC_Index(Q_obs, Q_sim, c_opt)

Arguments

Q_obs

Column with daily observed flows

Q_sim

Column with daily simulated flows

c_opt

Results, default 1 for indexes, 2 for duration curve values, 3 for DC plot and 4 for scatter plot

Details

If c_opt=1, it computes the numerical values of:
'BiasFMS','BiasFHV','BiasFLV','BiasFMM_log','BiasFMM' If c_opt=2, it computes the duration curve values:
'Pexc','Amount_Obs','Amount_Sim' If c_opt=3, it plots the duration curves. If c_opt=4, it plots scattered values and computes the numerical values of:
'r_pearson','MAE','rsq','NSE','KGE'

Value

BiasFMS: Diagnosis of vertical redistribution in the midsection of the duration curve.

BiasFHV: Bias in peak flows.

BiasFLV: Bias at low flows.

BiasFMM_log: Log_Mean Flow Bias.

BiasFMM: Mean Flow Bias.

Pexc: Exceedance probability computed with Weibull formula.

Amount_Obs: Observations in ascending order.

Amount_Sim: Simulations in ascending order.

r_pearson: Imported method from stats package, results in the Pearson product-moment correlation coefficient, -1 <= r_pearson <= 1

MAE: Imported method from ie2misc package, results in the Mean Absolute Error.

rsq: Results in the Coefficient of Determination, estimated as the squared of correlation, 0 <= R2 <= 1

NSE: Imported method from ie2misc package, results in the Nash-Sutcliffe Efficiency, -Inf <= NSE <= 1

KGE: Results in the Kling-Gupta Efficiency, 0 <= KGE <= 1

Examples

  Obs <- hydroDCindex::Q_obs
  Sim <- hydroDCindex::Q_sim
  # option 1 for indexes
  hydroDCindex::hydroDC_Index(Obs, Sim, 1)
  # option 2 for duration curve values
  hydroDCindex::hydroDC_Index(Obs, Sim, 2)
  # option 3 for duration curve plot
  hydroDCindex::hydroDC_Index(Obs, Sim, 3)
  # option 4 for scatter plot with goodness of fit tests
  hydroDCindex::hydroDC_Index(Obs, Sim, 4)