combine | R Documentation |
Combines BMA models having the same characteristics for different dates.
combine( x, y, ...)
x |
An |
y |
An |
... |
Other |
Input models are checked for compatibility, and entries from different
inputs having the same dates are eliminated. Dates are ordered
chronologically and ensemble members are ordered in the order in which
they occur in inout x
.
An ensembleBMA
model that merges the models from each input into
a single model for the common dates.
C. Fraley, A. E. Raftery, T. Gneiting and J. M. Sloughter,
ensembleBMA
: An R
Package for Probabilistic Forecasting
using Ensemble and Bayesian Model Averaging,
Technical Report No. 516R, Department of Statistics, University of
Washington, 2007 (revised 2010).
ensembleBMA
data(ensBMAtest) ensMemNames <- c("gfs","cmcg","eta","gasp","jma","ngps","tcwb","ukmo") obs <- paste("T2","obs", sep = ".") ens <- paste("T2", ensMemNames, sep = ".") tempTestData <- ensembleData( forecasts = ensBMAtest[,ens], dates = ensBMAtest[,"vdate"], observations = ensBMAtest[,obs], station = ensBMAtest[,"station"], forecastHour = 48, initializationTime = "00") ## Not run: # R check tempTestFit12 <- ensembleBMAnormal( tempTestData, trainingDays = 30, dates = c("2008010100","2008010200")) tempTestFit34 <- ensembleBMAnormal( tempTestData, trainingDays = 30, dates = c("2008010300","2008010400")) ## End(Not run) # for quick run only; use more training days for forecasting tempTestFit12 <- ensembleBMAnormal( tempTestData, trainingDays = 8, dates = c("2008010100","2008010200")) tempTestFit34 <- ensembleBMAnormal( tempTestData, trainingDays = 8., dates = c("2008010300","2008010400")) tempTestFit <- combine( tempTestFit12, tempTestFit34)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.