数据库代码公开：AUCell评估单细胞打分

#自定义基因集
genName = c("PCNA PSMD8 PSMD7 SET SNRPA1 RAN SRSF2 G3BP1 STARD7 NPM1 BUB3 EIF3D XPO1 FBL EIF4A1 CANX NAP1L1 CBX3 CCT3 C1QBP U2AF1 UBE2L3 SSBP1 SRSF1 TCP1 MCM2 EIF3B PSMD14 SNRPA PWP1 APEX1 TXNL4A HNRNPR PSMB2 HPRT1 MCM6 NME1 SNRPD1 EEF1B2 HSPD1 CAD RPL18 PGK1 DDX18 RPS2 LDHA RUVBL2 RNPS1 EIF2S1 RANBP1 MCM5 IARS1 UBE2E1 MYC")
genName = unlist(strsplit(genName," "))
#自定义数据集
CancerName <- "BRCA_GSE148673"
#引用R包
library(Seurat)
library(ggplot2)
library(patchwork)
library(dplyr)
library(Matrix)
library(AUCell)
#读取单细胞转录组
scRNA = readRDS(paste0(CancerName,".RDS"))
#自定义基因集
genName = intersect(genName,rownames(scRNA))
geneSets = list(Custom=genName)
cells_rankings = AUCell_buildRankings(GetAssayData(object = scRNA, assay = "RNA") %>%  as.matrix(), nCores=2)  #基因排序，使用10个核，加速计算
cells_AUC = AUCell_calcAUC(geneSets,cells_rankings,nCores = 1,aucMaxRank = nrow(cells_rankings)*0.05)
AUC_Exp = as.numeric(getAUC(cells_AUC)["Custom",,drop=F])
scRNA$AUC = AUC_Exp