Continuous rank probability score

Calculates continuous rank probability score (CRPS) for probabilistic and ensemble forecasts


Updated Thu, 11 Sep 2014 17:09:28 +0000

View License

The CRPS measures the closeness of forecast distribution (fcst) and corresponding observation (obs). This score is widely used in forecast verification.

[mean_CRPS] = crps(fcst,obs);
[mean_CRPS] = crps(fcst,obs,plot_pos);
[mean_CRPS,crps_values,num] = crps(fcst,obs);

obs: Vector of observations
fcst: Matrix of Ensemble forecast of size N x M. NB: N must equal length(obs), M equals the number of ensemble members
plot_pos: plotting positions that determine cumulative distribution function

mean_CRPS: Mean of non missing CRPS values
crps_values: A vector (length n) of CRPS values
num: number of non missing CRPS values used to compute mean_CRPS

fcst = rand(1000,1000);
obs = rand(1000,1);
[meanCRPS] = crps(fcst,obs);

Cite As

Durga Lal Shrestha (2023). Continuous rank probability score (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2014a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Version Published Release Notes