statstream.exact.streaming_mean_and_cov

statstream.exact.streaming_mean_and_cov(X, steps=None)

Mean and covariance matrix of a streaming dataset.

Computes the mean and the covariance matrix of a dataset from a stream of batches of samples. The data has to be provided by an iterator yielding batches of samples. Either a number of steps can be specified, or the iterator is assumed to be emptied in a finite number of steps. In the first case only the given number of batches is extracted from the iterator and used for the mean and covariance calculation, even if the iterator could yield more data.

Samples are given along the first axis. The mean has the same shape as the remaining axes, e.g. batches of shape [batch_size, d1, ..., dN] will produce a mean of shape [d1, ..., dN]. Covariances are arranged in the squared shape [d1, ..., dN, d1, ..., dN].

This function consumes an iterator, thus finite iterators will be empty after a call to this function, unless steps is set to a smaller number than batches in the iterator.

Parameters:
Xiterable

An iterator yielding batches of samples.

stepsint, optional

The number of batches to use from the iterator (all available batches are used if set to None). The defaul is None.

Returns:
array

The mean of the seen data samples.

array

The covariance matrix of the seen data samples.

Warning

Use this function only on data sets of reasonably small dimensions.

Full covariances matrices are costly to compute and require large amounts of memory. The shape of the covariance matrix is squared the size of each individual sample in the data set.

If your data is high dimensional and you do not need the exact covariance matrix, consider using streaming_mean_and_low_rank_cov or streaming_low_rank_cov from statstream.approximate instead.

See also

streaming_mean

get only the mean in a single pass.

streaming_cov

get only the covariance matrix in a single pass.

Notes

The streamed covariance calculation is a generalization of the streamed variance calculation as described in [1].

References

[1]

Tony F. Chan & Gene H. Golub & Randall J. LeVeque, “Updating formulae and a pairwise algorithm for computing sample variances”, 1979.