前几日偶然从科学网上看到了一副特别漂亮的统计图形,如下图(摘自博文:http://blog.sciencenet.cn/blog-397960-456174.html)所示:
感觉美翻了都,刚好在学习R,又遇到课题组有一些数据要处理,就花了大约一天的时间对其制作过程进行分析,现将其绘制过程简述如下:
1.数据准备
假设共有三组数据,分别为A,B,C,其数值如下:
A(
0.2114 | 0.201 | 0.001 | 0.0416 | 0.1238 | 0.1195 | 0.1928 | 0.2678 |
)
B(
0.2542 | 0.1367 | 0.1559 | 0.2009 | 0.1803 |
)
C(
0.1064 | 0.4412 | 0.1659 | 0.1596 | 0.1418 | 0.0975 | 0.2365 | 0.1755 |
)
首先需要按照下述格式编辑文档(data.txt)
0.2114 A
0.201 A
0.001 A
...
0.2542 B
0.1367 B
0.1559 B
...
0.1064 C
0.4412 C
...
将编辑好的文本文档文件保存于C:UsersAdministrator我的文档目录下
2.软件包准备
为完成上述作图,需要在安装R的基础上安装ggplot2软件包
打开R,
>install.packages('ggplot2')
这时R会提醒你选择要下载的镜像,本人人在海南,实验发现China(Xiamen)镜像速度最快。
待ggplot安装完成之后,便可以开始图形的绘制。
具体命令如下:
>data<-read.table("data.txt") #读取所有分析的文件
library(ggplot2) #调用ggplot软件包
p<-ggplot(data,aes(x=factor(V2),y=V1))
p+geom_boxplot(col=“blue",pch=16,cex=1)+geom_point(position="jitter"col=2,pch=16,cex=1)
作图完成,具体使用者可根据自身需要对图进行进一步优化。
效果图如下所示:
http://blog.sciencenet.cn/blog-1334016-796140.html
上一篇:关于Linux系统cd-hit软件和probcons软件的解析
下一篇:The question and solution of GEOquery install