Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

velocity_pseudotime: from scv.tl.latent_time vs scv.tl.velocity_pseudotime #3588

Closed
savagyan00 opened this issue Dec 10, 2023 · 2 comments
Closed

Comments

@savagyan00
Copy link

I am running the basic pipeline (scvelo 0.3.1), starting from an AnnData object:

  • Preprocessing
scv.pp.filter_and_normalize(
    adata_maca_10x, min_shared_counts=20, n_top_genes=2000, subset_highly_variable=False
)

sc.tl.pca(adata_maca_10x)
sc.pp.neighbors(adata_maca_10x, n_pcs=30, n_neighbors=30, random_state=0)
scv.pp.moments(adata_maca_10x, n_pcs=None, n_neighbors=None)
  • run velocity
scv.tl.recover_dynamics(adata_maca_smartseq, n_jobs=8)
scv.tl.velocity(adata_maca_10x, mode="dynamical")

after this step I take one two ways

  1. run latent_time : this generates a column "velocity pseudotime" in obs
scv.tl.latent_time(adata_maca_10x)

(when I run velocity_pseudotime after this step without velocity_graph, also different values are generated)

  1. run velocity pseudotime: this also generates a column "velocity pseudotime" in obs
scv.tl.velocity_graph(adata_maca_10x)
scv.tl.velocity_pseudotime(adata_maca_10x)

But the values from these two are completely different, which I am unsure how to interpret. Can you please explain or point to documentation where I can find descriptions for these two value outputs?

Thanks!

@mwaskom
Copy link
Owner

mwaskom commented Dec 10, 2023

I think you filed this in the wrong repository?

@savagyan00
Copy link
Author

So sorry, I confused the repository.

@mwaskom mwaskom closed this as not planned Won't fix, can't repro, duplicate, stale Dec 10, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants