matlab 多个热图,MATLAB实例：求相关系数、绘制热图并找到强相关对

MATLAB实例：求相关系数、绘制热图并找到强相关对

用MATLAB编程，求给定数据不同维度之间的相关系数，并绘制热图，保存图片，找到强相关的维度对。

1. demo.m

clear

clc

etea=0.8; %阈值

% 加载数据

data_load=dlmread('E:\scanplot\wine.data');

data=data_load(:,2:14);

[N, D]=size(data);

% 求维度之间的相关系数

rho = corr(data, 'type','pearson');

% 绘制热图

string_name={'Alcohol','Malic acid','Ash','Alcalinity of ash','Magnesium','Total phenols','Flavanoids','Nonflavanoid phenols','Proanthocyanins','Color intensity','Hue','Diluted wines','Proline'};

xvalues = string_name;

yvalues = string_name;

h = heatmap(xvalues,yvalues,rho, 'FontSize',10, 'FontName','Times New Roman');

h.Title = 'Correlation Coefficient';

colormap(jet)

saveas(gcf,sprintf('wine相关系数热图.jpg'),'bmp'); %保存图片

% 绝对值

rho=abs(rho);

rho_1=rho.*tril(ones(D,D),-1); %下三角

[row, col]=find(rho_1>etea); %找>etea的两个维度

[Num, ~]=size(row);

% A:存放相关系数>etea的两个维度及相关系数值

A=zeros(Num, 3);

for i=1:Num

A(i,:)=[row(i), col(i), rho_1(row(i), col(i))];

fprintf('强线性相关的两个维度是: 第%d个维度: %s与第%d个维度: %s, 两者的相关系数为：%f\n', row(i), string_name{row(i)}, col(i), string_name{col(i)}, rho_1(row(i), col(i)));

end

2. 结果

>> demo

强线性相关的两个维度是: 第7个维度: Flavanoids与第6个维度: Total phenols, 两者的相关系数为：0.864564

MATLAB查看并设置当前颜色图：colormap

pheatmap绘制“热图”，你需要的都在这

热图可以聚合大量的数据,并可以用一种渐进色来优雅地表现,可以很直观地展现数据的疏密程度或频率高低. 本文利用R语言 pheatmap 包从头开始绘制各种漂亮的热图.参数像积木,拼凑出你最喜欢的热图即可 ...

MATLAB实例：散点密度图

MATLAB实例:散点密度图 作者:凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ MATLAB绘制用颜色表示数据密度的散点图 数据来源:MATLAB中“fit ...

MATLAB实例：聚类网络连接图

MATLAB实例:聚类网络连接图 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 本文给出一个简单实例,先生成2维高斯数据,得到数据之后,用模糊C均值( ...

基于matplotlib的数据可视化 - 热图imshow

热图: Display an image on the axes. 可以用来比较两个矩阵的相似程度 mp.imshow(z, cmap=颜色映射,origin=垂直轴向) imshow( X, cma ...

扩增子图表解读3热图：差异菌、OTU及功能

热图是使用颜色来展示数值矩阵的图形.通常还会结合行.列的聚类分析,以表达实验数据多方面的结果.  热图在生物学领域应用广泛,尤其在高通量测序的结果展示中很流行,如样品-基因表达,样品-OTU相对丰度矩 ...

HTML5 WebAudioAPI(三)--绘制频谱图

HTML

用R包中heatmap画热图

一:导入R包及需要画热图的数据 library(pheatmap) data

python matplotlib imshow热图坐标替换/映射

今天遇到了这样一个问题,使用matplotlib绘制热图数组中横纵坐标自然是图片的像素排列顺序, 但是这样带来的问题就是画出来的x,y轴中坐标点的数据任然是x,y在数组中的下标, 实际中我们可能期望坐 ...

R语言学习 - 热图简化

绘制热图除了使用ggplot2,还可以有其它的包或函数,比如pheatmap::pheatmap (pheatmap包中的pheatmap函数).gplots::heatmap.2等.   相比于gg ...

随机推荐

编译libjpeg库

最近在写车牌识别软件,需要用到BMP转成JPG的功能,自然就想到借助libjpeg来实现 OS: win7 64位 编译器: VS2008 1. 下载源代码下载地址:http://www.ijg.or ...

Python 字典(Dictionary)

字典是另一种可变容器模型,且可存储任意类型对象. 字典的每个键值(key=>value)对用冒号(:)分割,每个对之间用逗号(,)分割,整个字典包括在花括号({})中 ,格式如下所示: d = ...

LeetCode——Valid Palindrome

Given a string, determine if it is a palindrome, considering only alphanumeric characters and ignori ...

Spring注解基本解读

在一个类中使用Spring对象,办法如下: 使用注解的形式注入 从ApplicationContext中获取. T t = new ApplicationContext.getBean("x ...

Ajax理解总结

前端开发拿数据页面实时更新是离不开Ajax这个技术的 AJAX即"Asynchronous Javascript And XML"(异步JavaScript和XML),是指一种创建 ...

cookieUtil

public class CookieUtil { /** * 设置cookie * @param name cookie名字 * @param value cookie值 * @param maxA ...

golang标准库 context的使用

本文索引 问题引入 context包简介 示例 问题引入 goroutine为我们提供了轻量级的并发实现,作为golang最大的亮点之一更是备受推崇. goroutine的简单固然有利于我们的开发,但 ...

hiho一下 第206周

题目1 : Guess Number with Lower or Higher Hints 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 There is a game ...

javaweb项目编译错误

Eclipse Maven 开发一个 jee 项目时,编译时遇到以下错误:Description Resource Path Location TypeDynamic Web Module 3.0 r ...

【基于初学者的SSH】struts2 的拦截器、令牌的简单应用及理解

一:拦截器与过滤器类似,但是它们的区别也很大: 01):过滤器理论上可以过滤任意内容,比如HTML,servlet,jsp,图片路径 02):拦截器只可以拦截action. 二:拦截器的原理  act ...