Joint statistical modeling to adjust for confounding of demographic factors (age, sex, and race) when testing for association of UMAP Components 1 and 2 with the NIH Ordinal Severity Scale was performed using ordinal logistic regression provided by the polr function of the R package MASS.