log_lik function

Evaluate the negative log-likelihood of an intercept-only Gaussian plmm model