1. library(GSVA)
    2. setwd("F:\\pancancer lncRNA immune\\全体集合结果\\单样本富集分析")
    3. data <- read.table("data_mrna.txt",stringsAsFactors = F)
    4. duiying <- read.table("ENSG-Gene symbol.txt",stringsAsFactors = F,sep = "\t")
    5. weizhi <- match(rownames(data),duiying[,2])
    6. data <- as.matrix(data)
    7. rownames(data) <- duiying[weizhi,1]
    8. gs <- read.csv("immune.csv",stringsAsFactors = F,row.names = 1)
    9. gs <- t(gs)
    10. gs <- as.data.frame(gs)
    11. gs <- as.list(gs)
    12. gs <- lapply(gs, function(x) x[which(x != "")])
    13. jieguo <- gsva(data,gs, method = "ssgsea", ssgsea.norm = TRUE, verbose = TRUE)
    14. #write.csv(jieguo,"immune_ssGSEA.csv",quote = F)