pyPDAF.PDAF.diag_compute_moments

pyPDAF.PDAF.diag_compute_moments(dim_p: int, dim_ens: int, ens: np.ndarray, kmax: int, bias: int) → np.ndarray

Computes the mean, the unbiased variance, the unbiased skewness, and the unbiased excess kurtosis from an ensemble.

Parameters:

  • dim_p (int) – local size of the state

  • dim_ens (int) – number of ensemble members/samples

  • ens (ndarray[np.float64, ndim=2]) – ensemble matrix Array shape: (dim_p,dim_ens)

  • kmax (int) – maximum order of central moment that is computed, maximum is 4

  • bias (int) – if 0 bias correction is applied (default) otherwise, no bias correction

Returns:

moments – The columns contain the moments of the ensemble - column 0: mean - column 1: variance - column 2: skewness - column 3: excess kurtosis Array shape: (dim_p, kmax)

Return type:

ndarray[np.float64, ndim=2]