icfit function

Fit a proportional hazards model with baseline hazard modeled by P-splines