在研究脑网络时,看到了用皮尔逊相关系数计算脑功能网络的,就简单的学了一下它的计算
这里用到的是Python的numpy库
- 首先,最前面要加import,先导入这个模块
import numpy as np
- 然后,调用corrcoef()
其中data是原始矩阵,pcc_matrix是生成的pcc矩阵
比如data为32*1280的矩阵,pcc_mattix则为32*32的矩阵
每个元素值的范围为[-1,1]
pcc_matrix = np.corrcoef(data)
- 如果要进一步处理这个矩阵,则使用下面的代码
作用是将元素值的范围映射到[0,1]
pcc_matrix = np.corrcoef(data) * 0.5 + 0.5
pcc_matrix矩阵是一个对称矩阵,下面是实验的部分截图: