整合多个DEGs 第一步 读入数据1234567891011121314151617# 配置数据路径root_path = "~/zlliu/R_data/TCGA" # 配置结果保存路径output_path = root_pathif (!file.exists(output_path)){dir.create(output_path)} # 设置工作目录,输出文件将保存 2021-10-11 单细胞下游分析
SigEMD (一) 简单教程 第一步 转换数据1234567891011121314151617181920212223242526272829303132library('biomaRt')mart <- useDataset("hsapiens_gene_ensembl", useMart("ensembl"))genes <- rownames(ge 2021-10-11 生物信息学
Seurat (六) 整合数据(二) 第一步 常规配置12345678910111213141516171819library(Matrix)library(Seurat)library(plyr)library(dplyr)library(patchwork)library(purrr)rm(list = ls())# 配置数据路径root_path = "~/zlliu/R_data/TCGA" # 配置结果保 2021-10-09 生物信息学
TCGA (四) 生存分析 第一步 导入程辑包1234library(cgdsr)library(survival)library(survminer)library(ggpubr) 12345678910# 配置数据路径root_path = "~/zlliu/R_data/TCGA" # 配置结果保存路径output_path = root_pathif (!file.exists(output_p 2021-10-05 组织测序分析
TCGA (三) 两组间DEGs分析 第一步 读取数据 TCGA 12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970library(TCGAbiolinks)library(plyr)library(limma)librar 2021-10-05 组织测序分析
Seurat (五) 简单总结 12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697989910010110210 2021-10-04 生物信息学
Seurat (四) 差异基因分析 第一步 火山图12345678910111213141516171819202122232425Idents(sc_Neuron) <- sc_Neuron[['hmca_class']]all_markers <- FindAllMarkers(sc_Neuron, min.pct = 0.25, logfc.threshold = 0.25)significan 2021-10-03 生物信息学
GEO (一) 绘图总结 第一步 火山图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 2021-10-03 组织测序分析
CSOmap (二) 简化流程 12345678910111213141516171819202122232425262728293031323334353637383940414243f_prepare4CSOmap <- function(lc_scRNA, lc_csomap_data_dir, lc_className){ lc_csomap_data_dir <- system(paste( 2021-10-03 单细胞下游分析
cellphoneDB (一) 使用流程 第一步 导出数据123456789101112f_prepare4cellphoneDB <- function(lc_scRNA, lc_dir, lc_className){ if (!file.exists(lc_dir)){dir.create(lc_dir)} # 生成 count.txt write.table(as.matrix 2021-10-02 单细胞下游分析