# 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
通过conda安装纯净环境的ggsurvplot Previous
scVelo运行1: 合并数据 Next

#### Search

Keyword