Я пытаюсь вычислить корреляцию Пирсона по столбцам, используя 2 двумерных массива размером 27495x92. I sh, чтобы получить результирующую матрицу корреляции размером 92х92. Я попробовал один способ, как показано ниже, но значения больше 1. Мне было интересно, если это проблема нормализации, и я был бы благодарен, если бы кто-нибудь мог предложить более простой способ сделать то же самое с функцией numpy. Любая помощь приветствуется. Спасибо!
trans_rsfMRI_corr_2d = (trans_rsfMRI_corr_2d - trans_rsfMRI_corr_2d.mean(axis = 0)) / trans_rsfMRI_corr_2d.std(axis = 0)
trans_mvfMRI_corr_2d = (trans_mvfMRI_corr_2d - trans_mvfMRI_corr_2d.mean(axis = 0)) / trans_mvfMRI_corr_2d.std(axis = 0)
rs_mv_corr = np.dot(trans_rsfMRI_corr_2d.T, trans_mvfMRI_corr_2d) / trans_rsfMRI_corr_2d.shape[0]`