compute_sufficient_statistics function

Compute sufficient statistics