Survival analysis

For continuous univariate and multivariate survival analyses, infiltration scores and the appropriate clinical covariates were fit to a Cox proportional hazards model using the “coxph” function from the R “survival” package. Survival distributions for different cell types were visualized using Kaplan-Meier curves created by the “survfit” function from the R “survival” package. The median infiltration score was used to stratify patients into “high” and “low” infiltration groups when performing univariate two-class comparisons, while an infiltration score of 0 was used to stratify the two groups when performing multivariate analyses. Differences between the survival distributions in each Kaplan-Meier plot were calculated using a log-rank test through the “survdiff” function from the R “survival” package.