library(EGRET)
eList <- Choptank_eList
nBoot <- 100
blockLength <- 200
repAnnualResults <- vector(mode = "list", length = nBoot)
for(n = 1:nBoot){
annualResults <- bootAnnual(eList, blockLength)
repAnnualResults[[n]] <- bootAnnual(eList, blockLength)
}
CIAnnualResults <- ciBands(eList, repAnnualResults)Run the code above in your browser using DataLab