[R语言可视化-精美图形绘制系列]--显著性箱线图

本期内容为[R语言可视化-精美图形绘制系列]--显著性箱线图

花费一点时间整理出来的，还是比较粗糙，主要的问题还是主题颜色，因此如有你好看的主题颜色，可以留言分享哦！！

是小杜好友的可以直接和小杜索要源代码和实例数据，我希望的是：交流（Communication）！！
欢迎来扰！！不要客气，要不你的加友费浪费了哦！

代码部分：
1. 此部分可看我们的个人喜好进行设置，是可选项。

## 自定义主题
top.mar=0.2
right.mar=0.2
bottom.mar=0.2
left.mar=0.2
## 合并上面的参数，将其合并成mythemel
mytheme1 <- theme(panel.background = element_blank(),
                  axis.ticks.length=unit(1.6,"mm"),
                  plot.margin=unit(x=c(top.mar,right.mar,botton.mar,left.mar),
                                   units="inches"))
#自定义主题2；
#隐藏纵轴，并对字体样式、坐标轴的粗细、颜色、刻度长度进行限定；
mytheme2<-theme_classic()+
  theme(text=element_text(family = "sans",colour ="gray30",size = 12),
        axis.line = element_line(size = 0.6,colour = "gray30"),
        axis.ticks = element_line(size = 0.6,colour = "gray30"),
        axis.ticks.length = unit(1.5,units = "mm"),
        plot.margin=unit(x=c(top.mar,right.mar,botton.mar,left.mar),
                         units="inches"))

2. 导入数据

librar(tidyr)
library(ggplot2)
library(ggthemr)
library(tidyverse)
library(dplyr)
library(ggpubr)
library(devEMF)
library(ggsignif)

##--
d2 <- read.csv("input_01.csv", header = T)
head(d2)
## 用one way anova计算 p value
pvalues <- sapply(d2$gene, function(x) {
  res <- aov(expr ~ stage, data = subset(d2, gene == x))
  summary(res)[[1]]$'Pr(>F)'[1] #
})
pv <- data.frame(gene = d2$gene, pvalue = pvalues)
##-----------------------------------------------------
gene      pvalue sigcode
1  mRNA_01 0.009434952      **
2  mRNA_01 0.009434952      **
3  mRNA_01 0.009434952      **
4  mRNA_01 0.009434952      **
5  mRNA_01 0.009434952      **

3. 绘图

ggplot(d2, aes(gene, expr, fill=stage)) + 
  geom_boxplot() + 
  geom_text(aes(gene, y=max(d2$expr) * 1.1, 
                label=paste("p = ",round(pvalue, 2))),
            data=pv, #
            inherit.aes=F) +
  xlab(NULL)+ylab("Relative expression (log2)")+
  ## 更改主题
  theme_classic()+ ##根据自己的喜好更改
  ## 添加颜色
  scale_fill_brewer(palette = "Dark2")+ # Accent, Dark2, Paired, Pastel1, Pastel2, Set1, Set2, Set3
## 更改字体大小
  theme(text = element_text(color = "black",size = 12),
    axis.text.x = element_text(color = "black",size = 12),
        axis.text.y = element_text(color = "black",size = 10))

往期文章（总汇）
01-[R语言可视化-精美图形绘制系列]--精美火山图
02-R语言可视化-精美图形绘制系列--柱状图
03-R语言可视化-精美图形绘制系列--功能富集分析
04-R语言可视化-精美图形绘制系列—多组GO富集可视化
05-[R语言可视化-精美图形绘制系列--堆积图]
06-[R语言可视化-精美图形绘制系列--组间相关性分析]
07-[R语言可视化-精美图形绘制系列]--Mental分析
08-[R语言可视化-精美图形绘制系列--复杂热图+两图渐变连线]-【转载】
09-[R语言可视化-精美图形绘制系列--桑基图（Sankey）]
10-[R语言可视化-精美图形绘制系列--柱状图误差线标记]
11-跟着NC学作图 | 柱状图与相关性图
12-[R语言可视化-精美图形绘制系列--GO、KEGG富集通路关联图]
13-[跟着“基迪奥生物学”作图]--截断图
--
小杜的生信筆記 ，主要发表或收录生物信息学的教程，以及基于R的分析和可视化（包括数据分析，图形绘制等）；分享感兴趣的文献和学习资料！