给出皮尔逊相关系数的python代码
时间: 2023-05-26 07:03:18 浏览: 157
import numpy as np
def pearson(x, y):
"""
计算皮尔逊相关系数
:param x: 第一个变量
:param y: 第二个变量
:return: 皮尔逊相关系数
"""
# 计算均值和标准差
x_mean = np.mean(x)
y_mean = np.mean(y)
x_stdev = np.std(x)
y_stdev = np.std(y)
# 计算协方差
covariance = np.sum((x - x_mean) * (y - y_mean))
# 计算皮尔逊相关系数
correlation = covariance / (len(x) * x_stdev * y_stdev)
return correlation
相关问题
皮尔逊相关系数python代码
当然,我可以为您提供皮尔逊相关系数的Python代码:
```python
import numpy as np
def pearson_correlation(x, y):
# 计算平均值
x_mean = np.mean(x)
y_mean = np.mean(y)
# 计算分子
numerator = np.sum((x - x_mean) * (y - y_mean))
# 计算分母
denominator = np.sqrt(np.sum((x - x_mean) ** 2)) * np.sqrt(np.sum((y - y_mean) ** 2))
# 防止分母为0
if denominator == 0:
return 0
# 计算皮尔逊相关系数
corr = numerator / denominator
return corr
```
此函数计算两个向量x和y的皮尔逊相关系数。 它使用numpy库来进行向量操作。
读取数据计算皮尔逊相关系数python代码
当然,我可以为您提供读取数据计算皮尔逊相关系数的Python代码。以下是一个简单的例子:
```
import pandas as pd
import numpy as np
from scipy.stats import pearsonr
#读取数据
data = pd.read_csv('data.csv')
#计算皮尔逊相关系数
corr, _ = pearsonr(data['column1'], data['column2'])
print('皮尔逊相关系数:', corr)
```
注:在此代码中,您需要用实际的列名替换“column1”和“column2”,并且将文件名“data.csv”替换为您实际的数据文件名。
阅读全文