selBestModel: Load Estimation

In USGS-R/rloadest: River Load Estimation

Description

Select the "best" predefined rating-curve (regression) model for river load estimation.

Usage

selBestModel(constituent, data, subset, na.action, flow, dates,
  flow.units = "cfs", conc.units = "", load.units = "kg",
  time.step = "day", station = "", criterion = c("AIC", "SPCC", "AICc"))

Arguments

constituent	a character string giving the name of the response variable for which loads are to be computed.
data	the data to search for the variables in formula.
subset	an expression to select a subset of the data.
na.action	what to do with missing values.
flow	character string indicating the name of the flow column.
dates	character string indicating the name of the date column.
flow.units	character string describing the flow unit.
conc.units	character string describing the concentration unit.
load.units	character string describing the load unit.
time.step	character string describing the time step of the calibration data.
station	character string description of the station.
criterion	the criterion to use for subset selection, must be one of "AIC," "SPCC," or "AICc."

Value

An object of class "loadReg."

References

will need some.

See Also

censReg

Examples

# Use the data from application 1 in the vignettes
data(app1.calib)
app1.lr <- selBestModel("Phosphorus", data = app1.calib, 
 flow = "FLOW", dates = "DATES", conc.units="mg/L",
 station="Illinois River at Marseilles, Ill.")
# Extract the fitted values
print(app1.lr)

USGS-R/rloadest documentation built on Oct. 2, 2020, 5:21 a.m.