使用R语言的正确姿势，R包干货奉献

生物信息学习的正确姿势

NGS系列文章包括NGS基础、在线绘图、转录组分析 （Nature重磅综述|关于RNA-seq你想知道的全在这）、ChIP-seq分析 （ChIP-seq基本分析流程）、单细胞测序分析 (重磅综述：三万字长文读懂单细胞RNA测序分析的最佳实践教程)、DNA甲基化分析、重测序分析、GEO数据挖掘（典型医学设计实验GEO数据分析 (step-by-step)）、批次效应处理等内容。

本文来源：R友舍公众号

先思考，后动手，能让你减少很多无谓的功夫，加快完成任务的效率。按照一定的逻辑思路来行事，更是如此，比如按照一定的通用流程在数据科学项目中使用R语言。

在实际工作中，每个数据科学项目各不相同，但基本都遵循一定的通用流程。具体如下：

• 数据导入

• 数据整理

• 数据可视化

• 数据转换

• 数据建模

• 统计推断

• 报告展示

• 自动化分析

下面列出每个步骤最有用的一些R包。

1. 数据导入

以下R包主要用于数据导入和保存数据：

• feather：一种快速，轻量级的文件格式；在R和python上都可使用

• readr：实现表格数据的快速导入

• readxl：读取Microsoft Excel电子表格数据

• openxlsx：读取Microsoft Excel电子表格数据

• googlesheets：读取google电子表格数据

• haven：读取SAS，SPSS和Stata统计软件格式的数据

• httr：从网站开放的API中读取数据

• rvest：网页数据抓取包

• xml2：读取HTML和XML格式数据

• webreadr：读取常见的Web日志格式数据

• DBI：数据库管理系统通用接口包

• RMySQL：用于连接MySQL数据库的R包

• RPostgres：用于连接PostgreSQL数据库的R包

• bigrquery：用于连接Google BigQuery的R包

• PivotalR：用于读取Pivitol(Greenplum)和HAWQ数据库中的数据

• dplyr：提供了一个访问常见数据库的接口

• data.table：data.table包的fread()函数可以快速读取大数据集

• git2r：用于访问git仓库

2. 数据整理

以下R包主要用于数据整理，以便于你后续建模分析：

• tidyr：用于整理表格数据的布局

• dplyr：用于将多个数据表连接成一个整齐的数据集

• purrr：函数式编程工具，在做数据整理时非常有用

• broom：用于将统计模型的结果整理成数据框形式

• zoo：定义了一个名zoo的S3类型对象，用于描述规则的和不规则的有序的时间序列数据

3. 数据可视化

以下R包用于数据可视化：

• ggplot2及其扩展包：ggplot2包提供了一个强大的绘图系统，并实现了以下扩展 (ggplot2高效实用指南 (可视化脚本、工具、套路、配色))

• ggthemes：提供扩展的图形风格主题

• ggmap：提供Google Maps、Open Street Maps等流行的在线地图服务模块

• ggiraph：绘制交互式的ggplot图形

• ggstance：实现常见图形的横向版本

• GGally：绘制散点图矩阵

• ggalt：添加额外的坐标轴、geoms等

• ggforce：添加额外geoms等

• ggrepel：用于避免图形标签重叠

• ggraph：用于绘制网络状、树状等特定形状的图形

• ggpmisc：光生物学相关扩展

• geomnet：绘制网络状图形

• ggExtra：绘制图形的边界直方图

• gganimate：绘制动画图

• plotROC：绘制交互式ROC曲线图

• ggspectra：绘制光谱图

• ggnetwork：网络状图形的geoms

• ggradar：绘制雷达图

• ggTimeSeries：时间序列数据可视化

• ggtree：树图可视化

• ggseas：季节调整工具

• lattice：生成栅栏图

• rgl：交互式3D绘图

• ggvis：交互式图表多功能系统

• htmlwidgets：一个专为R语言打造的可视化JS库

• leaflet：绘制交互式地图

• dygraphs：绘制交互式时间序列图

• plotly：交互式绘图包

• rbokeh：用于创建交互式图表和地图

• Highcharter：绘制交互式Highcharts图

• visNetwork：绘制交互式网状图

• networkD3：绘制交互式网状图 (network3D 交互式网络生成)

• d3heatmap：绘制交互式热力图

• DT：用于创建交互式表格

• threejs：绘制交互式3d图形和地球仪

• rglwidget：绘制交互式3d图形

• DiagrammeR：绘制交互式图表

• MetricsGraphics：绘制交互式MetricsGraphics图

• rCharts：提供了对多个javascript数据可视化库(highcharts/nvd3/polychart)的R封装

• coefplot：可视化统计模型结果

• quantmod：可视化金融图表

• colorspace：基于HSL的调色板

• viridis：Matplotlib viridis调色板

• munsell：Munsell调色板

• RColorBrewer：图形调色板

• igraph：用于网络分析和可视化

• latticeExtra：lattice绘图系统扩展包

• sp：空间数据工具

• amplicon： 扩增子分析

• YSX：一步可视化
  

4. 数据转换

以下R包用于将数据转换为新的数据类型：

• dplyr：一个用于高效数据清理的R包

• magrittr：一个高效的管道操作工具包

• tibble：高效的显示表格数据的结构

• stringr：一个字符串处理工具集

• lubridate：用于处理日期时间数据

• xts：xts是对时间序列数据(zoo)的一种扩展实现，提供了时间序列的操作接口。

• data.table：用于快速处理大数据集

• vtreat：一个对预测模型进行变量预处理的工具

• stringi：一个快速字符串处理工具

• Matrix：著名的稀疏矩阵包

5. 统计建模与推断

下述R包是统计建模最常用的几个R包，其中的一些R包适用于多个主题：

• car：提供了大量的增强版的拟合和评价回归模型的函数

• Hmisc：提供各种用于数据分析的函数

• multcomp：参数模型中的常见线性假设的同时检验和置信区间计算，包括线性、广义线性、线性混合效应和生存模型

• pbkrtest：用于线性混合效应模型的参数Bootstrap检验

• MatrixModels：用于稠密矩阵和稀疏矩阵建模

• mvtnorm：用于计算多元正态分布和t分布的概率，分位数，随机偏差等

• SparseM：用于稀疏矩阵的基本线性代数运算

• lme4：利用C 矩阵库 Eigen进行线性混合效应模型的计算。

• broom：将统计模型结果整理成数据框形式

• caret：一个用于解决分类和回归问题的数据训练综合工具包

• glmnet：通过极大惩罚似然来拟合广义线性模型

• gbm：用于实现随机梯度提升算法

• xgboost：全称是eXtremeGradient Boosting。是Gradient Boosting Machine的一个c 实现。目前已制作了xgboost工具的R语言接口

• randomForest：提供了用随机森林做回归和分类的函数

• ranger：用于随机森林算法的快速实现

• h2o：H2O是0xdata的旗舰产品，是一款核心数据分析平台。它的一部分是由R语言编写的，另一部分是由Java和Python语言编写的；用户可以部署H2O的R程序安装包，之后就可以在R语言环境下运行了

• ROCR：通过绘图来可视化分类器的综合性能

• pROC：用于可视化，平滑和对比ROC曲线

6. 沟通交流

以下R包用于实现数据科学结果的自动化报告，以便于你跟人们进行沟通交流：

• rmarkdown：用于创建可重复性报告和动态文档

• knitr：用于在PDF和HTML文档中嵌入R代码块

• flexdashboard：基于rmarkdown，可以轻松的创建仪表盘

• bookdown：以R Markdown为基础，用于创作书籍和长篇文档 (别人的电子书，你的电子书，都在bookdown)

• rticles：提供了一套R Markdown模板

• tufte：用于实现Tufte讲义风格的R Markdown模板

• DT：用于创建交互式的数据表

• pixiedust：用于自定义数据表的输出

• xtable：用于自定义数据表的输出

• highr：用于实现R代码的LaTeX或HTML格式输出

• formatR：通过tidy_source函数格式化R代码的输出

• yaml：用于实现R数据与YAML格式数据之间的通信

7. 自动化分析

以下R包用于创建自动化分析结果的数据科学产品：

• shiny：一个使用R语言开发交互式web应用程序的工具

• shinydashboard：用于创建交互式仪表盘

• shinythemes：给出了Shiny应用程序的常用风格主题

• shinyAce：为Shiny应用程序开发者提供Ace代码编辑器

• shinyjs：用于在Shiny应用程序中执行常见的JavaScript操作

• miniUI：提供了一个UI小部件，用于在R命令行中集成交互式应用程序

• shinyapps.io：为创建的Shiny应用程序提供托管服务

• Shiny ServerOpen Source：为Shiny应用程序提供开源免费的服务器

• Shiny Server Pro：为企业级用户提供一个Shiny应用程序服务器

• rsconnect：用于将Shiny应用程序部署到shinyapps.io

• plumber：用于将R代码转化为一个web API

• rmarkdown：用于创建可重复性报告和动态文档

• rstudioapi：用于安全地访问RStudio IDE的API

8. 程序开发

以下这些包主要用于开发自定义的R包：

• RStudio Desktop IDE：R的IDE。大家都懂，不用解释。

• RStudio Server Open Source：开源免费的RStudio服务器

• RStudio Server Professional：商业版RStudio服务器

• devtools：一个让开发R包变得简单的工具集

• packrat：创建项目的特定库，用于处理包的版本问题，增强代码重现能力。

• drat：一个用于创建和使用备选R包库的工具

• testthat：单元测试，让R包稳定、健壮，减少升级的痛苦

• roxygen2：通过注释的方式，生成文档，远离Latex的烦恼

• purrr：一个用于 提供函数式编程方法的工具

• profvis：用于可视化R代码的性能分析数据

• Rcpp：用于实现R与C 的无缝整合

• R6：R6是R语言的一个面向对象的R包，可以更加高效的构建面向对象系统。

• htmltools：用于生成HTML格式输出

• nloptr：提供了一个NLopt非线性优化库的接口

• minqa：一个二次近似的优化算法包

• rngtools：一个用于处理随机数生成器的实用工具

• NMF：提供了一个执行非负矩阵分解的算法和框架

• crayon：用于在输出终端添加颜色

• RJSONIO：rjson是一个R语言与json进行转的包，是一个非常简单的包，支持用 C类库转型和R语言本身转型两种方式

• jsonlite：用于实现R语言与json数据格式之间的转化

• RcppArmadillo：提供了一个Armadillo C Library（一种C 的线性代数库）的接口

9. 实验数据

以下R包给出了案例实战过程中可用的训练数据集：

• babynames：包含由美国社会保障局提供的三个数据集

• neiss：2009-2014年期间提供给美国急诊室的所有事故报告样本数据

• yrbss：美国疾病控制中心2009-2013年期间青年危险行为监测系统数据

• USAboundaries：2011年全年休斯顿机场的所有航班数据

• rworldmap：国家边界数据

• usdanutrients：美国农业部营养数据库

• fueleconomy：美国环保署1984-2015年期间的燃油经济数据

• nasaweather：包含了一个覆盖中美洲的非常粗糙的24*24格地理位置和大气测量数据。

• mexico-mortality：墨西哥死亡人数数据

• data-movies和ggplotmovies：来自互联网电影数据库imdb.com的数据

• pop-flows：2008年全美人口流动数据

• data-housing-crisis：经过清洗后的2008美国房地产危机数据

• gun-sales：纽约时报提供的有关枪支购买的每月背景调查统计分析数据

• stationaRy：从成千上万个全球站点收集到的每小时气象数据

• gapminder：摘自Gapminder的数据

• janeaustenr：简·奥斯丁小说全集数据

本文直接来源：R语言统计与绘图 公众号。

往期精品(点击图片直达文字对应教程)

后台回复“生信宝典福利第一波”或点击阅读原文获取教程合集