Plot reliability diagram

[Budusa]

I would like to use a reliability diagram for verification of my data organized in 16 ensemble members and one observation. Of course, all are in the netCDF format for June-August 2018. I was able to compute other scores such as PRS, CRPS, and others. provided all other skill scores and I would like to evaluate the results using the diagram. Thank you for your help in advance. Find also how I have tried to read and organize my data as follow:

require (ncdf4)
require(easyVerification)

fcstfiles <- list.files(pattern = 'wrf_t2m_')

fcst <- sapply(fcstfiles, function(f){

nc <- nc_open(f)

on.exit(nc_close(nc))

ncvar_get(nc, 't2m')

}, simplify = 'array')

library(cmsaf)
library(raster)

ncobs <- nc_open("era5_t2m_wrf_eth.nc")
obs <- ncvar_get(ncobs, 't2m')

corr <- veriApply("EnsCorr",
fcst = fcst,
obs = obs,
ensdim = 4, tdim = 3)

[rmanzanas]

Hi Budusa,

I see from your code that you're using the easyVerification package to compute skill scores others than reliability. Please note that whereas easyVerification requires matrix/array inputs, the reliabilityCategories function from visualizeR requires the specific data format of C4R.

Best regards.

[Budusa]

Hi Budusa,

I see from your code that you're using the easyVerification package to compute skill scores others than reliability. Please note that whereas easyVerification requires matrix/array inputs, the reliabilityCategories function from visualizeR requires the specific data format of C4R.

Best regards.

Hi rmanzamas,

Thank you for your comment. I will go through the material and come back to you if I fail to implement.

Best regards,