Added joint log likelihood calculation