当前位置:网站首页>皮尔森相关系数
皮尔森相关系数
2022-07-05 08:42:00 【python-码博士】
# encoding:utf-8
import pandas as pd
from pandas import DataFrame
import matplotlib.pyplot as plot
import math
#target_url = ("https://archive.ics.uci.edu/ml/machine-learning-databases/undocumented/connectionist-bench/sonar/sonar.all-data")
#rockVMines = pd.read_csv(target_url ,header=None,prefix="V") #prefix前缀
rockVMines = pd.read_csv('../rockdata.txt',header=None,prefix="V") #prefix前缀
row2 = rockVMines.iloc[1,0:60]
row3 = rockVMines.iloc[2,0:60]
n = len(row2)
mean2 = row2.mean()
mean3 = row3.mean()
t2=0 ; t3=0;t1=0
for i in range(n):
t2 += (row2[i] - mean2) * (row2[i] - mean2) / n
t3 += (row3[i] - mean3) * (row3[i] - mean3) / n
r23=0
for i in range(n):
r23 += (row2[i] - mean2)*(row3[i] - mean3)/(n* math.sqrt(t2 * t3))
print(r23)
corMat = DataFrame(rockVMines.corr()) #corr 求相关系数矩阵
print(corMat)
plot.pcolor(corMat)
plot.show()
边栏推荐
猜你喜欢
Shift operation of complement
C# LINQ源码分析之Count
Yolov4 target detection backbone
Low code platform | apaas platform construction analysis
STM32 summary (HAL Library) - DHT11 temperature sensor (intelligent safety assisted driving system)
实例010:给人看的时间
剑指 Offer 09. 用两个栈实现队列
猜谜语啦(4)
TypeScript手把手教程,简单易懂
319. 灯泡开关
随机推荐
实例008:九九乘法表
Business modeling | process of software model
Chapter 18 using work queue manager (1)
剑指 Offer 06. 从尾到头打印链表
319. 灯泡开关
U8g2 drawing
Typical low code apaas manufacturer cases
Search data in geo database
C# LINQ源码分析之Count
Business modeling of software model | object modeling
[daily training] 1200 Minimum absolute difference
Typescript hands-on tutorial, easy to understand
Guess riddles (6)
Example 004: for the day of the day, enter a day of a month of a year to judge the day of the year?
MATLAB小技巧(28)模糊綜合評價
【三层架构】
[daily training -- Tencent selected 50] 557 Reverse word III in string
The first week of summer vacation
Latex improve
2022.7.4-----leetcode.1200