Sunday, 22 October 2023

Lifelines Reference Group | Python Survival Analysis

Is there a way to calculate Hazard Ratios in Python's Lifeline package relative to a specified group?

What I have at the moment is:

from lifelines import CoxPHFitter

data=KM_DF[["Survival_Time","Event","30_40","40_50","50_60","60_70","70_80","Male"]]

cph=CoxPHFitter()
cph.fit(data,"Survival_Time",event_col="Event",show_progress=True)
cph.print_summary()

I would like to know the HR of the 40_50, 50_60, 60_70, and 70_80 group relative to the 30_40 group while also adjusting for sex. But at the moment, I get a HR of the group I want to be the reference ("30-40") as 3.02.

Is there a way to do this with lifelines?

For greater context/motivation, in R's workflow for Cox Proportional Hazard Ratio calculations, one can set variable "levels", and the subsequent Hazard Ratios are relative to the 1st level. However, I have been unable to sort out if this is possible in the popular "Lifelines" Python package.



from Lifelines Reference Group | Python Survival Analysis

No comments:

Post a Comment