title: "火山圖"

author: "劉大帥"

date: "2019/8/12"

output: html_document

火山圖事實(shí)上是一個(gè)以logFC為橫坐標(biāo)，-log的p值為縱坐標(biāo)畫(huà)出的散點(diǎn)圖，以說(shuō)明某個(gè)上調(diào)或下調(diào)的基因的狀態(tài)

可以試運(yùn)行該行代碼

 plot(tT$logFC,-log10(tT$P.Value))

更加美觀的火山圖需要用到ggpubr包來(lái)實(shí)現(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')))##對(duì)數(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)去查看