Change distance metric in PCA

#1

Hi,
I need to do a PCA using the Jaccard distance metric instead of euclidean distance. Is it possible to change the distance metric of pca() or hwe_normalized_pca() in any way?

Thanks!

#2

You want to do something like the algorithm proposed in this paper: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5860507/?

This unfortunately requires a totally different algorithm, which isn’t really PCA, so it’s not something we can tweak easily. Perhaps @jbloom could comment further.

#3

Yes, exactly that paper

#4

@jbloom what are your thoughts on this?

#5

@patrick-schultz and I were actually investigating Jaacaard similarity (and other generalizations of min hashing) around the same time as that paper, and he worked out some cool theoretical extensions that were challenging to implement in Hail using the infrastructure at the time. Patrick, any thoughts on where implementing this at the level of Python might fit into the current development arc?

#6

That would be awesome!