makeMonthPI: Make Monthly Prediction Intervals

In USGS-R/EGRETci: Exploration and Graphics for RivEr Trends Confidence Intervals

Make Monthly Prediction Intervals

Description

Month statistics using WRTDSKalman bootstrapping approach. The input to this function is the dailyBootOut matrix which contains nReplicate sets of daily flux values for the period of interest. The results are in the form of quantiles of concentration and of flux for each of these months.

Usage

makeMonthPI(dailyBootOut, eList, fluxUnit = 3)

Arguments

dailyBootOut	data frame returned from genDailyBoot
eList	named list with at least the Daily, Sample, and INFO dataframes. Created from the EGRET package, after running modelEstimation.
fluxUnit	number representing entry in pre-defined fluxUnit class array. printFluxUnitCheatSheet

Value

a list of 2 data frames, one for average concentration, in mg/L and one for flux (unit depends on fluxUnit argument) In each data frame the first column is monthSeq that corresponds to the months in the "MonthSeq" column in the eList$Daily data frame. The remaining columns are quantiles of the flux or concentration (depending on the data frame).

Examples

eList <- EGRET::Choptank_eList
# This example is only based on 4 iterations
# Actual prediction intervals should be calculated on
# a much larger number of iterations (several hundred). 
dailyBoot <- Choptank_dailyBootOut
monthPcts <- makeMonthPI(dailyBoot, eList)
head(monthPcts[["flux"]])
head(monthPcts[["conc"]])

USGS-R/EGRETci documentation built on April 17, 2023, 8 p.m.