ktindiana / sphinxval

SPHINX validation code for solar energetic particle models
MIT License
3 stars 3 forks source link

Deoverlapping for All Clear #93

Open ktindiana opened 5 months ago

ktindiana commented 5 months ago

Add deoverlapping techniques for All Clear validation.

For triggered forecasts:

As a first attempt for continuous forecasts:

rickyegeland commented 5 months ago

I like the approach for continuous forecasts, but would opt to make the quite time period bin width equal to the prediction window of the model – assuming the models are sane enough to keep that fixed. In any case, it would be good to be able to customize that bin width.

rickyegeland commented 5 months ago

Quiet time periods will generally not divide into integer N bin widths, so there will be 1 or 2 non-standard size bins for each quiet segment, depending on your choice for selecting the start time of the bin edges. I think this is OK given that N >> 2 most of the time.

rickyegeland commented 5 months ago

This approach should be extended to Probability forecasts as well, choosing to take the max probability in the quiet bins, as well as for SEP matches. The per-bin statistics (min, max, mean, std) of the model probability are also interesting quantities that indicate the model's robustness and self-consistency. A huge spread of probabilities in a bin indicates that the model produces a huge range of outputs that make it challenging for operators to trust any given single forecast.