From d06bf1b2a68dbec32c9d50838eec81af01a9187c Mon Sep 17 00:00:00 2001 From: Carlos Delgado Torres Date: Wed, 28 Jul 2021 12:54:13 +0000 Subject: [PATCH] Changed s2dv::MeanDims by multiApply(fun = mean) for not losing the dimnames --- R/CST_Calibration.R | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/R/CST_Calibration.R b/R/CST_Calibration.R index b58dd08e..2884712a 100644 --- a/R/CST_Calibration.R +++ b/R/CST_Calibration.R @@ -300,8 +300,8 @@ Calibration <- function(exp, obs, cal.method = "mse_min", #correct evaluation subset var.cor.fc[ , eval.dexes] <- .correct.crps.min.fc(fc.ev , mbm.par, na.rm = na.rm) } else if (cal.method == 'rpc-based') { - ens_mean.ev <- s2dv::MeanDims(data = fc.ev, dims = names(amt.mbr), na.rm = na.rm) - ens_mean.tr <- s2dv::MeanDims(data = fc.tr, dims = names(amt.mbr), na.rm = na.rm) ## Ensemble mean + ens_mean.ev <- multiApply::Apply(data = fc.ev, target_dims = names(amt.mbr), fun = mean, na.rm = na.rm)$output1 ## Ensemble mean + ens_mean.tr <- multiApply::Apply(data = fc.tr, target_dims = names(amt.mbr), fun = mean, na.rm = na.rm)$output1 ## Ensemble mean ens_spread.tr <- multiApply::Apply(data = list(fc.tr, ens_mean.tr), target_dims = names(amt.sdate), fun = "-")$output1 ## Ensemble spread exp_mean.tr <- mean(fc.tr, na.rm = na.rm) ## Mean (climatology) var_signal.tr <- var(ens_mean.tr, na.rm = na.rm) ## Ensemble mean variance -- GitLab