sufficient_statistics function

Compute sufficient statistics