scVelo运行2: 计算绘图

计算

1
2
3
4
5
6
7
8
9
scv.pp.moments(adata, n_pcs=30, n_neighbors=30)

scv.tl.recover_dynamics(adata, n_jobs=8)

scv.tl.velocity(adata, mode='dynamical')

scv.tl.velocity_graph(adata, n_jobs=8)

adata.write('hPB003.h5ad')

绘图

1
2
3
from matplotlib.pyplot import rc_context
with rc_context({'figure.figsize': (12, 12)}):
scv.pl.velocity_embedding_stream(adata, basis='umap', color=['cell_type_fig3'], save = "hPB003 velocity embedding stream.svg")

scVelo运行2: 计算绘图
https://occdn.limour.top/1817.html
Author
Limour
Posted on
May 8, 2022
Licensed under