Combine Compatible BMA Models
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).
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)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.