How to assess the similarity between programs across samples

[dylkot]

How to assess the similarity between programs across samples

[michelle-curtis]

Hi, thanks for your comment. It is effective to use either the Z-score or TPM spectra matrix to compare correlation across programs. However, if using the TPM matrix, you will want to variance-normalize it to allow genes expressed at different scales to contribute equally to the program correlations - gene variances have already been output in the tpm_stats file. We also suggest using either the union or intersection of the HVGs across samples for these comparisons, especially when using the var-normed TPM-spectra. The TPM spectra matrix tends to have more baseline correlation across the non-differentially expressed genes (whereas these genes will tend to have zero values in the Z-score matrix).

[dylkot]

Yes, like @michelle-curtis said, either can work. I think using the Z-score output by default seems to work the best. And also as Michelle said, good to subset to high variance genes