TonyLadson

#Packages

library(tidyverse)
library(stringr)
library(dplyr)
library(readxl)
library(readr)
library(purrr)
library(ggrepel)
library(hydroGOF)

TonyLadson

library(RCurl)
library(jsonlite)
library(stringr)
library(dplyr)
library(ggplot2)



rffe.data <- postForm("http://rffe.arr.org.au/",
                      catchment_name = "test1",

TonyLadson

# Figure 1
# DO - date of sample

WQ %>% 
  ggplot(aes(x = date, y = do.mg.l)) + geom_point() +
  labs(y = 'DO (mg/L)', x = 'Date') +
  theme_grey(base_size = 16)

# Figure 2
# DO - temperature

TonyLadson

library(ggplot2)
library(RColorBrewer)
library(grid)
library(cowplot)


ILb_hill <- function(d, MAR = 700){
  1 - 1/(1+142*sqrt(d)/MAR)
  }

TonyLadson

AEP2ARI <- function(AEP){
  1/-log(1-AEP)
}

ARI2AEP <- function(ARI){
  (exp(1/ARI) - 1)/exp(1/ARI)
}

EY2AEP <- function(EY){
  ARI2AEP(1/EY)

TonyLadson

library(tidyverse)
library(stringr)
library(RColorBrewer)
library(scales)
library(lubridate)
library(readxl)


# Avon River at Stratford, gauge 225201
# Stream gauging and rating data is available at http://www.bom.gov.au/waterdata/

TonyLadson

FDC <- function(Q, xnormal=TRUE, xlab='%time flow equalled or exceeded',
                ylab='flow',Qreturn=FALSE, plotFlow=FALSE, add=FALSE, alpha=0.5, ...) {
  
  # Q - flow data (vector of flow values)
  # xnormal - logical value to determine if the percentiles should be plotted in a normal (probit) scale
  # Qreturn - logical value indicating if quantiles should be returned
  # plotFlow - logical value indicating if a time series of flow should be overplotted on flow duration curve
  # add - should fdc be added to an existing plot?
  # alpha - transparency value for overplotting of flow time series
  

TonyLadson

#
# This is a Shiny web application. You can run the application by clicking
# the 'Run App' button above.
#
# Find out more about building applications with Shiny here:
#
#    http://shiny.rstudio.com/
#

library(shiny)

TonyLadson

---
title: "Events_blog"
author: "Tony Ladson"
date: "3 October 2016"
output:
  html_fragment:
    fig_caption: yes
    fig_height: 4
    fig_width: 6
---

TonyLadson

library(stringr)
library(tidyverse)

my.path <- c("...")
my.file <- '404216.rsrspf_org.tab.prn 2.txt'
fname <- str_c(my.path, my.file, sep = '/')

#  Read the file in as a character vector 
rating_table <- read_lines(fname)