今天使用Julia计算Pearson系数,本想用自己写的函数,结果发现自带的Statistics库能直接计算,如图所示。
如果我们有一个10X5的X矩阵,如何计算每列间的Pearson系数呢?直接使用cor(X, dims=1)。若想计算行之间的Pearson系数就将dims设置为2。
计算后得到相关系数矩阵Data,此时使用Makie中的heatmap即可绘制好看的图片了。因为Data中系数都在0~1之间,将colorrange设置为(0,0.3)画图要好看写,colormap设置为magma。