最全箱线图boxplot指南

前面介绍了用melt函数处理数据，得到的结果就可以画boxplot.

数据处理

原始数据为一个比较组的值，如下图

首先删除缺失值,并使用melt函数处理数据格式

data <- na.omit(data)
newdata <- melt(data,measure.vars = c("WT","control"),variable.name = "Sample")

处理后得到

画图

library(ggplot2)

ggplot(newdata,aes(y=value,x=Sample))+

stat_boxplot(geom="errorbar",width=0.15,aes(color=Sample))+

  geom_boxplot(aes(fill=Sample),outlier.shape=NA)+
  
  scale_fill_manual(values=c("#3D59AB", "#E3170D"))+
  
  theme(axis.text=element_text(size=15,color='black',face='bold'),
        axis.title=element_text(size=20,color='black',face='bold'),
        legend.title=element_text(size=18,color='black',face='bold'),
        legend.text=element_text(size=12,color='black',face='bold'),
        panel.background=element_rect(fill="white",colour="black"),
        panel.grid.major=element_line(color="gray"),
        panel.grid.minor=element_line(color="gray"))+ylim(0.5,4)

未完待续。。。