Walks user through the set-up for the WRTDS Bootstrap Test. Establishes start and end year for the test period. Sets the minimum number of bootstrap replicates to be run, the maximum number of bootstrap replicates to be run, and the block length (in days) for the block bootstrapping. The test is designed to evaluate the uncertainty about the trend between any pair of years.
Arguments
- eList
named list with at least the Daily, Sample, and INFO dataframes. Created from the EGRET package, after running
modelEstimation.- ...
additional arguments to bring in to reduce interactive options (year1, year2, nBoot, bootBreak, blockLength)
Value
caseSetUp data frame with columns year1, yearData1, year2, yearData2, numSamples, nBoot, bootBreak, blockLength, confStop. These correspond to:
| Column Name | Manuscript Variable |
| year1 | \(y_s\) |
| year2 | \(y_e\) |
| nBoot | \(M_max\) |
| bootBreak | \(M_min\) |
| blockLength | \(B\) |
Examples
eList <- EGRET::Choptank_eList
# Completely interactive:
# caseSetUp <- trendSetUp(eList)
# Semi-interactive:
# caseSetUp <- trendSetUp(eList, nBoot = 100, blockLength = 200)
# fully scripted:
caseSetUp <- trendSetUp(eList,
year1=1985,
year2=2005,
nBoot = 50,
bootBreak = 39,
blockLength = 200)
#> Sample set runs from 1979 to 2011
#> year1 = 1985 this is the first water year of trend period
#> year2 = 2005 this is the last water year of trend period
#> nBoot = 50 this is the maximum number of replicates that will be run
#> bootBreak = 39 this is the minimum number of replicates that will be run
#> blockLength = 200 this is the number of days in a bootstrap block