PBMC mass cytometry data from 10 subjects of cohort-1 or 8 subjects of cohort-2 was combined and batch normalized using harmony respectively. First, mass cytometry data from each cohort all subjects was combined into a single dataset. Second, harmony batch correction was performed using one of the samples. Third, mass cytometry data were lognormalized in the Seurat’s NormalizeData function across the aggregated dataset.