To assess how parameters chosen in the clustering algorithm affected final cluster membership, we re-ran the clustering workflow with a range of values for k-NN (5–50) and a range of values for the number of LSI components used (10–50). For each clustering, we calculated the adjusted Rand index between the clustering results obtained and the cell type labels derived from the gene expression assay using the mclust R package63.