最近在整理数据,写论文,需要画图。厌烦每次画图都要重新翻代码(几天不用就忘干净了),所以打算将每次画图的代码以及数据格式保存下来,以便以后查阅。这一次是画箱线图,即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")