added a version that updates sufficient stats to get rid of T dependent matrix multiplies