文章预览
经过前面 单细胞免疫组库VDJ|从数据下载开始完成cellranger vdj分析(1) 中的cellranger count 和 cellranger vdj的分析后,得到了单细胞转录组 和 单细胞TCR的结果数据 。考虑到样本数太多, 选择既有RNA又有TCR样本的数据进行后续分析。重点解决以下三个小问题 1,只看TCR时(二),不同样本,不同分组情况下的clone差异以及变化; 2,TCR结合转录组数据(三+四),展示clone的分布可视化 以及 不同celltype中clone的差异以及变化; 3,一些可能的报错:(1)结合单细胞数据时;(2)官网链接中一些group相关的函数,按照推文或者使用??查看报错函数的具体说明 一 准备R包,数据集 安装 scRepertoire包,如果出现如下的Error,按照提示安装gsl后再安装 scRepertoire 。 BiocManager:: install ( "scRepertoire" ) suppressMessages( library (scRepertoire)) #Error: package or namespace load failed for ‘scRe
………………………………