R语言两组变量特征相关关系热图绘制画法

 更新时间:2022年02月18日 09:19:56   作者:hxj7  
本文为大家介绍了如何画两组变量(特征)的相关关系热图的方法详解,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步

准备数据

两组变量的数据可以像下面这样处理,分别保存在两个csv文件中。

> # 导入数据及数据预处理
> setwd("D:/weixin/")
> rows <- read.csv("rows.csv")
> cols <- read.csv("cols.csv")
> str(rows)
'data.frame':   100 obs. of  6 variables:
 $ r1: num  476 482 640 452 308 ...
 $ r2: num  2059 1987 1952 1927 1854 ...
 $ r3: num  513 601 682 497 463 ...
 $ r4: num  2235 2114 2038 1945 1916 ...
 $ r5: num  433 376 525 395 238 ...
 $ r6: num  2028 1943 1802 1775 1748 ...
> str(cols)
'data.frame':   100 obs. of  5 variables:
 $ c1: num  2387 2437 2484 2349 2198 ...
 $ c2: num  540 535 706 509 359 ...
 $ c3: num  472 610 465 473 471 ...
 $ c4: num  74.4 57.3 49.5 51.8 47.6 ...
 $ c5: num  995 915 1038 794 652 ...

简单热图

> # 构建相关关系矩阵
> library(psych)
> data.corr <- corr.test(rows, cols, method="pearson", adjust="fdr")
> data.r <- data.corr$r  # 相关系数
> data.p <- data.corr$p  # p值
> 
> # 画热图
> library(pheatmap)
> pheatmap(data.r, clustering_method="average")

在这里插入图片描述

只对列进行聚类

> pheatmap(data.r, clustering_method="average", cluster_rows=F)

在这里插入图片描述

将相关系数显示在图上

> data.r.fmt <- matrix(sprintf("%.2f", data.r), nrow=nrow(data.p))  # 只保留小数点后两位
> pheatmap(data.r, clustering_method="average", cluster_rows=F, display_numbers=data.r.fmt)

在这里插入图片描述

在图上加上显著性标记

> getSig <- function(dc) {
+   sc <- ''
+   if (dc < 0.01) sc <- '***'
+   else if (dc < 0.05) sc <- '**'
+   else if (dc < 0.1) sc <- '*'
+   sc
+ }
> sig.mat <- matrix(sapply(data.p, getSig), nrow=nrow(data.p))
> str(sig.mat)
 chr [1:6, 1:5] "*" "***" "" "***" "***" "***" "***" "" "***" "**" ...
> pheatmap(data.r, clustering_method="average", cluster_rows=F, display_numbers=sig.mat)

在这里插入图片描述

如果想进一步改变图形效果,可以参考pheatmap函数的用法,修改相应的参数。比如:聚类方式改为complete,加上标题等。

> pheatmap(data.r, clustering_method="complete", cluster_rows=F, display_numbers=sig.mat, main="Corr Heatmap")

在这里插入图片描述

以上就是R语言两组变量特征相关关系热图绘制画法的详细内容,更多关于R语言绘制相关关系热图的资料请关注脚本之家其它相关文章!

相关文章

  • 利用R语言合并数据框的行与列实例代码

    利用R语言合并数据框的行与列实例代码

    实际操作中我们经常需要引入其他表中的列,即将其他表中列加入到表中,需要把两个或者更多的表合并成一个,下面这篇文章主要给大家介绍了关于利用R语言合并数据框的行与列的相关资料,需要的朋友可以参考下
    2022-07-07
  • R语言常用命令集合

    R语言常用命令集合

    这篇文章主要介绍了R语言常用命令集合,文中整理的都是常用的一些命令,有感兴趣的可以多学习下
    2021-03-03
  • R语言-有负下标里才能有零介绍

    R语言-有负下标里才能有零介绍

    这篇文章主要介绍了R语言-有负下标里才能有零的实例讲解,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2021-04-04
  • R语言绘图公式与变量对象混合拼接实现方法

    R语言绘图公式与变量对象混合拼接实现方法

    这篇文章主要为大家介绍了R语言绘图中的公式如何与变量对象混合拼接的实现方法,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步
    2021-11-11
  • R语言boxplot函数深入讲解

    R语言boxplot函数深入讲解

    这篇文章主要介绍了R语言boxplot函数深入讲解,文中图文讲解的很透彻,有感兴趣的同学可以研究下
    2021-03-03
  • R语言 小数点位数的设置方式

    R语言 小数点位数的设置方式

    这篇文章主要介绍了R语言 小数点位数的设置方式,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2021-03-03
  • R语言利用ggplot2绘制QQ图和箱线图详解

    R语言利用ggplot2绘制QQ图和箱线图详解

    这篇文章主要为大家介绍了R语言如何利用ggplot2绘制QQ图和箱线图,文中的示例代码讲解详细,对我们学习R语言有一定的帮助,需要的可以参考一下
    2022-06-06
  • R语言版本升级完全攻略基于Ubuntu

    R语言版本升级完全攻略基于Ubuntu

    这篇文章主要为大家介绍了在Ubuntu上针对不同的版本对R语言版本升级的完全攻略,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步
    2021-11-11
  • R语言-实现按日期分组求皮尔森相关系数矩阵

    R语言-实现按日期分组求皮尔森相关系数矩阵

    这篇文章主要介绍了R语言-实现按日期分组求皮尔森相关系数矩阵,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2021-04-04
  • R语言中c()函数与paste()函数的区别说明

    R语言中c()函数与paste()函数的区别说明

    这篇文章主要介绍了R语言中c()函数与paste()函数的区别说明,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
    2021-04-04

最新评论