log likelihood = -41100.25