McCullagh_compute_log_l function

Computes the log(likelihood) for the general nonlinear model.