#引用R包 library(clusterProfiler) library(org.Hs.eg.db) #富集的基因列表 geneList = c("PCNA PSMD8 PS
#引用R包 library(R.utils) R.utils::setOption("clusterProfiler.download.method","auto
#自定义基因集 genName = c("PCNA PSMD8 PSMD7 SET SNRPA1 RAN SRSF2 G3BP1 STARD7 NPM1 BUB3 EIF3D XPO1
#以PDCD1为例 geneName = "PDCD1" geneNamedata = as.data.frame(data[geneName,]) colnames(g
#step1 提取出TCGA列名前12位相同的样本 #step2 TCGA列名第14位是1代表是正常组或癌旁,新增列tissue命名为"Paired_Normal",0代表肿