python计算相关系数R

方法一:

python 复制代码
import numpy as np

# 计算相关系数R
def r(y_true, y_pred):
    y_true = np.array(y_true)
    y_pred = np.array(y_pred)
    corr = np.corrcoef(y_true, y_pred)[0][1]
    return corr

corr = r(yture, ypred)

方法二

python 复制代码
import scipy.stats
# 计算皮尔逊相关指数,并进行显著性检验
def get_pex(x,y):
    aaa = scipy.stats.pearsonr(x, y)
    print('皮尔逊系数是:{},显著性检测值是:{}'.format(aaa[0],aaa[1]))

pex_r = get_pex(x, y)[0]
pex_p = get_pex(x, y)[1]

参考:https://mp.weixin.qq.com/s/nABJf4WBR4Y4oEZiYKPK4Q

相关推荐
FRing031 分钟前
vue2干货总结7(自定义指令、插槽、单页应用程序、路由的使用(VueRouter))
开发语言·前端·vue.js·vue2
橙橙笔记2 分钟前
具身智能的学习 - 完成一个机械臂 第1节
python·学习·机械臂·具身智能
雪落漂泊5 分钟前
C++算法 二分查找篇(上)
开发语言·c++·算法
小小张说故事10 分钟前
Python 的 class 语句到底执行了什么?拆开 __prepare__、metaclass.__new__ 和那个被讲反的 __set_name__ 顺序
python
波力海苔夹心脆6751 小时前
C# WPF 资源字典详解:ResourceDictionary、MergedDictionaries 与自定义控件案例
开发语言·c#·wpf·入门教程·resource·资源字典
不会算法的笨小孩2 小时前
异常的使用以及会出现的问题
开发语言·c++·c++11·异常
天天被压力2 小时前
【A股分钟级高频数据完全手册 #07】北交所历史分钟K线接口:京市Pro版分时批量
java·人工智能·python
天天被压力2 小时前
【A股分钟级高频数据完全手册 #08】股票日、分时数据接口对比:普通分时与Pro版差异实测
java·人工智能·python
特立独行的猫A3 小时前
AI 时代,我为什么推荐用 cj-tauri 开发桌面应用
java·后端·python
小lo想吃棒棒糖3 小时前
把机器人训练从 GPU 搬到昇腾 NPU:启元 Q1 起身动作 PPO 训练全记录
python·机器人·昇腾·启元机器人