科研绘图——upset展示更复杂的交集情况

最新推荐文章于 2024-04-22 22:57:22 发布

生信小鹏

最新推荐文章于 2024-04-22 22:57:22 发布

阅读量1.1k

点赞数 3

分类专栏：科研绘图文章标签： python r语言数据分析数据挖掘

本文链接：https://blog.csdn.net/lijianpeng0302/article/details/128877527

版权

科研绘图专栏收录该内容

4 篇文章 2 订阅

订阅专栏

将几个数据集之间的交集进行可视化，除了可以利用Venn图外，还可以使用upset这种图形。这种方法是发表在2014年nature methods上面。Sets and intersections | Nature Methods

upset图形大致是这个样子

图形基本含义

那怎么观察这种图形，或者图形的意义是什么？
例如下图对这一部分进行了展示

这是UpSet图和Venn图之间的联系

这种图的好处在于，我们可以在矩阵旁边以柱状图的形式绘制交叉点的大小，正如下面的例子中看到

这样可以对交叉点的数量进行直观比较，也可以存储更多的信息。

并且，也可以通过交叉点的数量，来对柱状图进行排序

然后将图形按照水平或者竖直方式进行排列即可。大部分文章当中都是用的水平排列

怎么画一个简单的UpSet图形

UpSet图形可以使用R、python、Javascript都可以绘制，这里使用R绘制

载入数据

加载R包

install.packages("UpSetR")
library(UpSetR)

使用的是UpSetR中的数据，

mutations <- read.csv( system.file("extdata", "mutations.csv", package = "UpSetR"), header=T, sep = ",")

该突变数据集是来自TCGA，代表了多形性胶质母细胞瘤队列中100个突变最多的基因的突变。

行是样本名称，从第2列开始，是基因突变情况

现在就要表现选定的几个基因，有多少突变样本，并且之间的交集如何（也就是某个样本中，几个基因共同突变的情况）。

绘图


upset(mutations, sets = c("PTEN", "TP53", "EGFR", "PIK3R1", "RB1"), sets.bar.color = "#56B4E9",
      order.by = "freq", empty.intersections = "on")

绘制过程还是很简单

其中的解释，就是第一部分的内容。

当然，可以自行修改，让颜色变得艳丽起来

upset(mutations, sets = c("PTEN", "TP53", "EGFR", "PIK3R1", "RB1"),
      keep.order = F,
      sets.bar.color = "red",matrix.color = "blue", shade.color = "gray88",  main.bar.color = "brown",
      color.pal = 3,
      scale.intersections = "identity"
      )