Hello!
I'm having some problems with
pit function in ensembleBMA. When I run on mu data I have the
same warnings as below, when I type the example(pit).
Any suggestion on how to resolve it?
> library(ensembleBMA)
Loading required package: chron
> example(pit)
pit> data(ensBMAtest)
pit> ensMemNames <- c("gfs","cmcg","eta","gasp","jma","ngps","tcwb","ukmo")
pit> obs <- paste("T2","obs", sep = ".")
pit> ens <- paste("T2", ensMemNames, sep = ".")
pit> tempTestData <- ensembleData( forecasts = ensBMAtest[,ens],
pit+ dates = ensBMAtest[,"vdate"],
pit+ observations = ensBMAtest[,obs],
pit+ station = ensBMAtest[,"station"],
pit+ forecastHour = 48,
pit+ initializationTime = "00")
pit> ## Not run:
pit> ##D # R check
pit> ##D tempTestFit <- ensembleBMAnormal( tempTestData, trainingDays = 30)
pit> ## End(Not run)
pit>
pit> tempTestForc <- quantileForecast( tempTestFit, tempTestData)
pit> range(tempTestForc)
[1] 277.3138 281.1668
pit> tempTestPIT <- pit( tempTestFit, tempTestData)
Warning messages:
1: In matchDates(fitDates, ensembleValidDates(ensembleData), dates) :
some fit dates dates not represented in data
2: In matchDates(fitDates, ensembleValidDates(ensembleData), dates) :
some fit dates dates not represented in data
> sessionInfo()
R version 2.15.0 Patched (2012-04-10 r58965)
Platform: x86_64-unknown-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] ensembleBMA_5.0.3 chron_2.3-42
loaded via a namespace (and not attached):
[1] tools_2.15.0
Any suggestion?
Anna Freni Sterrantino
Ph.D Student
Department of Statistics
University of Bologna, Italy
via Belle Arti 41, 40124 BO.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
