title: "火山圖"
author: "劉大帥"
date: "2019/8/12"
output: html_document
火山圖事實上是一個以logFC為橫坐標(biāo)虑省，-log的p值為縱坐標(biāo)畫出的散點(diǎn)圖，以說明某個上調(diào)或下調(diào)的基因的狀態(tài)
可以試運(yùn)行該行代碼
plot(tT$logFC,-log10(tT$P.Value))
更加美觀的火山圖需要用到ggpubr包來實現(xiàn)
read.table('1.csv',sep=' ')
library(ggpubr)
attach(tT)
df=(tT)
df$v=-log10(P.Value)
ggscatter(df,x='logFC',y='v',size = 0.5)
df$g=ifelse(df$P.Value>0.01,'stable',
ifelse(df$logFC>1.5,'up',
ifelse(df$logFC<1.5,'down','stable')))##對數(shù)據(jù)進(jìn)行篩選分組
table(df$g)
df$NAME=rownames(df)
ggscatter(df,x='logFC',y='v',size = 0.5,color = 'g')
還有一些具體的參數(shù)可以到官網(wǎng)去查看
