makeMonthPI {EGRETci} | R Documentation |
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 |
eList |
named list with at least the Daily, Sample, and INFO
dataframes. Created from the EGRET package, after running |
fluxUnit |
number representing entry in pre-defined fluxUnit class array.
|
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"]])