Starlitnightly / omicverse

A python library for multi omics included bulk, single cell and spatial RNA-seq analysis.
https://starlitnightly.github.io/omicverse/
GNU General Public License v3.0
277 stars 32 forks source link

How to perform scVelo analysis on Seurat subtypes #101

Open asdfasf4safasg opened 3 days ago

asdfasf4safasg commented 3 days ago

When I was performing scVelo analysis, I first divided the single-cell data into main clusters, then extracted a group from the main clusters and further subdivided it into subtypes. However, when I directly used the UMAP analysis on the subset subtypes, I found that the analysis results were different from the pattern of the subset. Therefore, I had to map the subtypes back to the main clusters for analysis. I would like to ask how I can directly analyze the results of the subset. cur_celltypes = ['Astrocytes','Interneurons'] "adata_subset = adata[adata.obs['celltype'].isin(cur_celltypes)]"

Starlitnightly commented 3 days ago

Sry, i didn't understand what's mean of I found that the analysis results were different from the pattern of the subset.

Perhaps you could use a diagram to represent your question in analyse.

Regards, Zehua

asdfasf4safasg commented 3 days ago
image image image

After subsetting the colored figure from the main figure and reanalyzing it, how can I, in subsequent scVelo analyses, focus only on the subtype part but ensure the analysis is based on the subtype itself and not just on the main figure’s pattern?

Starlitnightly commented 3 days ago

image

All the analyses you are doing on this diagram are subtype related, not related to all the cells, you can re-run the clustering algorithms like leiden on subtype if you are not satisfied with raw clusters.

Zehua

asdfasf4safasg commented 3 days ago

不好意思,刚才微信只是想要留言,没有要你及时回复的意思。我能不能用中文回复,我可能英文没有描述清楚,就是scvelo的分析我想要的pattern是从main celltype中提取的subtype的pattern,而不是main celltype在scvelo分析的时候提取子集分析然后用的main celltype的pattern

从 Windows 版邮件https://go.microsoft.com/fwlink/?LinkId=550986发送


发件人: Zehua Zeng @.> 发送时间: Friday, July 5, 2024 11:14:08 PM 收件人: Starlitnightly/omicverse @.> 抄送: Lili ZHOU 周丽丽 @.>; Author @.> 主题: Re: [Starlitnightly/omicverse] How to perform scVelo analysis on Seurat subtypes (Issue #101)

你通常不会收到来自 @.*** 的电子邮件。了解这一点为什么很重要https://aka.ms/LearnAboutSenderIdentification

image.png (view on web)https://github.com/Starlitnightly/omicverse/assets/46667721/6512193e-ad15-47ab-8085-d54f65bb4f4d

All the analyses you are doing on this diagram are subtype related, not related to all the cells, you can re-run the clustering algorithms like leiden on subtype if you are not satisfied with raw clusters.

Zehua

― Reply to this email directly, view it on GitHubhttps://github.com/Starlitnightly/omicverse/issues/101#issuecomment-2211049179, or unsubscribehttps://github.com/notifications/unsubscribe-auth/BAC2M4TT7J4I7K64PZ2QHP3ZK2Z35AVCNFSM6AAAAABKNGIHEWVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDEMJRGA2DSMJXHE. You are receiving this because you authored the thread.Message ID: @.***>

Starlitnightly commented 3 days ago

没太看懂你的问题,你要计算subtype在subtype上单独运行你的分析就好,例如在subtype上运行velocity。如果英文都表述不清楚那中文更难表述清楚。例如,如果你是想保留main celltype的umap结果,那你就不要重新进行邻域图计算(ov.pp.neighbors),保留原来的邻域图的distance和connectivity。你告诉我这些单词我换成中文你能明白我在说什么吗?

asdfasf4safasg commented 3 days ago

是这样的,当我只用subtype的数据分析的时候,得到的stream的结果发生了改变,我用main celltype的时候不会出现这个问题

从 Windows 版邮件https://go.microsoft.com/fwlink/?LinkId=550986发送


发件人: Zehua Zeng @.> 发送时间: Saturday, July 6, 2024 12:06:44 AM 收件人: Starlitnightly/omicverse @.> 抄送: Lili ZHOU 周丽丽 @.>; Author @.> 主题: Re: [Starlitnightly/omicverse] How to perform scVelo analysis on Seurat subtypes (Issue #101)

你通常不会收到来自 @.*** 的电子邮件。了解这一点为什么很重要https://aka.ms/LearnAboutSenderIdentification

没太看懂你的问题,你要计算subtype在subtype上单独运行你的分析就好,例如在subtype上运行velocity。如果英文都表述不清楚那中文更难表述清楚。例如,如果你是想保留main celltype的umap结果,那你就不要重新进行邻域图计算(ov.pp.neighbors),保留原来的邻域图的distance和connectivity。你告诉我这些单词我换成中文你能明白我在说什么吗?

― Reply to this email directly, view it on GitHubhttps://github.com/Starlitnightly/omicverse/issues/101#issuecomment-2211119047, or unsubscribehttps://github.com/notifications/unsubscribe-auth/BAC2M4WV5PXHGFCNHM7AVMDZK3ABJAVCNFSM6AAAAABKNGIHEWVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDEMJRGEYTSMBUG4. You are receiving this because you authored the thread.Message ID: @.***>

Starlitnightly commented 3 days ago

What's mean of stream? What does `stream改变' mean?

asdfasf4safasg commented 2 days ago

image 我是想指这个速率的箭头结果的umap在提取子集分析的时候会发生改变,但是用main celltype的umap不会出现这个问题

Starlitnightly commented 1 day ago

image 我是想指这个速率的箭头结果的umap在提取子集分析的时候会发生改变,但是用main celltype的umap不会出现这个问题

我不觉得有什么区别?

image