## R code - Handy routines for hydrologists

### What is R?

R is an open source programming language and environment for data analysis. It has rich functionality for data processing, analysis and graphing. See the R home page for more information and to download the package. Be forewarned that R has a steep learning curve. However, once you gain proficiency in R, the time invested in learning R will be paid back many times over.

### R and hydrology

Data analysis is an integral part of hydrology. Hydrologists frequently use techniques, such as regression analysis, which are incorporated into conventional statistical packages and spreadsheet software. However, many hydrological analyses are not, including intensity-duration-frequency analysis and flood frequency analysis. These analyses are relatively simple to code in R.

### CSHShydRology

CSHShydRology is an R package of functions developed by Canadian hydrologists. The name is in recognition of the support provided by the Canadian Society for Hydrological Sciences (CSHS), which is an affiliated society of the Canadian Water Resources Association (CWRA). The CSHS website is http://www.cwra.org/en/branches/affiliates/cshs-a.

Many of my sample scripts will be superceded by functions within CSHShydRology.

### R primers and background material

In the following essay, Dave and I outline the advantages of using R for hydrologic analyses, particularly the fact you can generate integrated workflows that encompass all steps, from initial data cleaning to analysis to reporting. The following documents are based on materials that I developed for a graduate course on data analysis and modelling in hydrology and geomorphology. The first two (basics and programming) exclusively use base R. Later installments will incorporate functions from contributed packages including the `tidyverse` set of packages. The documents were generated using R Markdown.- R basics for hydrologists [pdf version] [web version] (2018-June-03)
- R programming for hydrologists [pdf version] [web version] (2018-June-03; updated 2018-Sept-12)

### R scripts for specific hydrologic analyses

What I have included here are bits of R code that hydrologists may find useful. If you have written some code that you would like to share, please send it along and I'll consider posting it here.

You are free to use, distribute and modify the code below. If you do so, please send me an email to let me know which scripts you have used.

**Disclaimer**: All of the code presented on this page is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. The code is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

**Note**: I am an old dog who generally avoids learning new tricks. I began my programming career with WATFIV and punch cards and then progressed through FORTRAN 77, Turbo Pascal and Fortran 90/95 before discovering R, which satisfies all of my current programming/data analysis requirements. Most of my older code uses base R, with more recent code making increasing use of *tidyverse* packages, especially *dplyr*. Although I use ggplot graphics occasionally for plots that involve grouping or conditioning, I tend to use base graphics for complex plots that require subtantial customization.

- Catchment delineation with RSAGA and mapping using Spatial* and raster objects in R (2016-Aug-20)
- Script (based in part on a script by Joel Trubilowicz, Headwater Analytics)
- DEM and shape files (zipped)

- Flood frequency and rainfall intensity-duration frequency analysis
- Intensity-duration-frequency graph paper
- Flood frequency plotting on extreme value graph paper
- Intensity-duration-frequency analysis for Nanaimo Airport, Vancouver Island (NB: needs Gumplot-idf.r; see next script)
- Function to generate Extreme Value I (Gumbel) plot (called by Nanaimo-idf.r)

- Various older scripts
- Fitting a power-law relation between stream discharge and stage
- Binomial filter for regularly sampled time series
- Baseflow separation using a recursive low-pass filter (Nathan and McMahon, 1991, WRR)
- Use of linear regression to fill in missing air temperature data, with separate regressions by month