log likelihood = -37130.71