利用R软件包ggplot2绘制boxplot图

数据集：三个阶段卵母细胞（NSN，NSN-SN，SN）的核内组蛋白集落（分大小，large foci，small foci）的数量。一般而言，ggplot2需要的数据，每个样本为一行，样本信息需要有观测值，以及样本的属性信息等。

文件名：mydata2.txt

*****************************

number    oocyteType       fociType

18     NSN large

18     NSN large

18     NSN large

15     NSN large

25     NSN large

26     NSN large

26     NSN large

11     NSN large

24     NSN large

16     NSN large

21     NSN large

18     NSN large

24     NSN large

19     NSN large

12     NSN large

20     NSN large

11     NSN large

17     NSN large

14     NSN large

18     NSN large

14     NSN large

23     NSN large

22     NSN large

16     NSN large

17     NSN large

15     NSN large

16     NSN large

19     NSN large

23     NSN large

19     NSN large

14     NSN large

3       SN    large

6       SN    large

9       SN    large

0       SN    large

11     SN    large

1       SN    large

3       SN    large

1       SN    large

2       SN    large

2       SN    large

5       SN    large

2       SN    large

4       SN    large

3       SN    large

4       SN    large

7       SN    large

7       SN    large

5       SN    large

15     NSN-SN    large

7       NSN-SN    large

2       NSN-SN    large

13     NSN-SN    large

15     NSN-SN    large

9       NSN-SN    large

8       NSN-SN    large

11     NSN-SN    large

13     NSN-SN    large

9       NSN-SN    large

1       NSN-SN    large

19     NSN-SN    large

8       NSN-SN    large

11     NSN-SN    large

13     NSN-SN    large

35     NSN small

23     NSN small

12     NSN small

15     NSN small

55     NSN small

28     NSN small

21     NSN small

38     NSN small

17     NSN small

18     NSN small

18     NSN small

37     NSN small

19     NSN small

14     NSN small

41     NSN small

21     NSN small

30     NSN small

7       NSN small

9       NSN small

9       NSN small

8       NSN small

29     NSN small

8       NSN small

15     NSN small

8       NSN small

13     NSN small

14     NSN small

6       NSN small

12     NSN small

5       NSN small

23     NSN small

8       SN    small

10     SN    small

6       SN    small

8       SN    small

9       SN    small

3       SN    small

8       SN    small

3       SN    small

9       SN    small

5       SN    small

9       SN    small

7       SN    small

7       SN    small

10     SN    small

10     SN    small

8       SN    small

8       SN    small

6       SN    small

12     NSN-SN    small

18     NSN-SN    small

16     NSN-SN    small

20     NSN-SN    small

15     NSN-SN    small

10     NSN-SN    small

15     NSN-SN    small

12     NSN-SN    small

11     NSN-SN    small

14     NSN-SN    small

16     NSN-SN    small

13     NSN-SN    small

9       NSN-SN    small

14     NSN-SN    small

11     NSN-SN    small

*****************************

（1）读入数据

mydata3<-read.table("mydata3.txt",header=T,sep="\t")

（2）初步绘图

> library(ggplot2)

> ggplot(data=mydata3, aes(x=oocyteType,y=number))+geom_boxplot(aes(fill=fociType))

得如下之图：

（3）加横轴纵轴标签

> p <- ggplot(data=mydata3, aes(x=oocyteType,y=number))+geom_boxplot(aes(fill=fociType))

> p1 <- p+labs(x = "oocyte developmental stage")

> p2 <- p1+labs(y = "XXX foci number")

> p2

（4）修改坐标轴标签之属性

> p3<-p2 + theme(axis.title.x=element_text(color="red", size="15", face="bold"))

> p4<-p3 + theme(axis.title.y=element_text(color="red", size="15", face="bold"))

> p4

（5）修改背景

> p4 + theme_bw()

背景应该在第一步开始设置，中途设置背景，将上面设置的标签掩盖了。

背景修改方式，可见

> ?theme_bw()

（6）更改Figure Legend位置

> p5 <- p4 + theme_bw()

> p5+theme(legend.position=c(0.8,0.8))

> p5 +theme(axis.title.x=element_text(color="red", size="15", face="bold"))+ theme(axis.title.y=element_text(color="red", size="15", face="bold"))

（7）修改颜色

> p6 <- p5 + scale_fill_manual(values=c("red","blue"))

> p6

（8）修改坐标轴刻度标签颜色大小

> p6<-p6+theme(axis.text.x=element_text(size=13,color="black"))

> p6<-p6+theme(axis.text.y=element_text(size=13,color="black"))

> p6

（9）修改legend

> p6+guides(fill=guide_legend(title=NULL))

（10）加图题

> p7<-p6+guides(fill=guide_legend(title=NULL))

> p7+annotate("text",label="my boxplot plot",x=2,y=50,color="Green",size=8)