R/sbi_statistics.R

In bcgov/bcsnowsbi: Package for Calculating Snow Basin Indices (SBI) Values for British Columbia Subbasins

#
# File for calculating statistics for SBI

# Calculate the return period based on the past SBI for the basin.
# Calculate for a specifed SBI value
# As a function that returns the return period and probability for a specific SBI percent.

# Calculate the probability for the SBI based on the past SBI values for that region
# Started 17Jan2019
# Ashlee Jollymore, RFC
# =============

# ===============================================================================
# Packages
# ===============================================================================
#library("RCurl")
#library("dplyr")
#library("tidyr")
#library("data.table")
#library("lubridate")
#library(gsubfn)
#library(RCurl)
#library(lubridate)
#library(tibble)
#library('kableExtra')
#library('knitr')
#library(dplyr)
#library(tinytex)
#options(knitr.table.format = "html")
#library(rmarkdown)
#library(lubridate)
#library(memoise)
#library(stringr)
#library(ggplot2)
#library(arsenal)
#library(magrittr)
#library(readr)

# ===================
# Load the necessary functions
# Should eventually just be replace with a bcsnow library call
# ===================
#source('R/SWE_getdata_function.R')
#source('R/snow.R')
#source('R/utils.R')
#source('R/SWE_function_wateryear.R')
#source('R/SWE_statisticalanalysis.R')
#source('R/SWE_statistics_function.R')
#source('SWE_statisticstable_function.R')
#source('R/SWE_getstationID_function.R')
#source('SWE_SBI_function.R')

# Gather the previously calculated SBI values

# Knit them together within an entire large dataset that you can then subset by the basin name

# Are SBI values normally distributed? Return distribution for the SBI for a specific basin

# calculate the return period for a specified value based on the past SBI values.

# Calculate the probability of a specific value based on past SBI values