ggplot2火山图展示RNAseq差异表达分析结果

结果图

完整代码

DEGs<-read.csv("volcano_plot_example_data.csv",header=T,stringsAsFactors = F)
dim(DEGs)
df<-DEGs[-log10(DEGs$padj)>25,]
dim(df)
library(ggplot2)
library(ggrepel)
ggplot(DEGs,aes(x=log2FoldChange,y=-log10(padj)))+
  geom_point(aes(color=change),size=2.5,alpha=1,na.rm = T)+
  geom_hline(yintercept =-log10(0.05),color="#990000",
             linetype="dashed")+
  geom_vline(xintercept = -2,color="#990000",linetype="dashed")+
  geom_vline(xintercept = 2,color="#990000",linetype="dashed")+
  theme_bw(base_size = 14)+
  scale_color_manual(values=c("red","#00B2FF","orange"))+
  xlab(expression(log[2]("FoldChange")))+
  ylab(expression(-log[10]("padj")))+
  theme(legend.title = element_blank())+
  ggtitle(label = "Volcano Plot", 
          subtitle = "Colored by fold-change direction")+
  geom_label_repel(data=df,aes(x=log2FoldChange,
                          y=-log10(padj),
                          label=name,fill=change),
                   color="white",size=5)+
  scale_fill_manual(values=c("red","orange"))+
  guides(fill=F)

需要示例数据的可以给我的公众号留言哈！

ggpubr包做火山图https://rpkgs.datanovia.com/ggpubr/reference/diff_express.html#examples

找到了一些使用R语言分析新型冠状病毒数据的例子

https://www.statsandr.com/blog/top-r-resources-on-covid-19-coronavirus/

Top 100 R resources on Novel COVID-19 Coronavirus

这个链接收集了很多使用R语言分析新型冠状病毒的例子，其中有很多shiny应用，有空的话多看看这些例子的代码。

https://www.r-bloggers.com/facts-about-coronavirus-disease-2019-covid-19-in-5-charts-created-with-r-and-ggplot2/

Facts About Coronavirus Disease 2019 (COVID-19) in 5 Charts created with R and ggplot2

这个链接有5幅漂亮的图片，但是没有发现代码和数据，抽空再找找