Multivariable analysis

output HR (95% CI) and P value