loglikssaturated function

Calculates the loglikelihoods of the saturated model