1.包的安装
代码语言：javascript复制install.packages(c("devtools", "data.table", "wesanderson", "Seurat", "devtools", "AUCell", "GSEABase", "GSVA", "ggplot2","rsvd"))
devtools::install_github("YosefLab/VISION@v2.1.0") #Please note that the version would be v2.1.
devtools::install_github("wu-yc/scMetabolism")2.计算代谢激活分数
subcell_sce是已经注释过的Seurat对象。
代码语言：javascript复制load("sce.Rdata")
countexp.Seurat<-sc.metabolism.Seurat(obj = subcell_sce,
method = "AUCell",
imputation = F, ncores = 2,
metabolism.type = "KEGG")method有4种：VISION, AUCell, ssgsea和gsva。作者推荐VISION。
代码语言：javascript复制MedBioInfoCloud: rownames(countexp.Seurat@assays[["METABOLISM"]][["score"]])[1:6]
[1] "Glycolysis / Gluconeogenesis"
[2] "Citrate cycle (TCA cycle)"
[3] "Pentose phosphate pathway"
[4] "Pentose and glucuronate interconversions"
[5] "Fructose and mannose metabolism"
[6] "Galactose metabolism" 3.可视化
代码语言：javascript复制DimPlot.metabolism(obj = countexp.Seurat,
pathway = "Pentose phosphate pathway",
dimention.reduction.type = "umap",
dimention.reduction.run = F, size = 1)
DimPlot.metabolism(obj = countexp.Seurat,
pathway = "Citrate cycle (TCA cycle)",
dimention.reduction.type = "tsne",
dimention.reduction.run = F, size = 1)
input.pathway<-rownames(countexp.Seurat@assays[["METABOLISM"]][["score"]])[1:6]
DotPlot.metabolism(obj = countexp.Seurat, pathway = input.pathway,
phenotype = "seurat_clusters", norm = "y")
相关文章：
单细胞专题 | 1.单细胞测序（10×genomics技术）的原理 单细胞专题 | 2.如何开始单细胞RNASeq数据分析 单细胞专题 | 3.单细胞转录组的上游分析-从BCL到FASTQ 单细胞专题 | 4.单细胞转录组的上游分析-从SRA到FASTQ 单细胞专题 | 5.单细胞转录组的上游分析-从FASTQ到count矩阵 单细胞专题 | 6.单细胞下游分析——不同类型的数据读入 单细胞专题 | 7.单细胞下游分析——常规分析流程案例一 单细胞专题 | 8.单细胞类型注释之SingleR包详解 单细胞专题 | 9.如何人工注释单细胞类群？ 单细胞专题 | 10.细胞周期分析 单细胞专题 | 11.最新版cellphoneDB细胞通讯分析 单细胞专题 | 12.cellChat细胞通讯代码案例
参考：https://github.com/wu-yc/scMetabolism
