wrtdsK: WRTDS-Kalman

In EGRET: Exploration and Graphics for RivEr Trends

WRTDS-Kalman

Description

This function uses an autoregressive model to produce more accurate estimates of concentration and flux

Usage

WRTDSKalman(eList, rho = 0.9, niter = 200, seed = NA, verbose = TRUE)

Arguments

eList	named list with the INFO, Daily, and Sample dataframes and surfaces matrix
rho	numeric the lag one autocorrelation. Default is 0.9.
niter	number of iterations. Default is 200.
seed	integer value. Defaults to NA, which will not change the current seed. Setting the seed to any given value can be used to create repeatable output.
verbose	logical specifying whether or not to display progress message

Details

This function takes an existing eList Including the estimated model (the surfaces object in the eList) And produces the daily WRTDSKalman estimates of concentration and flux These generated estimates are called genConc and genFlux

Examples

eList <- Choptank_eList
eList <- WRTDSKalman(eList, niter = 10)
summary(eList$Daily)

#All flux values in AnnualResults are expressed as a rate in kg/day
AnnualResults <- setupYears(eList$Daily)
head(AnnualResults)

EGRET documentation built on April 18, 2023, 5:09 p.m.