EGRET (version 3.0.2)

modelEstimation: Estimation process for the WRTDS (Weighted Regressions on Time, Discharge, and Season)

Description

This one function does three things. 1) a jack-knife cross-validation of a WRTDS model in which it augments the Sample data frame in the eList 2) fits the WRTDS model creating the, fits the surfaces matrix and places it in the eList (the surfaces matrix expresses the estimated concentration as a function of discharge and time), 3) estimates the daily values of concentration and flux, and flow normalized concentration and flux and places these in the Daily data frame in the eList values. It returns a named list with the following dataframes: Daily, INFO, Sample, and the matrix: surfaces.

Usage

modelEstimation(eList, windowY = 7, windowQ = 2, windowS = 0.5,
  minNumObs = 100, minNumUncen = 50, edgeAdjust = TRUE,
  verbose = TRUE, run.parallel = FALSE)

Arguments

eList

named list with at least the INFO, Daily, and Sample dataframes

windowY

numeric specifying the half-window width in the time dimension, in units of years, default is 7

windowQ

numeric specifying the half-window width in the discharge dimension, units are natural log units, default is 2

windowS

numeric specifying the half-window with in the seasonal dimension, in units of years, default is 0.5

minNumObs

numeric specifying the miniumum number of observations required to run the weighted regression, default is 100

minNumUncen

numeric specifying the minimum number of uncensored observations to run the weighted regression, default is 50

edgeAdjust

logical specifying whether to use the modified method for calculating the windows at the edge of the record. The edgeAdjust method tends to reduce curvature near the start and end of record. Default is TRUE.

verbose

logical specifying whether or not to display progress message

run.parallel

logical to run WRTDS in parallel or not

Value

eList named list with INFO, Daily, and Sample dataframes, along with the surfaces matrix. Any of these values can be NA, not all EGRET functions will work with missing parts of the named list eList.

Examples

# NOT RUN {
eList <- Choptank_eList
# }
# NOT RUN {
 
#Run an estimation adjusting windowQ from default:
eList <- modelEstimation(eList, windowQ=5)

library(doParallel)
nCores <- parallel::detectCores() - 1
cl <- makePSOCKcluster(nCores)
registerDoParallel(cl)
eList <- modelEstimation(eList, windowQ=5, run.parallel = TRUE)
stopCluster(cl)
# }