R语言的重要性我也就不再赘述了,它不仅在生物信息数据处理中发挥着重要作用,也是其它主流数据处理人士的首选工具。现在非常多的半路出家自学生物信息学的小伙伴必须学而且有可能学的就是R,所以写一个R的系统性入门指导是非常有必要的。这本来应该是我下面的生信分析人员如何系统入门编程语言的姊妹篇的,但是因为时隔太久,我的感悟可能发生了变化,所以这个R跟前面的两个看起来总结指引模式不太一样的。
生信分析人员如何系统入门perl?
生信分析人员如何系统入门linux?
我作为老一辈的生信工程师,所以喜欢perl一点,排斥python,其实呢,我也稍微看过一些python的语法,个人认为R和python几乎是一模一样的。R的特点就是内置了大量的函数,基本上你认识的英文单词都可以是一个函数,即使不是,你也可以自定义为函数。搞清楚了函数和变量,就可以看懂大部分的R代码了。
下面是生信菜鸟团QQ群管理员赵云对这3种编程语言的心得体会!
python跟perl都是高级语言, 两个开发的目的不同, perl更面向过程一些,优势是严谨,快。 python主流面向对象编程, 这个跟R类似, 数据结构等方面有些不同,但可以互相调用。 实际上以上三者之间可以互相调用部分功能。python的语法并不是很严谨,个人感觉,越偏向自然语言的编程语言越通俗但不严谨,以上,是跟C比较的。
R本身起源于S语言,是主要针对统计的, 也是面向对象的, 本质上,,是把一个比excel功能强大的软件归零化成了命令行吧.excel高级应用也是要编程的, 所以R的初级应用可以当成是没有用户交互界面的excel,细心一点, 把示例代码都打对,当功能强大但不好使版的excel吧, 这样至少心理上不会畏难跟抵触.
内部集成的越多,用户需要做的越少, 你用C画个图累死你, 用python得写几行, R一行就行了!
PS:菜鸟发言,如有误导,概不负责!
接下来是我写的了:
必须要做的:
你首先肯定得安装R语言软件,然后推荐也安装一些R studio,在R studio里面写代码会比较方便。
下载R语言的软件:
•下载Rstudio这个R编辑器:
然后我就要强调help函数了,你必须要把help函数用一百次以上,不然你不可能入门的!你必须要自学R语言基础,或者看书,或者看视频,或者有人手把手教你,书的话,我推荐:
《R in Action》
《The Art of_R Programming》
R的特性就是有着大量的包,所以你必须学会安装包:
安装包 install.packages(" xxxxxx ")
加载包 library( xxxxx )
查看包的帮助文档help("xxxxx") 或?xxxxx
获取当前工作区间getwd()
更改工作区间 setwd( "xxxxxx")
清除当前对象rm()
安装包你一定会遇到错误,请参考: R包终极解决方案! ,还有R的包(package)
上面的书籍都会提供一些简单的测试代码,你跟着傻瓜式的敲代码就好,但是实践的过程中,请务必注意一些英文单词file 文件路径 Description简述 Usage用法 Arguments参数 Details详细 value 数值 Examples例子 header 表标题 logical_value 逻辑值 delimiter 分隔符 object 对象 col列 row 行 vector向量 dimensions维度 data数据 。
首先你要明白的是变量:
向量和因子,我懒得说了。(向量特简单,没什么好说的,因子太复杂了,我说不清楚,你们慢慢理解。)
数据框–就像我们的表格,第一行就是每一列的名字,我们称之为字段,或者变量名.那么对应每列下面的数据就叫做记录或者观测.用data.frame( 字段1,字段2,…. )创建
列表–与数据框类似,区别就是每一列向量类型和长度可以不一致.用list( 字段1, 字段2,….. )创建
数组–其形式就像我们玩的模方,每一个面都是一个矩阵数据,用array(数据,各维度的最大值,各维度的名称)
然后要了解对这些变量的基础操作函数:
变量怎么来,对它们处理什么?
我们处理生物信息学数据一般很少会手动创建这些对象,都是从文本里面读取,比如kegg数据库文件,差异分析结果,RNA-seq的表达量矩阵,但是读入之后,我们的重点就是知道它们变成了什么,该如何去一步步的转换它们。
数据的特性函数也必须要知道,无非就是一些英文单词而已,你经常的玩一下,就慢慢的熟练了。str,class,names,row.names,col.names,length,unique,view,min,max,summay,table
可视化你的变量:
了解了R里面的基础变量和对象,也学会了对它们进行简单的转换,接下来就可以尝一尝R的甜头了,对任何数据都可以可视化,简简单单的就可以画一大堆的图。
plot,boxplot,barplot,pie,hist,pair,它们每个绘图函数都有自己要求的输入数据,特定的可视化结果,请务必在还没熟练使用之前help一下它们,自己主动查看它们好玩的地方,好好自学。
dev.new()新建画板
plot()绘制点线图,条形图,散点图.
barplot( ) 绘制条形图
dotchart( ) 绘制点图
pie( )绘制饼图.
pair( )绘制散点图阵
boxplot( )绘制箱线图
hist( )绘制直方图
scatterplot3D( )绘制3D散点图.
低级绘图函数:
par() 可以添加很多参数来修改图形
title( ) 添加标题
axis( ) 调整刻度
rug( ) 添加轴密度
grid( ) 添加网格线
abline( ) 添加直线
lines( ) 添加曲线
text( ) 添加标签
legend() 添加图例
它们还有一系列的绘图参数,坐标轴,图例,颜色,性状,大小,空白,布局,非常繁琐,想掌握,花费的时间会非常多,但是很多人直接跳到ggplot的绘图世界了,不想搞那么多底层绘图代码。
但是我看过一个底层R绘图集大成者,就是Combining gene mutation with gene expression data improves outcome prediction in myelodysplastic syndromes文章的作者的github里面有。
但是对大部分人来说,生信的绘图,都是有套路的,其实都被别人包装成函数了,做好数据,一个函数就出了所有复杂的图。比如热图,cluster等等
数据对象的一些高级操作
前面我们对向量,数据框,数组,列表都了解了,也知道如何查看数据的特性,但是要进行高级转换,就需要一些时间来学习apply系列函数,aggregate,split等函数的用法,这是一个分水岭,用好了你就算是R入门了。也可以用一些包,比如reshape2,dplyr
当然,R里面的字符串对象是另外完全不一样的操作模式,建议大家自行搜索学习。
R的bioconductor世界
这个是生物信息学特有的,也是为什么我要求搞生物信息学数据处理的人必须学习R,就是为了应用大量的bioconductor包。在这里面所有的对象都不在是基础的向量,数据框,数组,列表了,而是S3,S4对象,这个高级知识点我就不推荐了,你学会了前面的东西,就有了自己的学习经验了,后面的分分钟就搞定了。(其实你永远也搞不定的)
每学一个bioconductor的包,都是自己R水平的提升,大家可以参考我的博客:http://www.bio-info-trainee.com/tag/bioconductor 我就是这样学习过来的。我还创建了bioconductor中国这个社区,可惜效果不好,有志者可以继续联系我,我们看看有没有可能做起来。
当然R肯定不只是应用在生物信息学啦,其实它在非常多的地方都有应用,尤其是金融和地理,在如何一个方向学习R,就不仅仅是R本身的语法了,你需要学习的东西太多了,我简单列出几个我接触过的方向吧:
统计,科学计算,数据挖掘,文本挖掘,基础绘图,ggplot绘图,高级编程,都有着丰富的书籍和视频资料,尤其推荐炼数成金的R七种武器系列。(全套视频很容易找到)
《A Handbook of Statistical Analyses_Using_R》
《Modern Applied Statistics With S》
《Introduction to Scientific Programming and Simulation Using R》
《Mastering Scientific Computing with R》
《Practical Data Science with R》
《Data Mining explain using R》
《ggplot2 Elegant Graphics for Data Analysis》
《R Graphics Cookbook》
《R Cookbook》
《R in a Nutshell》
《R Programming for Bioinformatics》
《software for data analysis programming with R》
看完这些,你就是R大神了,当然,前提是你看懂了也会灵活应用。
# 网络资源#
- [R语言官方站](http://www.r-project.org/)
- [R-blogger](http://www.r-bloggers.com/)
- [R语言资源汇总](https://github.com/qinwf/awesome-R)
- [R语言搜索引擎](http://www.rseek.org/)
- [R函数在线帮助](http://www.rdocumentation.org/)
- [关于R的问答网站](http://stackoverflow.com/questions/tagged/r)
- [一个入门级的R在线教程](http://tryr.codeschool.com/)
- [交互式的R在线教程](https://www.datacamp.com)
- [统计之都](http://cos.name/)
- [我的博客](http://xccds1977.blogspot.com)
- [美国计算机世界杂志提供的R语言初学者入门资料](http://www.computerworld.com/s/article/9239625/Beginner_s_guide_to_R_Introduction)
- 各种cheatsheet
最后借鉴一个生物信息学习方法:
知识和耐心,是成为强者的唯一方法。
- 通过阅读来学习。
包括了阅读经典的教材、代码、论文、学习公开课。
- 通过牛人来学习。
包括同行的聚会、讨论、大牛的博客、微博、twitter、RSS。
- 通过练习来学习。
包括代码练习题、参加kaggle比赛、解决实际工作中的难题。
- 通过分享来学习。
包括自己写笔记、写博客、写书、翻译书,和同伴分享交流、培训新人。