New parameter for 'abs_thresholds' (!168) · Merge requests · Earth Sciences / s2dv

Merged Nuria Pérez-Zanón requested to merge dev-rpss-probs into master Oct 04, 2023

I have just pushed the changes I have for the GetProbs function. It now includes two new parameters: 'abs_thresholds' and 'prob_dim'.

Could you help me to test it whenever you have time?

It seems to work with the example:

source("/esarchive/scratch/nperez/git2/s2dv/R/GetProbs.R")
res <- GetProbs(data = data, time_dim = 'sdate', memb_dim = 'ensemble',
                 indices_for_quantiles = 4:17)

And you can also see my sample data in "/esarchive/scratch/nperez/sample_full_cross.RDS" to run the test:

source("/esarchive/scratch/nperez/git2/s2dv/R/GetProbs.R")
test <- GetProbs(ano_hcst_ev, time_dim = 'loop', prob_thresholds = NULL,
                 prob_dim = 'probs', indices_for_quantiles = NULL,
                 memb_dim = 'member', abs_thresholds = lims_ano_hcst_tr)

Could I keep adding other functions developments to this branch?

Thanks for your help!!

Núria

FYI @cdelgado