log likelihood = -41113.94