Skip to content

Instantly share code, notes, and snippets.

View TonyLadson's full-sized avatar

Tony Ladson TonyLadson

20 followers · 2 following
View GitHub Profile
@TonyLadson
TonyLadson / Calc_slope_ea.R
Last active August 16, 2021 17:04
Formula to calculate the equal area slope. See https://tonyladson.wordpress.com/2017/03/04/on-the-calculation-of-equal-area-slope/
# Description
# Calc_slope_ea(x,y) calculates teh equal area slope
# Arguments
#
# x = distance along the stream in km
# y = elevation along the stream in m
# Details
@TonyLadson
TonyLadson / PercentilePlot.R
Last active June 16, 2021 20:11
Plotting percentiles see: https://tonyladson.wordpress.com/2016/02/22/plotting-percentiles/
library(quantreg)
library(dplyr)
library(lubridate)
library(stringr)
library(readr)
library(ggplot2)
library(devtools)
library(splines)
@TonyLadson
TonyLadson / gist:286077221cc07022a77aa74d0a6c211c
Created October 4, 2020 23:49
Regression examples and figures for the blog https://tonyladson.wordpress.com/2020/10/05/errors-in-variables-regression/
library(tidyverse)
library(deming)
# Functions ---------------------------------------------------------------
# given a point (x1, y1) and a line defined by a slope m and intercept c1
# function to return the point x on the line where a line drawn perpendicular to x1, y1
# will intercept the line
@TonyLadson
TonyLadson / Flood_sampling_distribution.R
Last active June 2, 2020 11:10
Sampling distribution of floods. See https://tonyladson.wordpress.com/2020/06/02/sampling-distribution-of-the-1-flood/
library(tidyverse)
#pbinom(q, size, prob, lower.tail = TRUE, log.p = FALSE)
# Probability of 1, 1% flood in 100 years
dbinom(1, size = 100, prob = 0.01, log = FALSE)
choose(100,1) *0.01^1*0.99^99
@TonyLadson
TonyLadson / Issues_column_graph.R
Created June 1, 2020 01:08
Code to draw the column graph in the blog at https://tonyladson.wordpress.com/2020/04/29/comments-on-the-report-impact-of-lower-inflows-on-state-shares-under-the-murray-darling-basin-agreement/
library(tidyverse)
issues <- tribble(~issue, ~pc,
'Inflows', 6,
'Demand changes', 6,
'System operation', 18,
'Equity of water sharing arrangements', 10,
'State entitlements and policies', 10,
'Environmental water', 11,
'Reliability of different licences', 4,
'Miscellaneous', 4,
@TonyLadson
TonyLadson / ARF_poly_interp.R
Last active May 23, 2020 12:14
Polynomial interpolation of ARF values between 12 and 24 hours. See https://tonyladson.wordpress.com/2020/05/23/smooth-interpolation-of-arf-curves/
# Smooth interpolation between long and short duration ARFs using a cubic.
library(tidyverse)
library(pracma)
devtools::source_url("https://gist.githubusercontent.com/TonyLadson/fc870cf7ebfe39ea3d1a812bcc53c8fb/raw/d8112631a92a32be749cabe334a22931c035711e/ARF2019.R?raw=TRUE")
@TonyLadson
TonyLadson / ARF_map.R
Created May 12, 2020 00:16
Code to produce Figure 2 for the Blog https://tonyladson.wordpress.com/2020/04/05/arr2019-areal-reduction-factors/
library(tidyverse)
xdf <- tribble(~area_min, ~area_max, ~dur_min, ~dur_max, ~type,
10, 1000, 0, 12, 'short',
10, 30000, 24, 168, 'long',
10, 30000, 12, 24, 'interpolation',
0, 1, 0, 168, 'ARF = 1',
1, 10, 0, 12, 'Interp3',
1, 10, 12, 24, 'Interp2',
1, 10, 24, 168, 'Interp1')
@TonyLadson
TonyLadson / ARF_interp.R
Created May 12, 2020 00:03
Code to reproduce the figures in the blog: ARR2019 – Areal Reduction Factors: interpolating between short and long duration ARFs
# Load the functions we need
devtools::source_url("https://gist.githubusercontent.com/TonyLadson/fc870cf7ebfe39ea3d1a812bcc53c8fb/raw/d8112631a92a32be749cabe334a22931c035711e/ARF2019.R?raw=TRUE")
#source(file.path('ARR2019_ARF', "ARF_2019.R"))
#source(file.path('ARR2019_ARF', "ARF_tests.R")) # Check that we pass tests
# Functions and constants ------------------------------------------------------
@TonyLadson
TonyLadson / ARF_sd.R
Created April 19, 2020 23:40
Investigating short duration Areal Reduction Factors https://tonyladson.wordpress.com/2020/04/19/arr2019-areal-reduction-factors-wobbles-in-short-duration-arfs/
library(tidyverse)
library(optimx)
devtools::source_url("https://gist.githubusercontent.com/TonyLadson/fc870cf7ebfe39ea3d1a812bcc53c8fb/raw/d8112631a92a32be749cabe334a22931c035711e/ARF2019.R?raw=TRUE")
#source(file.path('ARR2019_ARF', "ARF_2019.R"))
#source(file.path('ARR2019_ARF', "ARF_tests.R")) # Check that we pass tests
@TonyLadson
TonyLadson / ARF_edge_cases.R
Last active April 17, 2020 00:35
Areal reduction factors - some edge cases https://tonyladson.wordpress.com/2020/04/14/arr2019-areal-reduction-factors-some-edge-cases/
library(tidyverse)
devtools::source_url("https://gist.githubusercontent.com/TonyLadson/fc870cf7ebfe39ea3d1a812bcc53c8fb/raw/d8112631a92a32be749cabe334a22931c035711e/ARF2019.R?raw=TRUE")
#source(file.path('ARR2019_ARF', "ARF_2019.R"))
#source(file.path('ARR2019_ARF', "ARF_tests.R"))
# Functions and data ------------------------------------------------------