计算两个多维数组之间的相关系数

Question

计算两个多维数组之间的相关系数

我有两个具有形状N X T和数组的阵列M X T。我想计算T每对可能的行n和m（分别来自N和M）之间的相关系数。

什么是最快，最pythonic的方式来做到这一点？（循环过来N，M在我看来既不快也不是pythonic。）我期待答案涉及numpy和/或scipy。现在我的数组是numpy arrays，但我愿意将它们转换为不同的类型。

我期待我的输出是一个形状的数组N X M。

NB当我说“相关系数”时，我的意思是Pearson乘积矩相关系数。

以下是一些需要注意的事项：

• 该numpy函数correlate要求输入数组为一维。

• 该numpy函数corrcoef接受二维数组，但它们必须具有相同的形状。

• 该scipy.stats函数pearsonr要求输入数组为一维。