np.corrcoef返回相关系数矩阵而非单个数值：传入两个一维数组得2×2矩阵，[0,1]位置为皮尔逊相关系数；需预处理NaN并注意输入形状（按行视为变量）。
corrcoef 返回的不是单个数值，而是相关系数矩阵
很多人调用 np.corrcoef 后看到一个 2×2 数组就懵了——它默认把两个输入当作两组变量，返回的是完整的协方差归一化矩阵。比如传入两个一维数组 a 和 b，结果 np.corrcoef(a, b) 的形状是 (2, 2)，其中 [0, 1] 和 [1, 0] 才是你想要的皮尔逊相关系数（对称），[0, 0] 和 [1, 1] 恒为 1。
- 想直接拿到标量：取
np.corrcoef(a, b)[0, 1]，别写成[0][1]（虽等价但可读性差） - 如果只传一个数组，比如
np.corrcoef(a)，它会把a当作“单变量多观测”，返回 1×1 矩阵[[1.]]，毫无意义 - 传入二维数组时，
np.corrcoef默认按行组织变量（即每行是一组变量，列是样本），和 pandas 默认按列不同，容易搞反
NaN 值会让 corrcoef 直接返回全 NaN 矩阵
np.corrcoef 对缺失值零容忍：只要输入里有 np.nan，整个输出矩阵就全是 nan，不会跳过或警告。这不是 bug，是设计如此——它底层调用的是 np.cov，而后者也不处理 NaN。
- 必须提前清理：用
a = a[~np.isnan(a)]和b = b[~np.isnan(b)]同步剔除对应位置的 NaN（注意不是分别剔除，否则长度不一致） - 更稳妥的做法是构造掩码：
mask = ~(np.isnan(a) | np.isnan(b))，再用a[mask]和b[mask] - 别依赖
np.ma.corrcoef：它行为不一致，且对非 masked array 输入可能报错或静默失败
和 scipy.stats.pearsonr 的关键区别在哪
如果你只需要两个向量间的相关系数和 p 值，scipy.stats.pearsonr 更直接。它返回的是元组 (r, p-value)，不是矩阵，也不强制要求同长度（会自动对齐）。
-
np.corrcoef快、纯 NumPy、支持批量多变量（比如一次算 10 个变量两两相关），但接口笨重 -
pearsonr自动处理 NaN（跳过含 NaN 的配对），还附带统计显著性，但一次只能算一对 - 性能上：对百万级数据，
np.corrcoef比循环调用pearsonr快 10 倍以上；但若只算一两次，用pearsonr更省心
用 corrcoef 算多变量相关矩阵时，行/列顺序决定变量含义
当你传入 shape 为 (n_features, n_samples) 的二维数组（比如 X），np.corrcoef(X) 返回的是 n_features × n_features 矩阵，其中第 i 行第 j 列是第 i 个变量和第 j 个变量的相关系数。
- 常见错误：把原始数据按列存（pandas 默认），却直接传给
np.corrcoef，结果得到的是“每个样本和其他样本的相关性”，完全不是你想要的变量间关系 - 正确做法：传入前转置，
np.corrcoef(X.T)，或者确保输入是 (n_vars, n_obs) 形状 - 验证方法：检查对角线是否全为 1；如果不是，说明维度理解错了
实际用的时候，最常卡住的不是公式，而是数据形状和 NaN 处理——这两个地方出错，结果看起来“算出来了”，但数字完全不可信。
Python免费学习笔记（深入）：立即使用
在学习笔记中，你将探索 Python 的核心概念和高级技巧！
