R语言绘制花瓣图flower plot

R语言中有很多现成的R包，可以绘制venn图，但是最多支持5组，当组别数大于5时，venn图即使能够画出来，看上去也非常复杂，不够直观；

在实际的数据分析中，组别大于5的情况还是经常遇到的，这是就可以考虑用花瓣图来进行数据的可视化

比如下面这个例子：

来源于该链接 https://www.researchgate.net/figure/235681265_fig3_The-pan-genome-of-Sinorhizobium-The-flower-plots-and-Venn-diagrams-illustrate-the-number

A和D是我们常见的venn图，B和C 就是花瓣图了

在花瓣图中，我们能够看到两种信息；

1）所有样本共有的信息；

2）每个样本独有的信息；

花瓣图既美观，展示信息也很直观，那么这样的图如何画呢？

我在网上找了半天，也没找到现成的工具，只能自己写代码来画!

在写代码之前，首先来分析下这张图，每一片花瓣就是一个椭圆型，整幅图片可以看做有1个椭圆通过旋转得到

通过以上分析，我们只需要先画一个椭圆，然后循环旋转即可

中间的调试过程就不细讲了，直接看写好的代码

flower_plot <- function(sample, value, start, a, b, ellipse_col = rgb(135, 206, 235, 150, max = 255), circle_col = rgb(0, 162, 214, max = 255),circle_text_cex = 1.5) {par( bty = "n", ann = F, xaxt = "n", yaxt = "n", mar = c(1,1,1,1))plot(c(0,10),c(0,10),type="n")n   <- length(sample)deg <- 360 / nres <- lapply(1:n, function(t){draw.ellipse(x = 5 + cos((start + deg * (t - 1)) * pi / 180), y = 5 + sin((start + deg * (t - 1)) * pi / 180), col = ellipse_col,border = ellipse_col,a = a, b = b, angle = deg * (t - 1))text(x = 5 + 2.5 * cos((start + deg * (t - 1)) * pi / 180),y = 5 + 2.5 * sin((start + deg * (t - 1)) * pi / 180),value[t])if (deg * (t - 1) < 180 && deg * (t - 1) > 0 ) {text(x = 5 + 3.3 * cos((start + deg * (t - 1)) * pi / 180),y = 5 + 3.3 * sin((start + deg * (t - 1)) * pi / 180),sample[t],srt = deg * (t - 1) - start,adj = 1,cex = circle_text_cex)} else {text(x = 5 + 3.3 * cos((start + deg * (t - 1)) * pi / 180),y = 5 + 3.3 * sin((start + deg * (t - 1)) * pi / 180),sample[t],srt = deg * (t - 1) + start,adj = 0,cex = circle_text_cex)}           })draw.circle(x = 5, y = 5, r = 1.3, col = circle_col, border = circle_col)
}

这里我写成了1个函数，函数的调用方式如下;

flower_plot(c("WSM419", "A321", "M1", "M2", "M22", "M58", "M102", "M161", "KH36b", "KH36c", "KH36d", "KH53a", "KH53b"),c(519, 556, 83, 62, 415, 425, 357, 441, 22, 41, 33, 44, 43), 90, 0.5, 2)

第一个参数为样本名字构成的向量，第二个参数为每个样本独有的数目，第三个参数为起始椭圆的角度，第四个参数为椭圆的短轴的长度，第五个参数为椭圆的长轴的长度

效果图如下：

R语言绘制花瓣图flower plot相关推荐

R语言绘制火山图（volcano plot）实战：为差异表达基因（DEGs）添加颜色、基于显著性阈值进行点的颜色美化、为选定基因添加标签
R语言绘制火山图(volcano plot)实战:为差异表达基因(DEGs)添加颜色.基于显著性阈值进行点的颜色美化.为选定基因添加标签目录 R语言绘制火山图(volcano plot)实战 #导入 ...
R语言绘制棒棒糖图（火柴杆图）
本博客介绍几种利用R语言绘制棒棒糖图(火柴杆图)的方法. 2. 使用原生ggplot方法最容易也是最简单想到的方法是直接使用ggplot2包进行更新,这里需要使用ggplot本身的特性,通过图层叠加 ...
r语言绘制雷达图_用r绘制雷达蜘蛛图
r语言绘制雷达图 I've tried several different types of NBA analytical articles within my readership who are ...
R语言绘制线图（line）实战
R语言绘制线图(line)实战目录 R语言绘制线图(line)实战 #仿真数据 #基础线图
R语言绘制空白图实战
R语言绘制空白图实战目录 R语言绘制空白图实战 #绘制空白图1 #绘制空白图2 #绘制空白图3
运用R语言绘制小提琴图
运用R语言绘制小提琴图一.概念小提琴图是一种绘制连续型数据的方法,可以认为是箱形图与核密度图的结合体,与此同时,还可使用核密度图展示数据分布的'轮廓'效果,'轮廓'越大,即意味着数据越集中于该处, ...
matlab 画qq图,科学网—[转载]R语言绘制QQ图 - 刘朋的博文
R语言绘制QQ图实例1: #############加载数据 data R R=apply(R,2,as.numeric) #R语言将字符串矩阵转化为数值型矩阵,apply()函数里面的第2个值,如 ...
R语言绘制热图（其实是相关系数图）实践(二）corrplot包
目录前言 corrplot包简介语法和常用参数介绍函数语法参数介绍实践 summary 参考资料前言在我的上一篇的内容中(R语言绘制热图实践(一)pheatmap包 ),我以绘制相关系数 ...
R语言绘制QQ图实战（qqplot函数、qqnorm函数、qqline函数）
R语言绘制QQ图实战(qqplot函数.qqnorm函数.qqline函数) 目录 R语言绘制QQ图实战(qqplot函数.qqnorm函数.qqline函数)
R语言绘制小提琴图geom_violin，如何给每个小提琴添加自定义直线？
R语言绘制小提琴图geom_violin,如何给每个小提琴添加自定义直线? 效果如下图所示: 数据按照下面的方式进行存储下面提供所需要的代码: #加载必要的包 rm(list = ls()) lib ...

R语言绘制花瓣图flower plot

R语言绘制花瓣图flower plot相关推荐

最新文章

热门文章