We propose a novel algebraic framework for treating probability distributions
represented by their cumulants such as the mean and covariance matrix. As an
example, we consider the unsupervised learning problem of finding the subspace
on which several probability distributions agree. Instead of minimizing an
objective function involving the estimated cumulants, we show that by treating
the cumulants as elements of the polynomial ring we can directly solve the
problem, at a lower computational cost and with higher accuracy.