Combination of CST_SplitDim and CST_QuantileMapping
As reported by @silvia:
There is an issue when trying to use CST_SplitDim() in combination with the CST_QuantileMapping() function with seasonal forecast data at daily frequency.
I use the CST_SplitDim to split the time dimension in two, 'month' and 'day'. The dimension 'day' is of length 31 and it is filled with NA values for months with number of days < 31 .
The problem is that CST_QuantileMapping can not take directly the output of CST_SplitDim() because it can not manage NA values. The function fails with the following error:
exp.qm <- CST_QuantileMapping(expm,obsm, method = "QUANT", wet.day=F,sample_dims = c('member', 'sdate', 'ftime')) Error in QuantileMapping(exp = expdata, obs = obs
data, exp_cor = exp_cor$data, : Parameter 'exp' contains NA values.
I guess that we can not avoid producing NA when applying the SplitDim() function .. But is there a way to enable CST_QuantileMapping() to manage (or omit) NA values? Or to remove NAs from a CSTool object? In this way we could fully exploit the potential of the two functions and use them sequentially.