環(huán)境初始化

rm(list=ls())
Sys.setenv(LANGUAGE = "en") #顯示英文報(bào)錯(cuò)信息
options(stringsAsFactors = FALSE) #禁止chr轉(zhuǎn)成factor
options(BioC_mirror="https://mirrors.ustc.edu.cn/bioc/") 
options("repos" = c(CRAN="https://mirrors.tuna.tsinghua.edu.cn/CRAN/")) 
### 安裝包
cranpackage<-c("ggplot2","ggpubr")


for (pkg in cranpackage){
  if (! require(pkg,character.only=T) ) {
    install.packages(pkg,ask = F,update = F)
    require(pkg,character.only=T)
  }
}


# library(ggplot2)
# library(ggpubr)

讀入差異分析后結(jié)果

load("04_火山圖/volcano.Rdata")#通過(guò)將數(shù)據(jù)單獨(dú)存儲(chǔ)為R可以直接讀取的Rdata文件，方便調(diào)用

data<-deg1
data$gene <- rownames(data)#在data變量中新創(chuàng)建一個(gè)列，名稱為gene
colnames(data)
ml<-"adj.P.Val"
# ml<-"P.Value"
up<-length(data$gene[data[,ml]<0.05&data$logFC > 0.05])#注意，這種大小比較返回的是TRUE和FALSE,length()求的是長(zhǎng)度
down<-length(data$gene[data[,ml]<0.05&data$logFC< c(-0.05) ])#adj.P.Val小于0.05而且logFC< c(-0.05)

畫(huà)圖

ggplot(data=data, aes(x=logFC, y =-log10(adj.P.Val))) +
  ## 畫(huà)出變化倍數(shù)少的點(diǎn)
  geom_point(data=subset(data,abs(data$logFC) <= 0.05),color="#8B8B83",alpha=0.3) +
  ## 畫(huà)出logFC > 0.05的點(diǎn)
  geom_point(data=subset(data,data[,ml]<0.05 & data$logFC > 0.05),color="#CD5555",alpha=0.4) +
  ## 畫(huà)出1 < logFC < 2的點(diǎn)
  # geom_point(data=subset(data,data[,ml]<0.05 & data$logFC > 1&data$logFC < 2),aes(size=abs(logFC)),color="#FA8072",alpha=0.4) +
  ## 畫(huà)出logFC < -2的點(diǎn)
  geom_point(data=subset(data,data[,ml]<0.05 & data$logFC < c(-0.05)),color="darkgreen",alpha=0.4) +
  ## 畫(huà)出-2< logFC < -1的點(diǎn)
  # geom_point(data=subset(data,data[,ml]<0.05 & data$logFC >c(-2)& data$logFC <c(-1)),aes(size=abs(logFC)),color="#6fae6f",alpha=0.4) +
  ## 畫(huà)豎線
  geom_vline(xintercept = c(0.05,-0.05),lty=2,lwd=0.6,alpha=0.8)+
  geom_hline(yintercept = c(-log10(0.05)),lty=2,lwd=0.6,alpha=0.8)+
  ## 主題(這里面不設(shè)置主題)，ylim是人為規(guī)定的
  theme_bw()+
  ylim(-2,80)+
  scale_x_continuous(breaks=seq(-1, 1, 0.5)) +
  theme(panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(),
        axis.text.x = element_text(size = 10,colour = "black"),
        axis.text.y = element_text(size = 10,colour = "black"),
        axis.title.x = element_text(size = 10,colour = "black"),
        axis.title.y = element_text(size = 10,colour = "black"),
        axis.line = element_line(colour = "black"))+
  ## 修改坐標(biāo)軸
  labs(x="Log2 (fold change)",y="-log10 (q-value)")+
  ## 去掉注釋
  theme(legend.position='none')
ggsave("04_火山圖/vocalno.pdf",width = 4,height = 4.5)