All covariates except tumor size were modeled as categorical variables using dummy variables, whereas tumor size was modeled as continuous variable using a 3-knot restricted cubic spline[12] to obtain flexible fits and to allow for its prognostic effect as a variable that was not equivalent in every part of the range.