Description Usage Arguments Details Value
Takes output from compile.vars
and determines
the best average model for each regime.
1 2 3 4 5 6 | best.models(top_n_list, nvmax, n, unfilled_FDCs, expl, comp.wys,
save.dir = NULL, cens_level = 0.005, zero.val = 0.001, WY.lim = 10,
regimes = list(lowflow = c("0.0002", "0.0005", "0.001", "0.002", "0.005",
"0.01", "0.02", "0.05", "0.1"), medflow = c("0.2", "0.25", "0.3", "0.4",
"0.5", "0.6", "0.7", "0.75", "0.8", "0.9"), highflow = c("0.95", "0.98",
"0.99", "0.995", "0.998", "0.999", "0.9995", "0.9998")))
|
top_n_list |
Output from |
nvmax |
The maximum number of variables that will be
considered in a regression. Prior to 9/2016, the default was |
n |
The maximum number of models subsetted.
Prior to 9/2016, the default was |
unfilled_FDCs |
A matrix of the raw FDC quantiles for each site.
This is derived from the output of |
expl |
A data frame of the potential explanatory variables,
derived from the output of |
comp.wys |
A vector indicating how many complete water years
were present for each site. This is derived from the output of
|
save.dir |
(optional) A directory to which results are to be saved.
The default behavior |
cens_level |
(optional) The value establishing left censoring.
The default is |
zero.val |
(optional) The value to which zeroes or negative quantiles
will be set to. The default is |
WY.lim |
(optional) The minimum number of water years required to be
included in the formation of regional regressions. Only reference-quality
sites with at least this many complete water years will be used.
The default is |
regimes |
(optional) A list of two or more elements ( |
This function outputs the top 3 models for each flow regime by both adjusted R^2 and AIC (for a total of 6 models). It takes the average adjusted R^2 and AIC across quantiles within a flow regime by each number of variables, and whichever number of variables has the highest R^2 or lowest AIC is the "best" type of model for the flow regime.
Top three models (by AIC and adjusted R2) for each regime.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.