cox.sat.dev function

Calculate the CoxPH saturated log-likelihood