From 1b12e3c30da405ee0637d9041c49ca3abefad5cc Mon Sep 17 00:00:00 2001 From: nperez Date: Mon, 16 Mar 2020 11:41:56 +0100 Subject: [PATCH 1/2] Calibratin exposed to user and example --- NAMESPACE | 1 + R/CST_Calibration.R | 26 +++++++++++++------------- man/CST_Calibration.Rd | 4 ++-- man/Calibration.Rd | 12 ++++++++++-- 4 files changed, 26 insertions(+), 17 deletions(-) diff --git a/NAMESPACE b/NAMESPACE index bd5d0f16..a319f567 100644 --- a/NAMESPACE +++ b/NAMESPACE @@ -21,6 +21,7 @@ export(CST_RFWeights) export(CST_RainFARM) export(CST_SaveExp) export(CST_SplitDim) +export(Calibration) export(EnsClustering) export(MergeDims) export(MultiEOF) diff --git a/R/CST_Calibration.R b/R/CST_Calibration.R index a7f1924b..ca290397 100644 --- a/R/CST_Calibration.R +++ b/R/CST_Calibration.R @@ -34,12 +34,9 @@ #'str(a) #'@export -CST_Calibration <- function(exp, obs, - cal.method = "mse_min", - eval.method = "leave-one-out", - multi.model = F, - na.fill = T, - ncores = 1) { +CST_Calibration <- function(exp, obs, cal.method = "mse_min", + eval.method = "leave-one-out", multi.model = FALSE, + na.fill = TRUE, ncores = 1) { if (!inherits(exp, "s2dv_cube") || !inherits(obs, "s2dv_cube")) { stop("Parameter 'exp' and 'obs' must be of the class 's2dv_cube', ", "as output by CSTools::CST_Load.") @@ -92,13 +89,16 @@ CST_Calibration <- function(exp, obs, #' #'@seealso \code{\link{CST_Load}} #' - -Calibration <- function(exp, obs, - cal.method = "mse_min", - eval.method = "leave-one-out", - multi.model = F, - na.fill = T, - ncores = 1) { +#'@examples +#'mod1 <- 1 : (1 * 3 * 4 * 5 * 6 * 7) +#'dim(mod1) <- c(dataset = 1, member = 3, sdate = 4, ftime = 5, lat = 6, lon = 7) +#'obs1 <- 1 : (1 * 1 * 4 * 5 * 6 * 7) +#'dim(obs1) <- c(dataset = 1, member = 1, sdate = 4, ftime = 5, lat = 6, lon = 7) +#'a <- Calibration(exp = mod1, obs = obs1) +#'str(a) +#'@export +Calibration <- function(exp, obs, cal.method = "mse_min", eval.method = "leave-one-out", + multi.model = FALSE, na.fill = TRUE, ncores = 1) { dim.exp <- dim(exp) amt.dims.exp <- length(dim.exp) diff --git a/man/CST_Calibration.Rd b/man/CST_Calibration.Rd index 891e2e5f..76812a43 100644 --- a/man/CST_Calibration.Rd +++ b/man/CST_Calibration.Rd @@ -9,8 +9,8 @@ CST_Calibration( obs, cal.method = "mse_min", eval.method = "leave-one-out", - multi.model = F, - na.fill = T, + multi.model = FALSE, + na.fill = TRUE, ncores = 1 ) } diff --git a/man/Calibration.Rd b/man/Calibration.Rd index 9f884671..64452279 100644 --- a/man/Calibration.Rd +++ b/man/Calibration.Rd @@ -9,8 +9,8 @@ Calibration( obs, cal.method = "mse_min", eval.method = "leave-one-out", - multi.model = F, - na.fill = T, + multi.model = FALSE, + na.fill = TRUE, ncores = 1 ) } @@ -37,6 +37,14 @@ Four types of member-by-member bias correction can be performed. The \code{bias} Both in-sample or our out-of-sample (leave-one-out cross validation) calibration are possible. } +\examples{ +mod1 <- 1 : (1 * 3 * 4 * 5 * 6 * 7) +dim(mod1) <- c(dataset = 1, member = 3, sdate = 4, ftime = 5, lat = 6, lon = 7) +obs1 <- 1 : (1 * 1 * 4 * 5 * 6 * 7) +dim(obs1) <- c(dataset = 1, member = 1, sdate = 4, ftime = 5, lat = 6, lon = 7) +a <- Calibration(exp = mod1, obs = obs1) +str(a) +} \references{ Doblas-Reyes F.J, Hagedorn R, Palmer T.N. The rationale behind the success of multi-model ensembles in seasonal forecasting-II calibration and combination. Tellus A. 2005;57:234-252. doi:10.1111/j.1600-0870.2005.00104.x -- GitLab From 5b73c14290d1f133912547bdab6946b001b3ba38 Mon Sep 17 00:00:00 2001 From: nperez Date: Fri, 29 May 2020 19:43:18 +0200 Subject: [PATCH 2/2] fix list in NEWS --- NEWS.md | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/NEWS.md b/NEWS.md index 04f986fb..3a846b1f 100644 --- a/NEWS.md +++ b/NEWS.md @@ -3,7 +3,8 @@ - Fixes + CST_Anomaly handles exp, obs or both - + PlotForecastPDF vignette displays figures correctly. + + PlotForecastPDF vignette displays figures correctly + + Calibration function is exposed to users ### CSTools 3.0.0 -- GitLab