I need to test, which cox model fits the best based on partial AIC. Below you see the code, that I am using to test the combination of 5 parameters. How can I adjust my code, to test all possible variables combinations?
parameters = ['v1', 'v3', 'v4','v5','CRD','followup']
p = []
for parameter in parameters:
p.append(parameter)
print(p)
cph = CoxPHFitter()
cph.fit(df[p], duration_col="followup", event_col="CRD", show_progress=True, step_size = .5)
cph.print_summary(decimals = 5)