Title: | Duration Curve Hydrological Model Indexes |
---|---|
Description: | Compute duration curves of daily flow series, both real and modeled, to be compared through indexes of flow duration curves. The package functions include comparative plots and goodness of fit tests. Flow duration curve indexes are based on: Yilmaz et al., (2008) <DOI:10.1029/2007WR006716>. |
Authors: | Alonso Arriagada [aut, cre] |
Maintainer: | Alonso Arriagada <[email protected]> |
License: | GPL-3 |
Version: | 1.0.0 |
Built: | 2024-11-23 05:04:14 UTC |
Source: | https://github.com/alobondo/hydrodcindex |
Duration Curve Hydrological Model Indexes
hydroDC_Index(Q_obs, Q_sim, c_opt)
hydroDC_Index(Q_obs, Q_sim, c_opt)
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 |
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'
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
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)
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)
A column containing values of observed flow.
Q_obs
Q_obs
A object with 13514 rows and 1 variable:
observations, runoff in mm per day
A column containing values of simulated flow.
Q_sim
Q_sim
A object with 13514 rows and 1 variable:
simulations, runoff in mm per day