脚本更新----空间转录组系列分析(cell2location + 细胞聚类 + 空间通讯)

作者,Evil Genius

我们今天更新一个内容,那就是visium的通讯分析。

最近也接触了很多分析空转的人,从我的角度来看,我觉得很多人还没有达到可以分析空转的水准,不过个人感觉啊,也可能是错觉。

我还是建议大家先学一学,把基础做好,有没有钱,课题组支不支持,我们无法决定,但是呢,培训班也没有说完全不透风,盗版很多的,大家学一学再去分析会大有裨益。

今天我们需要分析的目标



首先我们看一下需要分析的内容

1、空间转录组的数据比对
2、空间转录组解卷积分析
3、空间转录组细胞聚类
4、空间转录组共定位分析和差异分析
5、空间转录组通讯分析
6、空间转录组细胞类型与通讯的关联分析。

第一部分我们就不看了,直接上第二部分,空间转录组解卷积分析

之前讲过,单细胞空间联合,cell2location是高分标配,一般来讲,我们cell2location、RCTD二选一。
分析前需要处理好两个问题

1、单细胞数据定义好,允许的范围内越细越好
2、样本配对,优先一一配对样本,次选同类型样本。

cell2location分享过很多了,脚本大家自己找一份,解卷积分析即可,但是有几个问题大家要考虑。
问题1、估计细胞类型特征的时候,基因该如何选择?
问题2、关键性的参数怎么写,detection_alpha和N_cells_per_location
# create and train the model
mod = cell2location.models.Cell2location(
    adata_vis, cell_state_df=inf_aver,
    # the expected average cell abundance: tissue-dependent
    # hyper-prior which can be estimated from paired histology:
    N_cells_per_location=5,
    # hyperparameter controlling normalisation of
    # within-experiment variation in RNA detection:
    detection_alpha=200
)
第三部分:空间转录组细胞聚类
这部分也分享了很多了,课程上也都讲过,最近也分享了RCTD的解卷积分析 + 层次聚类。
那么通常的做法就是cell2location + leiden分析为主。
直接采用cell2location解卷积的结果 + leiden聚类分析(其中需要注意resolution的参数)
还有 88% 的精彩内容
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
禁止转载,如需转载请通过简信或评论联系作者。
支付 ¥200.00 继续阅读

推荐阅读更多精彩内容