The UMAP generated by `reduceDimension` is rather "weak"

Hi my dear friends,

I noticed the umap generated by dyn.tl.reduceDimension is not as good as that from other tools (seurat, scanpy, etc.)

Here is a umap generated by dyn.tl.reduceDimension (preprocessed by recipe_monocle). We can see the clusters formed a condensed round mass. These cells are on a developmental trajectory. It can form a line structure according to our previous experiences.

However, the umap generated by scanpy is better structured - the clusters correctly formed along some axises. (the same loom file as used above, the same preprocessed adata.X but fed to scanpy package for downstream umap analysis ) In this case, I used size-factor normalized and log1p transformed adata.X generated by recipe_monocle (I didn’t use scanpy’s preprocessing steps, I know there are normalization differences), used hvg genes to calculated pca, neighbor graphs, and umap in scanpy. Since the input are the same (recipe_monocle normalized data), why do we got quite different umap? (dynamo’s umaps seem less informative)

Moreover, it seems that even scanpy’s pca plot has a better structure preservation ability. we all know umap is non-linear after all.

These situations are not limited to this single dataset. All dataset I tried (up to now) hints that the UMAP produced by dynamo is rather “weak” (not as good as a non-linear method’s expectation)

Looking forward to your suggestions and guidance.

p.s. I raised a lot of questions/problems recently. And hope I didn’t bother your new-manuscript preparations too much.

[chansigit]

For your information, I deleted the umap embedding and recomputed umap with dyn.tl.reduceDimension(adata2,kwargs={'min_dist':0.01},enforce=True) only to find it yields exactly the same umap

The default min_dist=0.5 in dynamo https://github.com/aristoteleo/dynamo-release/blob/1d1f5c521d0b8763ab866dcc3c3d96be0ba3f8f3/dynamo/tools/utils_reduceDimension.py#L203,

The default min_dist of scanpy is also 0.5,see https://scanpy.readthedocs.io/en/latest/generated/scanpy.tl.umap.html

1. it seems that min_dist is not the cause.
2. it seems that passing a min_dist in kwargs has no impact at all.

[github-actions[bot]]

This issue is stale because it has been open 90 days with no activity. Remove stale label or comment or this will be closed in 14 days