基于voerlap的相似性度量-单细胞克隆亚型分析

基于voerlap的相似性度量-单细胞克隆亚型分析

一、背景

基于overlap的相似性度量,有多种多样的方法,其应用也是多种多样。这里,我们主要讲解在单细胞免疫组库中的应用。

基于单细胞的免疫组库的测序,我们可以获得每一个细胞类群中的不用克隆亚型,基于这些克隆亚型的overlap,我们可以推断不同细胞亚群的生命体的状态变化,或者是进化轨迹分析。

二、基于overlap的相似性度量

2.1 Jaccard 系数

Jaccard index 衡量有限样本集合之间的相似度,定义为交集的大小除以样本集并集的大小。
J a c c a r d   i n d e x = A ∩ B A ∪ B Jaccard \ index = \frac{A \cap B}{A \cup B} Jaccard index=ABAB

2.2 Cosine similarity

Cosine similarity是一种余弦相似性指标,适用于两个非零向量(只要你的数据能够转化)

基于两个向量: a ⃗ = ( x 1 , y 1 ) \vec a = (x_1, y_1) a =(x1,y1) b ⃗ = ( x 2 , y 2 ) \vec b = (x_2, y_2) b =(x2,y2)
C o s i n e   s i m i l a r i t y = a ⋅ b ∣ ∣ a ∣ ∣ × ∣ ∣ b ∣ ∣ = x 1 x 2 + y 1 y 2 x 1 2 + y 1 2 × x 2 2 + y 2 2 Cosine \ similarity = \frac{a \cdot b}{||a|| \times ||b||} = \frac{x_{1}x_{2}+y_{1}y_{2}}{\sqrt{x_1^2+y_1^2} \times \sqrt{x_2^2 + y_2^2}} Cosine similarity=∣∣a∣∣×∣∣b∣∣ab=x12+y12 ×x22+y22 x1x2+y1y2

2.3 Morisita-Horn similarity

M o r i s i t a   H o r n   s i m i l a r i t y = ∑ i = 1 n 2 × X i × Y i ∑ i = 1 n X i 2 + ∑ i = 1 n Y i 2 Morisita \ Horn \ similarity = \frac{\sum\limits_{i=1}^{n}{2 \times X_i \times Y_i}}{\sum\limits_{i=1}^{n}{X_i^2} + \sum\limits_{i=1}^{n}{Y_i^2}} Morisita Horn similarity=i=1nXi2+i=1nYi2i=1n2×Xi×Yi

其中,

  • X i X_i Xi表示在 X X X 集合中,第 i i i 类别的个数(或相对个数)
  • Y i Y_i Yi表示在 Y Y Y 集合中,第 i i i 类别的个数(或相对个数)

比如,在两种组织的单细胞亚型的相似性度量:
M H   i n d e x = 2 × ∑ i = 1 n T i × B i ∑ i = 1 n ( T i 2 + B i 2 ) T i = t i T N B i = b i B N MH \ index = \frac{2 \times \sum\limits_{i=1}^{n}{T_i \times B_i}}{\sum\limits_{i=1}^{n}{(T_i^2 + B_i^2)}} \\ T_i = \frac{t_i}{T_N} \\ B_i = \frac{b_i}{B_N} MH index=i=1n(Ti2+Bi2)2×i=1nTi×BiTi=TNtiBi=BNbi
其中,

  • b i b_i bi t i t_i ti 分别表示在血液和肿瘤组织内,第 i i i 个克隆亚型的细胞的数量
  • B N B_N BN T N T_N TN 分别表示在血液和肿瘤组织内,所有克隆亚型的细胞数量

三、个人理解

  • Jaccard index通常是指标签unique之后,进行重叠相似性分析
  • Cosine index既可以是unique之后的,也可以是每个样本类型中,重叠的数量,看你如何构建两个向量了
  • Morisita-Horn index通常是基于重叠的类型在各个样本组织中的数量或者比值进行计算的
  • 可以通过R包 immunarch 实现更多的分析
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值