11.13.2. astroML.resample.jackknife

astroML.resample.jackknife(data, user_statistic, kwargs=None, return_raw_distribution=False, pass_indices=False)

Compute first-order jackknife statistics of the data.

Parameters

dataarray_like

A 1-dimensional data array of size n_samples

user_statisticfunction

The statistic to be computed. This should take an array of data of size (n_samples, n_samples - 1) and return an array of size n_samples or tuple of arrays of size n_samples, representing the row-wise statistics of the input.

kwargsdictionary (optional)

A dictionary of keyword arguments to be passed to the user_statistic function.

return_raw_distributionboolean (optional)

if True, return the raw jackknife distribution. Be aware that this distribution is not reflective of the true distribution: it is simply an intermediate step in the jackknife calculation

pass_indicesboolean (optional)

if True, then the indices of the points rather than the points themselves are passed to user_statistic

Returns

mean, stdevfloats

The mean and standard deviation of the jackknifed distribution

raw_distributionndarray

Returned only if return_raw_distribution is True The array containing the raw distribution (length n_samples) Be aware that this distribution is not reflective of the true distribution: it is simply an intermediate step in the jackknife calculation

Notes

This implementation is a leave-one-out jackknife. Jackknife resampling is known to fail on rank-based statistics (e.g. median, quartiles, etc.) It works well on smooth statistics (e.g. mean, standard deviation, etc.)