【scMetabolism】一站式代谢通路基因集打分

阅读 117

2022-04-26

install.packages(c("devtools", "data.table", "wesanderson", "Seurat", "devtools", "AUCell", "GSEABase", "GSVA", "ggplot2","rsvd"))

devtools::install_github("YosefLab/VISION")
devtools::install_github("wu-yc/scMetabolism")
load(file = "pbmc_demo.rda")

library(scMetabolism)
library(ggplot2)
library(rsvd)

quantify single cell metabolism in seurat

countexp.Seurat<-sc.metabolism.Seurat(obj = countexp.Seurat, method = "VISION", imputation = F, ncores = 2, metabolism.type = "KEGG")

metabolism.matrix <- countexp.Seurat@assays$METABOLISM$score

精彩评论(0)

0 0 举报