log likelihood = -4803.26