给出皮尔逊相关系数的python代码

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