system.file("extdata","run1_Hydrographs.csv", package="RavenR") %>%
rvn_hyd_read(.) %>%
rvn_hyd_extract(subs="Sub36",.) ->
hyd_data
# Calculate quantiles for the simulated hydrograph
qdat <- rvn_annual_quantiles(hyd_data$sim)
head(qdat)
# Plot
p <- rvn_annual_quantiles_plot(qdat)
p # view plot
# Add a second hydrograph to compare
qdat_sim <- rvn_annual_quantiles(hyd_data$sim)
p1 <- rvn_annual_quantiles_plot(qdat_sim, mediancolor = 'blue', ribboncolor = 'red', explot = p)
p1 # view plot
Run the code above in your browser using DataLab