Computes confidence intervals for Flow-Normalized Concentration and Flow-Normalized Flux for a WRTDS model.
Usage
ciBands(eList, repAnnualResults, probs = c(0.05, 0.95))Arguments
- eList
named list with at least the Daily, Sample, and INFO dataframes. Created from the EGRET package, after running
modelEstimation.- repAnnualResults
named list returned from bootstrapping process.
- probs
numeric vector low and high confidence interval frequencies, default = c(0.05, 0.95) (which results in a 90% confidence interval).
Examples
library(EGRET)
eList <- Choptank_eList
nBoot <- 100
blockLength <- 200
if (FALSE) {
repAnnualResults <- vector(mode = "list", length = nBoot)
for(n in 1:nBoot){
annualResults <- bootAnnual(eList, blockLength, startSeed = n)
repAnnualResults[[n]] <- annualResults
}
CIAnnualResults <- ciBands(eList, repAnnualResults)
plotConcHistBoot(eList, CIAnnualResults)
}