csep.core.poisson_evaluations.number_test

csep.core.poisson_evaluations.number_test(gridded_forecast, observed_catalog)

Computes “N-Test” on a gridded forecast. author: @asim

Computes Number (N) test for Observed and Forecasts. Both data sets are expected to be in terms of event counts. We find the Total number of events in Observed Catalog and Forecasted Catalogs. Which are then employed to compute the probablities of (i) At least no. of events (delta 1) (ii) At most no. of events (delta 2) assuming the poissonian distribution.

Parameters:

• observation – Observed (Gridded) seismicity (Numpy Array): An Observation has to be Number of Events in Each Bin It has to be a either zero or positive integer only (No Floating Point)

• forecast – Forecast of a Model (Gridded) (Numpy Array) A forecast has to be in terms of Average Number of Events in Each Bin It can be anything greater than zero

Returns:

(delta_1, delta_2)

Return type:

out (tuple)