Track down and fix two hidden bugs in mvnormal and multinomial
Use sum(x) and sum((x - mu)^2) as Normal's sufficient statistics