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 #1168

Open
savagyan00 opened this issue Dec 10, 2023 · 1 comment

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

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)

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!

@sbombin
Copy link

sbombin commented Apr 12, 2024

I have the same issue. Did you find the answer?

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