第一步 读入数据1234567891011121314151617# 配置数据路径root_path = "~/zlliu/R_data/TCGA" # 配置结果保存路径output_path = root_pathif (!file.exists(output_path)){dir.create(output_path)} # 设置工作目录，输出文件将保存

第一步 转换数据1234567891011121314151617181920212223242526272829303132library('biomaRt')mart <- useDataset("hsapiens_gene_ensembl", useMart("ensembl"))genes <- rownames(ge

第一步 常规配置12345678910111213141516171819library(Matrix)library(Seurat)library(plyr)library(dplyr)library(patchwork)library(purrr)rm(list = ls())# 配置数据路径root_path = "~/zlliu/R_data/TCGA" # 配置结果保

第一步 导入程辑包1234library(cgdsr)library(survival)library(survminer)library(ggpubr) 12345678910# 配置数据路径root_path = "~/zlliu/R_data/TCGA" # 配置结果保存路径output_path = root_pathif (!file.exists(output_p

第一步 读取数据 TCGA 12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970library(TCGAbiolinks)library(plyr)library(limma)librar

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697989910010110210

第一步 火山图12345678910111213141516171819202122232425Idents(sc_Neuron) <- sc_Neuron[['hmca_class']]all_markers <- FindAllMarkers(sc_Neuron, min.pct = 0.25, logfc.threshold = 0.25)significan

第一步 火山图123456789101112131415161718192021library(ggplot2)library(ggrepel)f_DEG_Volcano <- function(lc_logFC, lc_p, lc_gene, Threshold_logFC = 1, Threshold_p = 0.05, lc_rep=1:10){ col_vector

12345678910111213141516171819202122232425262728293031323334353637383940414243f_prepare4CSOmap <- function(lc_scRNA, lc_csomap_data_dir, lc_className){ lc_csomap_data_dir <- system(paste(

第一步 导出数据123456789101112f_prepare4cellphoneDB <- function(lc_scRNA, lc_dir, lc_className){ if (!file.exists(lc_dir)){dir.create(lc_dir)} # 生成 count.txt write.table(as.matrix