R语言交互式可视化包CanvasXpress

作者简介

刘永鑫，博士。2008年和2011年毕业东北农业大学微生物学和作物遗传育种专业。2014年在中科院遗传发育所获生物信息学博士学位，2016年博士后出站留所工作，任宏基因组学实验室工程师，目前主要研究方向为宏基因组学数据分析方法、培养组学方法优化。2017年7月创办“宏基因组”公众号，目前关注人数1.5万，累计阅读超百万。

公众号：宏基因组(meta-genome)

往期回顾

扩增子图表解读-理解文章思路(8种图型讲解)
扩增子分析流程-把握分析细节
扩增子统计绘图-冲击高分文章(R绘制8种图型)
微生物组入门必读+宏基因组实操课程
宏基因组分析教程
你想要的生信知识全在这—生信宝典
肠道菌群在人体中的作用
看完此片我想把身上的细菌寄生虫供起来

CanvasXpress包简介

感谢“宏基因组”群友李海敏推荐CanvasXpress包。

CanvasXpress核心是一个JavaScript语言编写的库，主要用于可重复研究。

在R中可以安装canvasXpress包，使用canvasXpress命令绘制各种各样的交互式图形，非常高效方便，实现如下数百种图形，令你的结果瞬间B格爆满。

更多介绍和实例，访问官网：https://canvasxpress.org/html/index.html

源代码、安装及使用示例详见 Github: https://github.com/neuhausi/canvasXpress

canvasXpress安装及加载

安装包位于github上，需要devtools工具安装

# 没有devtools也需要安装
install.packages(devtools)
library(devtools)# 安装canvasXpress
devtools::install_github('neuhausi/canvasXpress')
library(canvasXpress)

绘图实战

1. 三维散点图

# 读数据文件
y <- read.table("http://www.canvasxpress.org/data/cX-irist-dat.txt", header=TRUE, sep="\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 读取分组信息
z <- read.table("http://www.canvasxpress.org/data/cX-irist-var.txt", header=TRUE, sep= "\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 绘制三维散点图，主要参数为数据、分组、分组列、置信椭圆列、图表类型以及相关标签
canvasXpress(data      = y,varAnnot  = z,colorBy   = "Species",ellipseBy = "Species",graphType = "Scatter3D",title     = "Iris Data Set",xAxis     = list("Sepal.Length"),yAxis     = list("Petal.Width"),zAxis     = list("Petal.Length"))

看到按分组绘制了三维散点图，而且添加了分组颜色和椭球形置信区间，我们可以鼠标托动变换角度，也可点选查看点的坐标

2. 矩阵散点图

# 读取数据表
y <- read.table("http://www.canvasxpress.org/data/cX-irist-dat.txt", header=TRUE, sep="\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 分组信息
z <- read.table("http://www.canvasxpress.org/data/cX-irist-var.txt", header=TRUE, sep= "\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 绘制矩阵散点图
canvasXpress(data              = y,varAnnot          = z,graphType         = "Scatter2D",scatterPlotMatrix = TRUE,colorBy           = "Species",showTransition    = TRUE)

二维矩阵散点图展示多条件或因子数据的相关性，点选显示坐标数据

3. 箱线图

# 数据文件
y <- read.table("http://www.canvasxpress.org/data/cX-iris-dat.txt", header=TRUE, sep="\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)# 分组信息
x <- read.table("http://www.canvasxpress.org/data/cX-iris-smp.txt", header=TRUE, sep= "\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)# 绘制箱线图，按物种分大类，再按变量类型分小类
canvasXpress(data              = y,smpAnnot          = x,graphType         = "Boxplot",graphOrientation  = "vertical",title             = "Iris flower data set",smpTitle          = "Species",smpLabelFontStyle = "italic",smpLabelRotate    = 90,xAxis2Show        = FALSE,afterRender       = list(list("groupSamples", list("Species"))))

箱线图比较组间整体数据分布，组内不同因子间比较

4. 多维热图

# 热图数据文件
y  <- read.table("http://www.canvasxpress.org/data/cX-multidimensionalheatmap-dat.txt", header=TRUE, sep="\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 另一个数据矩阵文件
y2 <- read.table("http://www.canvasxpress.org/data/cX-multidimensionalheatmap-dat2.txt", header=TRUE, sep="\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 第三个数据属性矩阵，非连续型
y3 <- read.table("http://www.canvasxpress.org/data/cX-multidimensionalheatmap-dat3.txt", header=TRUE, sep="\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 第四个数据属性矩阵，非连续型
y4 <- read.table("http://www.canvasxpress.org/data/cX-multidimensionalheatmap-dat4.txt", header=TRUE, sep="\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 样品列表
x  <- read.table("http://www.canvasxpress.org/data/cX-multidimensionalheatmap-smp.txt", header=TRUE, sep= "\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 变量列表
z  <- read.table("http://www.canvasxpress.org/data/cX-multidimensionalheatmap-var.txt", header=TRUE, sep= "\t", quote="", row.names=1, fill=TRUE, check.names=FALSE, stringsAsFactors=FALSE)
# 四维热图绘制：y为热图主体，y2为边框颜色，y3为形状，y4为大小
canvasXpress(data          = list(y = y, data2 = y2, data3 = y3, data4 = y4),smpAnnot      = x,varAnnot      = z,graphType     = "Heatmap",guides        = TRUE,outlineBy     = "Outline",outlineByData = "data2",shapeBy       = "Shape",shapeByData   = "data3",sizeBy        = "Size",sizeByData    = "data4")

四个矩阵绘制的热图，分别代表图中颜色、边框颜色、形状和大小，鼠标悬停显示属性值

5. 韦恩图

canvasXpress(vennData   = data.frame(AC=456, A=340, ABC=552, ABCD=148, BC=915, ACD=298, BCD=613, B=562, CD=143, ABD=578, C=620, D=592, AB=639, BD=354, AD=257),graphType  = "Venn",vennLegend = list(A="List 1", D="List 4", C="List 3", B="List 2"),vennGroups = 4)

四组比较韦恩图，需要已知15种情况下的信息。鼠标划过显示信息，双击还可显示固定信息窗口，窗口也可以托动

交互图表导出和交互

以第一个三维散图点图为例：

导出网页格式图片

导入网页：选择Export中的Save as web page，弹窗中输入或选择文件名，如果询问是否覆盖，点yes即可；结束后会显示是否网页显示，选择Try again即网页中打开显示(windows中可能没有)。

我们发现右上角存在一组工具箱，它存在非常多的交互功能。

右上角工具箱有4个按扭，分别为保存图片、移动、探索数据和全屏。其中探索数据中在复杂的操作面板，可以按组、样品选择、查找等操作，功能强大，用户可以自行尝试挖掘更多功能。

Shiny小程序

更强大的交互实例，使用cxShinyExample查看

# 展示所有Shiny示例
cxShinyExample()# 运行其中一个
cxShinyExample(example = "example1")

把数据编写成shiny小程序，可以网页中进入坐标轴选择组合展示数据，图中可用滚轮缩放图片区域。

更多学习资源和最新的软件教程，详见其官网：http://www.canvasxpress.org

学习16S扩增子、宏基因组科研思路和分析实战，关注“宏基因组”

点阅读原文，跳转最新文章目录阅读 https://mp.weixin.qq.com/s/5jQspEvH5_4Xmart22gjMA