Я пытаюсь выяснить, как рассчитать ковариацию с помощью функции Python Numpy cov. Когда я передаю ему два одномерных массива, я получаю матрицу результатов 2х2. Я не знаю, что с этим делать. Я не силен в статистике, но я считаю, что ковариация в такой ситуации должна быть единым числом.Этот это то, что я ищу. Я написал свой собственный:

def cov(a, b):

    if len(a) != len(b):
        return

    a_mean = np.mean(a)
    b_mean = np.mean(b)

    sum = 0

    for i in range(0, len(a)):
        sum += ((a[i,] - a_mean) * (b[i] - b_mean))

    return sum/(len(a)-1)

Это работает, но я полагаю, что версия Numpy гораздо эффективнее, если бы я мог понять, как ее использовать.

Кто-нибудь знает, как заставить функцию Numpy cov работать так, как я написал?

Спасибо,

Дейв