ggplot火山图的绘制（volcano plot）
library(ggplot2)
volcano<-subset(data,select = c(Pvalue,Foldchange))
threshold<-as.factor((log2(volcano$Foldchange)>1.5|log2(volcano$Foldchange)<(-1.5))&volcano$Pvalue<0.05)
r03=ggplot(volcano,aes(log2(Foldchange),-log2(Pvalue),colour=threshold))+geom_point()
r04=r03+labs(title="Volcanoplot")+theme(plot.title = element_text(hjust = 0.5))+xlim(-10,10)
r05=r04+geom_vline(xintercept=c(-1.5,1.5),linetype="dotted",size=1)+geom_hline(yintercept=-log2(0.05),col="blue")
文章图片
火山图
| 基因名 | Pvalue | Foldchange |
|---|---|---|
| ENSG00000000003 | 0.868134121987662 | 1.2 |
| ENSG00000000457 | 0.256006872721537 | 1.47826086956522 |
| ENSG00000000971 | 0.932736306872995 | 0.961486486486487 |
| ENSG00000001084 | 0.512134516499327 | 0.555555555555556 |
| ENSG00000001167 | 0.175213431516265 | 0.291666666666667 |
| ENSG00000001617 | 0.331530822974447 | 1.58441558441558 |
| ENSG00000001630 | 0.407563817673791 | 0.588235294117647 |
【ggplot火山图的绘制（volcano plot）】
推荐阅读
- 笔记|R语言画图 | 如何看已知基因list的细胞类型特异性表达（）
- 拓端tecdat|拓端tecdat|Python多项式Logistic逻辑回归进行多类别分类和交叉验证准确度箱线图可视化
- 线性代数在前端中的应用（二）（实现鼠标拖拽旋转元素、Canvas图形）
- C++版图书管理系统
- OpenCV实现抠图工具
- C++项目开发实现图书管理系统
- 超好用的两款作图工具，用起来~~~
- mysql|文末送书 | WAF 那些事
- 转载--vlan的基本概念
- JAVA的自学之路|【Java课程设计】【图形界面程序】【计算器】