R/ensembleBMA.R

In ensembleBMA: Probabilistic Forecasting using Ensembles and Bayesian Model Averaging

`ensembleBMA` <-
function(ensembleData, trainingDays, dates = NULL, control = NULL, 
         model = NULL, exchangeable = NULL, minCRPS = NULL)
{
#
# copyright 2006-present, University of Washington. All rights reserved.
# for terms of use, see the LICENSE file
#
 if (!inherits(ensembleData,"ensembleData")) stop("not an ensembleData object")
 mc <- match.call()   
 mc$model <- NULL

## if (!is.null(model)) {
##   MODELS <- c("normal","gamma0")
##   m <- pmatch( model, MODELS, nomatch=0)
##   model <- if (m) MODELS[m] else "?"
##}
##else stop("unspecified model")

  switch( model,
        "normal" = {
             mc[[1]] <- as.name("ensembleBMAnormal")
          },
        "gamma0" = {
             mc[[1]] <- as.name("ensembleBMAgamma0")
             if (!is.null(minCRPS)) {
               if (minCRPS) warning("minCRPS not available for gamma0 model")
               mc$minCRPS <- NULL
             }
         },
        "gamma" = {
             mc[[1]] <- as.name("ensembleBMAgamma")
             if (!is.null(minCRPS)) {
               if (minCRPS) warning("minCRPS not available for gamma model")
               mc$minCRPS <- NULL
             }
         },
       stop("unrecognized model")
    )

 if (length(attr(ensembleData, "class")) > 2) {
   attr(ensembleData, "class") <- attr(ensembleData, "class")[-1]
   mc$ensembleData <- ensembleData
 }

 eval(mc, parent.frame())
}