log likelihood = -4877.65