# NOT RUN {
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 checl
tempTestFit <- ensembleBMA( tempTestData, trainingDays = 30,
model = "normal")
## equivalent to
## tempTestFit <- ensembleBMAnormal( tempTestData, trainingDays = 30)
# }
# NOT RUN {
# for quick run only; use more training days for forecasting
tempTestFit <- ensembleBMA( tempTestData[1:20,], trainingDays = 8,
model = "normal")
# }
Run the code above in your browser using DataCamp Workspace