historical_daily_range {SWMPrExtension} | R Documentation |
Historical Daily Range Timeseries
Description
Compare daily averages for a target year to historical highs and lows
Usage
historical_daily_range(swmpr_in, ...)
## S3 method for class 'swmpr'
historical_daily_range(
swmpr_in,
param = NULL,
hist_rng = NULL,
target_yr = NULL,
criteria = NULL,
free_y = FALSE,
log_trans = FALSE,
converted = FALSE,
criteria_lab = "WQ Threshold",
plot_title = FALSE,
plot = TRUE,
...
)
Arguments
swmpr_in |
input swmpr object |
... |
not used |
param |
chr string of variable to plot |
hist_rng |
numeric vector, if historic range is not specified then the min/max values of the data set will be used. |
target_yr |
numeric, the target year that should be compared against the historic range. If target year is not specified then dot will not be plotted |
criteria |
numeric, a numeric criteria that will be plotted as a horizontal line |
free_y |
logical, should the y-axis be free? Defaults to |
log_trans |
logical, should y-axis be log? Defaults to |
converted |
logical, were the units converted from the original units used by CDMO? Defaults to |
criteria_lab |
chr, label for the threshold criteria defined in |
plot_title |
logical, should the station name be included as the plot title? Defaults to |
plot |
logical, should a plot be returned? Defaults to |
Details
This function compares the average daily minimums/maximums and absolute daily minimums/maximums from a historical range to the average daily value from a target year. If hist_rng
is not specified then the minimum and maximum years within the data set will be used. If target_yr
is not specified then only the results for the hist_rng
will be returned.
The user also has the option to add a threshold line using the criteria
argument. Typically, this value is a water quality threshold, which is why criteria_lab
defaults to 'WQ Threshold'
. However, the user has the option to specify any other type of threshold they wish. when doing so, the value for criteria_lab
should be changed accordingly.
Value
Returns a ggplot
object
Author(s)
Julie Padilla, Kimberly Cressman
See Also
Examples
data(apacpwq)
dat <- apacpwq
dat <- qaqc(apacpwq, qaqc_keep = c('0', '3', '5'))
# with criteria
y <- historical_daily_range(dat, param = 'do_mgl', target_yr = 2013, criteria = 2)
# w/o criteria
z <- historical_daily_range(dat, param = 'do_mgl', target_yr = 2013)
# add a y label
zz <- z + labs(x = NULL, y = "Dissolved Oxygen (mg/L)")