log likelihood = -36991.25