火山图在转录组分析中运用较多,常用于展示显著表达的差异基因,以下是绘制基础火山图的一个小示例。
1.调取ggplot2包,读取数据文件;
#####火山图
library(ggplot2)
#读取数据
dataset <- read.table('volcano_data.txt',header = TRUE)
2.设置阈值;
# 设置pvalue和logFC的阈值
cut_off_pvalue = 0.01
cut_off_logFC = 1
dataset$change = ifelse(dataset