Skip to content

Commit e6a28ac

Browse files
authored
Update velocity.md
1 parent 66d7080 commit e6a28ac

File tree

1 file changed

+2
-22
lines changed

1 file changed

+2
-22
lines changed

scRNAseq/scRNAseq_analysis_tutorial/lessons/velocity.md

Lines changed: 2 additions & 22 deletions
Original file line numberDiff line numberDiff line change
@@ -88,26 +88,6 @@
8888
import scvelo as scv
8989
adata = scv.read("mouseBM.h5ad")
9090
adata
91-
92-
# This may not be the best function to use - normalization procedure not as robust as scTransform
93-
scv.pp.filter_and_normalize(adata, min_shared_counts=20, n_top_genes=3000)
94-
95-
scv.pp.moments(adata, n_pcs=30, n_neighbors=30)
96-
97-
scv.tl.velocity(adata)
98-
99-
scv.tl.velocity_graph(adata)
100-
101-
scv.pl.velocity_embedding_stream(adata, basis="umap", color="seurat_clusters")
102-
103-
scv.pl.velocity_embedding(adata, basis="umap", color="seurat_clusters", arrow_length=3, arrow_size=2, dpi=120)
104-
105-
scv.tl.recover_dynamics(adata)
106-
107-
scv.tl.latent_time(adata)
108-
109-
scv.pl.scatter(adata, color="latent_time", color_map="gnuplot")
110-
111-
top_genes = adata.var["fit_likelihood"].sort_values(ascending=False).index[:300]
112-
scv.pl.heatmap(adata, var_names=top_genes, sortby="latent_time", col_color="seurat_clusters", n_convolve=100)
11391
```
92+
93+
The scVelo code used for scVelo analysis can be found [here](https://www.dropbox.com/s/jxlworc6td3mfcy/velocity_jupyter_notebook.pdf?dl=1).

0 commit comments

Comments
 (0)