返回
python计算相关系数（python实践统计学中的三大相关性系数）
python计算相关系数（python实践统计学中的三大相关性系数）
2024-09-29 09:31:32  作者:赤道等你爱  网址:https://m.xinb2b.cn/tech/bjn527747.html

本文首发地址：

https://yishuihancheng.blog.csdn.net/article/details/83547648

欢迎关注我的博客【Together_CZ】，我是沂水寒城！


今天我简单地使用了scipy模块进行了统计学中三大相关性分析方法（皮尔森相关性系数、斯皮尔曼相关性系数、肯德尔相关性系数）的使用，这里的主要工作是将相关性计算结果进行可视化展示，以便于更直观地来对不同因子之间的相关性进行分析，下面是具体的实践：

#!usr/bin/env python#encoding:utf-8 '''__Author__:沂水寒城功能： python实践统计学中的三大相关性系数，并绘制相关性分析的热力图''' import numpy as npimport pandas as pdimport seaborn as snsimport matplotlib.pyplot as plt from scipy.stats import pearsonr,spearmanr def randomPlot(): ''' 构造随机数矩阵来绘制热力图 ''' data=np.random.rand(8,8) print(data) fig,ax=plt.subplots(figsize=(10,10)) key_list=['A','B','C','D','E','F','G','H'] sns.heatmap(pd.DataFrame(np.round(data,4),columns=key_list,index=key_list),annot=True,vmax=1,vmin=0,xticklabels=True, yticklabels=True,square=True,cmap="YlGnBu") ax.set_title(' Heat Map ', fontsize=18) ax.set_ylabel('Y', fontsize=18) ax.set_xlabel('X', fontsize=18) plt.savefig('Random.png') def dataPlot(): ''' 基于相关性系数计算结果来绘制 ''' data1=[[1.0, 0.679138216972867, 0.8723347117585479, 0.8678591195124671, 0.9091565423610072, -0.681888271075189, -0.7228630059512009, -0.681888271075189], [0.679138216972867, 1.0, 0.7808137626700863, 0.45749648218056993, 0.5348626411851899, -0.44326149230111495, -0.39558814315535173, -0.44326149230111495], [0.8723347117585479, 0.7808137626700863, 1.0, 0.7020628568851344, 0.7658545326255977, -0.7217424383082518, -0.5875065812285509, -0.7217424383082518], [0.8678591195124671, 0.45749648218056993, 0.7020628568851344, 1.0, 0.9903366506522474, -0.5874500887255807, -0.8382341966002466, -0.5874500887255807], [0.9091565423610072, 0.5348626411851899, 0.7658545326255977, 0.9903366506522474, 1.0, -0.6439505593700681, -0.8244259230387567, -0.6439505593700681], [-0.681888271075189, -0.44326149230111495, -0.7217424383082518, -0.5874500887255807, -0.6439505593700681, 1.0, 0.4821038714558963, 1.0], [-0.7228630059512009, -0.39558814315535173, -0.5875065812285509, -0.8382341966002466, -0.8244259230387567, 0.4821038714558963, 1.0, 0.4821038714558963], [-0.681888271075189, -0.44326149230111495, -0.7217424383082518, -0.5874500887255807, -0.6439505593700681, 1.0, 0.4821038714558963, 1.0]] data2=[[0.9999999999999999, 0.8345554576233073, 0.9243854856819417, 0.9330027648545496, 0.9482209505391218, -0.6482501837911674, -0.8818546914603946, -0.6482501837911674], [0.8345554576233073, 0.9999999999999999, 0.9003585423988755, 0.7649830893396171, 0.7979521001324441, -0.4691190784521079, -0.6480284658216039, -0.4691190784521079], [0.9243854856819417, 0.9003585423988756, 1.0, 0.8506386730833658, 0.8830994621736679, -0.7006298478257242, -0.7676395373967932, -0.7006298478257242], [0.9330027648545496, 0.7649830893396171, 0.8506386730833659, 0.9999999999999999, 0.9951726440650351, -0.5278605167290854, -0.9273089460158745, -0.5278605167290854], [0.9482209505391219, 0.7979521001324442, 0.8830994621736679, 0.9951726440650353, 1.0, -0.5724194583833783, -0.9112198965249181, -0.5724194583833783], [-0.6482501837911674, -0.46911907845210793, -0.7006298478257242, -0.5278605167290855, -0.5724194583833783, 1.0, 0.49515922475047763, 1.0], [-0.8818546914603947, -0.6480284658216039, -0.7676395373967932, -0.9273089460158744, -0.9112198965249181, 0.49515922475047763, 1.0, 0.49515922475047763], [-0.6482501837911674, -0.46911907845210793, -0.7006298478257242, -0.5278605167290855, -0.5724194583833783, 1.0, 0.49515922475047763, 1.0]] data3=[[1.0, 0.6447733742285494, 0.7648449086941359, 0.7781871747188993, 0.8077705036447606, -0.45314582645069557, -0.7063398069618333, -0.45314582645069557], [0.6447733742285494, 1.0, 0.723107089412363, 0.5724058545087037, 0.6071562414937144, -0.3151362329600627, -0.47203186715609546, -0.3151362329600627], [0.7648449086941359, 0.7231070894123629, 1.0, 0.6636168802979515, 0.7036503531839524, -0.5094145396347782, -0.5801992046091554, -0.5094145396347782], [0.7781871747188993, 0.5724058545087037, 0.6636168802979514, 1.0, 0.9417388451103401, -0.35110754201106087, -0.818036667286979, -0.35110754201106087], [0.8077705036447606, 0.6071562414937144, 0.7036503531839524, 0.9417388451103401, 1.0, -0.3871446641295609, -0.7851151107221745, -0.3871446641295609], [-0.45314582645069557, -0.3151362329600627, -0.5094145396347782, -0.35110754201106087, -0.3871446641295609, 1.0, 0.331222719795913, 1.0], [-0.7063398069618333, -0.4720318671560954, -0.5801992046091554, -0.818036667286979, -0.7851151107221744, 0.331222719795913, 1.0, 0.331222719795913], [-0.45314582645069557, -0.3151362329600627, -0.5094145396347782, -0.35110754201106087, -0.3871446641295609, 1.0, 0.331222719795913, 1.0]] data1,data2,data3=np.array(data1),np.array(data2),np.array(data3) fig,ax=plt.subplots(figsize=(10,10)) key_list=['A','B','C','D','E','F','G','H'] sns.heatmap(pd.DataFrame(np.round(data1,4),columns=key_list,index=key_list),annot=True,vmax=1,vmin=0,xticklabels=True, yticklabels=True,square=True,cmap="YlGnBu") ax.set_title(' Heat Map ', fontsize=18) ax.set_ylabel('Y', fontsize=18) ax.set_xlabel('X', fontsize=18) plt.savefig('data1.png') if __name__=='__main__': randomPlot() dataPlot()

其中，随机生成数矩阵得到的热力图如下：


基于相关性计算结果得到的三个热力图分别如下：

data1.png


data2.png


data3.png


从上面三个热力图可以看到：虽然三者的计算结果有些偏差，但是整体的趋势都是相同的。

  • 柘城起台镇舞之韵舞蹈（每晚演出为居民带来欢乐）
  • 2024-09-29每晚演出为居民带来欢乐每当夜幕降临，新落成的陵城区糜镇吉祥社区陆乡公园内灯火璀璨人头攒动，广场上一支广场舞队伍伴随着优美的音乐翩翩起舞，她们个个脸上洋溢着幸福、快乐和自信快乐的旋律、优美的舞姿，形成了社区夜晚一道最靓丽的风。
  • 韩国街拍女人漂亮的穿搭（韩国女孩初秋街拍合集）
  • 2024-09-29韩国女孩初秋街拍合集韩国女生有多会穿衣不用我说，相信很多小姐姐都深有体会了，在穿衣上她们特别擅长搭配，巧用基础单品来搭配出高级感秋天到了，你的战衣准备好了么？为了帮助大家变得更美，这次给大家整理了一些韩国街头穿搭示范，很。
  • 汽车自动启停系统关闭（如何关掉让人闹心的汽车自动启停系统）
  • 2024-09-29如何关掉让人闹心的汽车自动启停系统自动启停系统就是在堵车和等红灯时，发动机可以自动熄火，起步时再自动起动发动机，可以减少燃油消耗，降低废气排放自动启停系统在一定程度上还能规范驾驶，不系好安全带自动启停系统就不能工作，停车后解开安全带，。
  • 薛家燕资料简介（薛家燕资料）
  • 2024-09-29薛家燕资料薛家燕，1950年3月30日出生于中国香港，香港资深女艺人，七十年代的粤语片女星，港星中受人尊敬的女前辈之一1959年，作为童星开始拍电影，处女作是《七儿八女九状元》；1969年，拍摄由龙刚导演的成名。
  • 美国最著名的10大野鸡大学（美国顶尖大学纷纷开放线下访校）
  • 2024-09-29美国顶尖大学纷纷开放线下访校这两年由于疫情影响，给原本热闹的实地探校按下了暂停键，但是随着疫情得到良好的控制，越来越多的美国顶尖大学纷纷开放了线下访校，你是不是迫不及待地想去学校一探究竟了？目前美国TOP30综合性大学提供线下访。
  • 中通快递的技术分析（看好中通快递的四个原因）
  • 2024-09-29看好中通快递的四个原因金融界网5月7日消息，中通快递去年少赚24％，绩后股价跌穿上市价并低见200.5港元，然而时隔一个多月股价已反弹超过2成，以下是可继续睇好公司前景的原因1.内地经济复苏势头依然强劲中国今年首季GDP增。
  • 怎么做营养的披萨（披萨做法）
  • 2024-09-29披萨做法用高筋面粉、水、食盐、白糖、黄油、酵母和成光滑的面团后等待发酵将发酵好的面团排气，擀成盘子大小的圆饼放入相应尺寸的披萨盘中，用手按成边缘厚，中间薄的饼坯，用叉子扎孔，以防烤制的时候膨胀鼓起把蔬菜改刀，。
  • 鲍鱼汁怎么调好吃（鲍鱼汁的调法）
  • 2024-09-29鲍鱼汁的调法先将干鲍鱼用冷水浸泡，大概泡两天左右用冷水泡开之后，再用温水泡一会，直到鲍鱼涨起来为止将鲍鱼清洗干净后，放入水中蒸10个小时左右将火腿、猪瘦肉、鸡爪、猪皮及鲍鱼，注入剩余的二汤，用猛火烧沸后，调入蚝油。
  • 自动挡挂什么挡可以前进（自动挡上）
  • 2024-09-29自动挡上自动挡车上“”“-”档到底有什么用？老司机告诉你，答案很明了人们生活水平的提高，在买车这一块都会选择自动挡车，毕竟自动挡车在操作的时候更加的方便然而对于自动挡车来说档位这块操作和方便，但是在实际过程中。
  • 全新大众速腾是否还有越级表现（新款大众速腾减配增车型）
  • 2024-09-29新款大众速腾减配增车型大众速腾作为“A级”这个概念的开创者，一直以来都称得上是这个级别标杆车型在其进入国内市场至今，累计销量已经突破320万台因此不论是口碑还是产品力，都是这个级别中首屈一指的存在纵观国内A级轿车市场，基本。