ods output ParameterEstimates=HR;
Proc PHReg Data=analysis;
   Model diag_time*diag(0) = X / risklimits=wald; *risklimits adds CI to ods ds*;
run;



*fill in HR if missing*;   *usually if there are interactions*;
Data HR;
   Set HR;
   if missing(HazardRatio) and not missing(Estimate) then do;
      HazardRatio = exp(estimate);
      HRLowerCL = exp(estimate - 1.96*stderr);
      HRUpperCL = exp(estimate + 1.96*stderr);
   end;
run;