《使用python对两组被试功能连接相关矩阵做t检验》有一个想法

  1. 使用scipy.spatial.distance.squareform函数,可以简单地取其上三角矩阵拉直后的vector:
    from scipy.spatial.distance import squareform as sqf
    cc_vec = sqf(cc_mat)

    但是使用np.corrcoef()函数直接算出的矩阵cc_mat,sqf函数可能会将其识别为非对称:
    ValueError: Distance matrix ‘X’ must be symmetric.

    此时,可以使用cc_mat.astype(float32),减少一些精度,就可以正常使用sqf函数了。
    比文中操作上三角矩阵的方法简单很多。

发表评论

电子邮件地址不会被公开。 必填项已用*标注