You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Hi, I have a anndata converted by seurat data which adata.X was scaled and have a umap. I found adata.X in tutorial was only normalised, not scaled. I am not sure whether the expression matrix will be used except dimensionality reduction. I want to know if I can run it as below:
adata=scv.read('adata.h5ad') #adata.X is normalised and scaled and have a umap
loom=scv.read('my.loom', cache=True) #velocyto: run10x output
rep2=scv.utils.merge(adata,loom)
scv.pp.filter_and_normalize(rep2, min_shared_counts=20, n_top_genes=3000) # normalised expression matrix was skipped. Only for spliced and unspliced matrix.
scv.pp.moments(rep2, n_pcs=30, n_neighbors=30)
scv.tl.recover_dynamics(rep2)
scv.tl.velocity(rep2, mode='dynamical')
scv.tl.velocity_graph(rep2)
scv.pl.velocity_embedding_stream(rep2, basis='umap',color='celltype',save='umap-velocity_embedding.png')
reacted with thumbs up emoji reacted with thumbs down emoji reacted with laugh emoji reacted with hooray emoji reacted with confused emoji reacted with heart emoji reacted with rocket emoji reacted with eyes emoji
-
Hi, I have a anndata converted by seurat data which adata.X was scaled and have a umap. I found adata.X in tutorial was only normalised, not scaled. I am not sure whether the expression matrix will be used except dimensionality reduction. I want to know if I can run it as below:
Thanks!
Beta Was this translation helpful? Give feedback.
All reactions