我有两个具有形状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要求输入数组为一维.