log likelihood = -36941.38