These plots use the jack-knife estimates from WRTDS to investigate the potential flux bias problem. It can also be used for estimates constructed by other methods (such as LOADEST) if the results are stored in a data frame organized like the Sample data frame. It allows additional label information to indicate what method is used. The use of this plot is described in Hirsch, Robert M., 2014. Large Biases in Regression-Based Constituent Flux Estimates: Causes and Diagnostic Tools. Journal of the American Water Resources Association (JAWRA) 1-24. DOI: 10.1111/jawr.12195
Although there are a lot of optional arguments to this function, most are set to a logical default.
Data come from named list, which contains a Sample dataframe with the sample data, a Daily dataframe with the daily flow data, and an INFO dataframe with metadata.
fluxBiasMulti(eList, qUnit = 2, fluxUnit = 3, moreTitle = "WRTDS",
cex = 0.7, cex.axis = 1.1, cex.main = 1.1, randomCensored = FALSE,
col = "black", lwd = 1, concLab = 1, monthLab = 1, ...)
named list with at least Sample, Daily, and INFO dataframes
object of qUnit class. printqUnitCheatSheet
, or numeric represented the short code, or character representing the descriptive name.
object of fluxUnit class. printFluxUnitCheatSheet
, or numeric represented the short code, or character representing the descriptive name.
character specifying some additional information to go in figure title, typically some information about the specific estimation method used, default is no additional information
numerical value giving the amount by which plotting symbols should be magnified
magnification to be used for axis annotation relative to the current setting of cex
magnification to be used for main titles relative to the current setting of cex
logical, if TRUE plot a random value for censored data. Default is FALSE.
color of points on plot, see ?par 'Color Specification'
number line width
object of concUnit class, or numeric represented the short code,
or character representing the descriptive name. By default, this argument sets
concentration labels to use either Concentration or Conc (for tiny plots). Units
are taken from the eList$INFO$param.units. To use any other words than
"Concentration" see vignette(topic = "units", package = "EGRET")
.
object of monthLabel class, or numeric represented the short code, or character representing the descriptive name.
arbitrary graphical parameters that will be passed to genericEGRETDotPlot function (see ?par for options)
eList <- Choptank_eList
# Water year:
fluxBiasMulti(eList)
fluxBiasMulti(eList, fluxUnit = 2)
# Graphs consisting of Jun-Aug
eList <- setPA(eList,paStart=6,paLong=3)
fluxBiasMulti(eList)
Run the code above in your browser using DataLab