​最近有小伙伴问到带mantel test检验的相关性矩阵怎么画,于是小编就在网上找到了相关的绘图示例,整理了网上各位大神关于ggcor包的下载以及该包自带的示例代码供各位小伙伴参考使用~

首先是要安装ggcor这个包。由于该包的作者已经下架了这个包,所以只能自己下载压缩包本地安装,小编主要参考的是这篇文章:
https://www.zhihu.com/question/409189621/answer/1842210647

示例代码
ggcor包自带的几个绘图示例

Correlation plot

library(ggplot2)
library(ggcor)
set_scale()
quickcor(mtcars) + geom_square()
quickcor(mtcars, type = "upper") + geom_circle2()
quickcor(mtcars, cor.test = TRUE) +geom_square(data = get_data(type = "lower", show.diag = FALSE)) +geom_mark(data = get_data(type = "upper", show.diag = FALSE), size = 2.5) +geom_abline(slope = -1, intercept = 12)


Mantel test plot

library(dplyr)
library(ggplot2)
library(ggcor)data("varechem", package = "vegan")
data("varespec", package = "vegan")mantel <- mantel_test(varespec, varechem,spec.select = list(Spec01 = 1:7,Spec02 = 8:18,Spec03 = 19:37,Spec04 = 38:44)) %>% mutate(rd = cut(r, breaks = c(-Inf, 0.2, 0.4, Inf),labels = c("< 0.2", "0.2 - 0.4", ">= 0.4")),pd = cut(p.value, breaks = c(-Inf, 0.01, 0.05, Inf),labels = c("< 0.01", "0.01 - 0.05", ">= 0.05")))quickcor(varechem, type = "upper") +geom_square() +anno_link(aes(colour = pd, size = rd), data = mantel) +scale_size_manual(values = c(0.5, 1, 2)) +scale_colour_manual(values = c("#D95F02", "#1B9E77", "#A2A2A288")) +guides(size = guide_legend(title = "Mantel's r",override.aes = list(colour = "grey35"), order = 2),colour = guide_legend(title = "Mantel's p", override.aes = list(size = 3), order = 1),fill = guide_colorbar(title = "Pearson's r", order = 3))


Circular heatmap

rand_correlate(100, 8) %>% ## require ambient packagesquickcor(circular = TRUE, cluster = TRUE, open = 45) +geom_colour(colour = "white", size = 0.125) +anno_row_tree() +anno_col_tree() +set_p_xaxis() +set_p_yaxis()


General heatmap

d1 <- rand_dataset(20, 30) %>% gcor_tbl(cluster = TRUE)
p <- matrix(sample(LETTERS[1:4], 90, replace = TRUE), nrow = 30,dimnames = list(paste0("sample", 1:30), paste0("Type", 1:3))) %>% gcor_tbl(name = "Type", row.order = d1) %>% qheatmap(aes(fill = Type)) + coord_fixed() + remove_y_axis()
d2 <- data.frame(x = sample(paste0("var", 1:20), 200, replace = TRUE))set_scale()
quickcor(d1) +geom_colour(aes(fill = value)) +anno_hc_bar(width = 1) +anno_row_custom(p) +anno_row_tree() +anno_hc_bar(pos = "top") +anno_bar(d2, aes(x = x), height = 0.12) +anno_col_tree(height = 0.12)


“作图帮”微信公众号同步更新~

图图即将上线一个永久免费的云绘图工具,小伙伴可关注“作图帮”公众号后扫码进入“生信作图交流群”,群内届时会提供绘图账号供各位小伙伴使用~

ggcor包的安装与绘图示例相关推荐

  1. R语言安装ggcor包

    找了很多方法,一直安装不上去,后来把包下载了本地,用本地安装的方法,也还是报错.后来兜兜转转找到了这个方法,非常好用!参考了看见你啦,自己根据实际情况,也稍作修改. 1.下载ggcor包到本地 ggc ...

  2. linux泰语语言包,linux安装中文语言包(示例代码)

    相关配置如下: yum install fonts-chinese.noarch yum install m17n-db-common-cjk yum install m17n-db-chinese ...

  3. linux gnuplot yum源,在CentOS上编译安装gnuplot绘图工具

    gnuplot是一个命令行的交互式绘图工具.用户通过输入命令,可以逐步设置或修改绘图环境,并以图形描述数据或函数,使我们可以借由图形做更进一步的分析. gnuplot是由Colin Kelly和Tho ...

  4. R语言-包的安装、载入及使用方法

    一.原理简述 包是R函数.数据.预编译代码以一种定义完善的格式组成的集合.计算机上存储包的目录称为库(library).函数.libPaths()能够显示库所在的位置,函数library()则可以显示 ...

  5. matplotlib 和 pandas 两个包的安装

    matplotlib 和 pandas 两个包的安装 matplotlib是强大的python 绘图包.pandas 是强大的python分析工具包.numpy是强大的python统计包. 都超级好用 ...

  6. linux平台下rpm方式和源码包方式安装mysql5.7

    博主QQ:819594300 博客地址:http://zpf666.blog.51cto.com/ 有什么疑问的朋友可以联系博主,博主会帮你们解答,谢谢支持! 一.下载mysql的rpm包 Mysql ...

  7. Rpm包的安装与yum的配置

    RPM包管理器的常用命令和示例: 1.安装程序包: rpm [option]/path/to/package_file –i:安装程序包 –v:显示安装的过程 -h :显示安装进度 -vv:显示详细信 ...

  8. debian 安装_如何在 Debian/Ubuntu 系统中检查程序包是否安装?

    有很多种方法可以检查程序包是否已安装,我找到了 7 种命令,你可以从中选择你喜欢的使用. -- Magesh Maruthamuthu(作者) 我们近期发布了一篇关于批量程序包安装的文章.在此同时,关 ...

  9. 文本分析常用R包的安装(Rweibo、wordcloud、tm、tmcn、Rwordseg、Rcharts、xlsx、XLConnect)

    为了使用文本分析挖掘,必须要先安装一些R包,比如Rwordseg包常用于分词.tmcn用于词频统计.wordcloud用于绘制文字云.在R中,有些包是需要额外设置安装的,且设置起挺麻烦的,因此,以下总 ...

最新文章

  1. Nginx vs Apache--reference
  2. ADO.NET Entity Framework Extensions 简单应用
  3. Vs2010删除空白行
  4. SegmentFault 创始人祁宁对话 C# 之父 Anders Hejlsberg
  5. .net core 基于Dapper 的分库分表开源框架(core-data)
  6. springboot启动报错LoggerFactory is not a Logback LoggerContext but Logback...
  7. Windows Linux Mac 路由添加删除
  8. POJ读书笔记2.1 —— 鸡兔笼带
  9. Java日志系统学习之log4j!
  10. 进销存软件网络版连接不上,如何设置IP地址
  11. 保姆级的接口自动化教程,不会写代码也能2小时学会
  12. LIBOR-OIS息差利率 伦敦银行同业拆息与隔夜指数掉期
  13. matlab对图片边缘化处理
  14. java nurbs几何库_OpenGL超级宝典笔记——NURBS与曲面细分
  15. 打通零售新渠道,实现线上线下业务一体化
  16. umi首屏加载速度优化
  17. powershell创建新员工AD账号
  18. WPF 设置属性“System.Windows.FrameworkElement.Style”时引发了异常
  19. 利用Composer搭建企业内部仓库
  20. 特斯拉新款Model Y曝光 续航仅449公里

热门文章

  1. Redis 日志篇:无畏宕机实现高可用的杀手锏
  2. 选择大于努力?某程序员吐槽:自己毕业去了百度,技术不如自己的室友去了快手,如今股票1500万!...
  3. 年轻人不讲武德,where 1=1 是什么鬼?
  4. 基于Spring Boot+Cloud构建微云架构
  5. MQTT:物联网的神经系统
  6. MongoDb优化指南
  7. 强烈推荐10个新媒体运营必备工具,极大提高工作效率
  8. 如何打造企业的组织能力?战略决定组织,而组织决定成败!
  9. 基于linux的驱动设计,《基于LINUX的虚拟驱动设计》-毕业论文.doc
  10. QT OpenCV Linux