title: "火山图"

author: "刘大帅"

date: "2019/8/12"

output: html_document

火山图事实上是一个以logFC为横坐标，-log的p值为纵坐标画出的散点图，以说明某个上调或下调的基因的状态

可以试运行该行代码

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

更加美观的火山图需要用到ggpubr包来实现

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')))##对数据进行筛选分组

table(df$g)

df$NAME=rownames(df)

ggscatter(df,x='logFC',y='v',size = 0.5,color = 'g')

还有一些具体的参数可以到官网去查看