log likelihood = -36995.11