其实不同的打分算法偏好不太一样 应该视数据而定
1.制作基因集合:是一个list,并且list下是基因集的命名
## 制作
geneset=read.table('./GTs.txt',header = T,check.names = F,sep = '\t')
geneset=list(geneset$GT)
names(geneset)='GTs'
gc()
2.R包
library(Seurat)
library(SeuratData)
library(UCell)
library(irGSEA)
library(AUCell) ## 版本需要大于1.14
3.打分 以单细胞数据为例
## 1. AUCell
cells_rankings <- AUCell_buildRankings(scRNA@assays$RNA@data,
nCores=1,
plotStats=TRUE)
gc()
cells_AUC <- AUCel