R/ensembleData.R

In ensembleBMA: Probabilistic Forecasting using Ensembles and Bayesian Model Averaging

ensembleData <-
function(forecasts, dates = NULL, observations = NULL, ...,
         forecastHour, initializationTime,
         startupSpeed = NULL, exchangeable = NULL)
{
#
# copyright 2006-present, University of Washington. All rights reserved.
# for terms of use, see the LICENSE file
#
 object <- match.call()
 object[[1]] <- as.name("cbind.data.frame")
 object$forecastHour <- object$initializationTime <- object$startupSpeed <- object$exchangeable <- NULL
 object <- object[!unlist(lapply(object, is.null))]
 nam <- names(object)
 nam[2] <- nam[1]
 names(object) <- nam
 object <- eval(object, parent.frame())
 attr(object, "ensembleSize") <- if (is.null(dim(forecasts))) 1 else ncol(forecasts)
 if (missing(dates)) {
   warning("no dates specified")
 }
 else if (length(dates) == 1) {
   nObs <- if (is.null(dim(forecasts))) length(forecasts) else nrow(forecasts)
   dates <- as.factor(rep( as.character(dates), nObs))
 }
 else dates <- as.factor(dates)
 if (missing(forecastHour)) {
   warning("data cannot be used for modeling unless forecast hour is specified")
  }
 else attr(object, "forecastHour") <- forecastHour
 if (missing(initializationTime)) {
   warning("initialization time not specified")
 }
 else  attr(object, "initializationTime") <- initializationTime
 attr(object, "exchangeable") <- exchangeable
 attr(object, "startupSpeed") <- startupSpeed
 class(object) <- c("ensembleData", "data.frame")
 object
}