使用R包ggplot2画箱线图（boxplot）

  最近在整理数据，写论文，需要画图。厌烦每次画图都要重新翻代码（几天不用就忘干净了），所以打算将每次画图的代码以及数据格式保存下来，以便以后查阅。这一次是画箱线图，即boxplot。

数据输入格式如下：

Expression_level        Organ_type      Count
unexpressed     vegetative_organs       699
unexpressed     vegetative_organs       700
unexpressed     vegetative_organs       660
unexpressed     vegetative_organs       666
unexpressed     vegetative_organs       672
unexpressed     vegetative_organs       731
unexpressed     vegetative_organs       725
unexpressed     vegetative_organs       712
unexpressed     vegetative_organs       720
unexpressed     reproductive_organs     624
unexpressed     reproductive_organs     617
unexpressed     reproductive_organs     566
unexpressed     reproductive_organs     693
unexpressed     reproductive_organs     634
unexpressed     reproductive_organs     697
unexpressed     reproductive_organs     674
unexpressed     reproductive_organs     703
unexpressed     reproductive_organs     711
unexpressed     reproductive_organs     556
unexpressed     reproductive_organs     743
unexpressed     reproductive_organs     760
low_expressed   vegetative_organs       1187
low_expressed   vegetative_organs       1182
low_expressed   vegetative_organs       1217
low_expressed   vegetative_organs       1220
low_expressed   vegetative_organs       1213
low_expressed   vegetative_organs       1145
low_expressed   vegetative_organs       1162
low_expressed   vegetative_organs       1178
low_expressed   vegetative_organs       1159
low_expressed   reproductive_organs     1243
low_expressed   reproductive_organs     1245
low_expressed   reproductive_organs     1286
low_expressed   reproductive_organs     1149
low_expressed   reproductive_organs     1244
low_expressed   reproductive_organs     1144
low_expressed   reproductive_organs     1175
low_expressed   reproductive_organs     1148
low_expressed   reproductive_organs     1131
low_expressed   reproductive_organs     1298
low_expressed   reproductive_organs     1131
low_expressed   reproductive_organs     1114
medium_expressed        vegetative_organs       20
medium_expressed        vegetative_organs       26
medium_expressed        vegetative_organs       32
medium_expressed        vegetative_organs       23
medium_expressed        vegetative_organs       23
medium_expressed        vegetative_organs       33
medium_expressed        vegetative_organs       23
medium_expressed        vegetative_organs       19
medium_expressed        vegetative_organs       25
medium_expressed        reproductive_organs     31
medium_expressed        reproductive_organs     42
medium_expressed        reproductive_organs     45
medium_expressed        reproductive_organs     51
medium_expressed        reproductive_organs     21
medium_expressed        reproductive_organs     50
medium_expressed        reproductive_organs     43
medium_expressed        reproductive_organs     40
medium_expressed        reproductive_organs     49
medium_expressed        reproductive_organs     41
medium_expressed        reproductive_organs     21
medium_expressed        reproductive_organs     19
 high_expressed vegetative_organs       6
 high_expressed vegetative_organs       4
 high_expressed vegetative_organs       3
 high_expressed vegetative_organs       3
 high_expressed vegetative_organs       4
 high_expressed vegetative_organs       3
 high_expressed vegetative_organs       2
 high_expressed vegetative_organs       3
 high_expressed vegetative_organs       8
 high_expressed reproductive_organs     14
 high_expressed reproductive_organs     8
 high_expressed reproductive_organs     15
 high_expressed reproductive_organs     19
 high_expressed reproductive_organs     13
 high_expressed reproductive_organs     21
 high_expressed reproductive_organs     20
 high_expressed reproductive_organs     21
 high_expressed reproductive_organs     21
 high_expressed reproductive_organs     17
 high_expressed reproductive_organs     17
 high_expressed reproductive_organs     19

R 代码

library(ggplot2)

Data<-read.table("/Users/mashengwei/PycharmProjects/common/1.txt", head=T)

p<-ggplot(data=Data, aes(x=Expression_level,y=Count))+geom_boxplot(aes(fill=Organ_type))

p+ facet_wrap(~ Expression_level, scales="free")

png("/Users/mashengwei/PycharmProjects/common/1.png") ＃保存图片

图片展示

添加坐标轴标题

p+ facet_wrap(~ Expression_level, scales="free") + xlab(" Expression level of Triticeae-specific genes") + ylab("Count of Triticeae-specific genes")