log likelihood = -41075.62