|
|
|
|
@ -22,19 +22,19 @@ okz2.py
|
|
|
|
|
```
|
|
|
|
|
import math
|
|
|
|
|
def correlation(l1, l2):
|
|
|
|
|
n = min(len(l1), len(l2))
|
|
|
|
|
n = min(len(l1), len(l2))#определяю длину для расчёта — берётся минимальная длина из двух списков
|
|
|
|
|
sum1 = sum(l1)
|
|
|
|
|
sum2 = sum(l2)
|
|
|
|
|
sum1sq = sum(now ** 2 for now in l1)
|
|
|
|
|
sum1sq = sum(now ** 2 for now in l1)#суммы квадратов элементов каждого списка
|
|
|
|
|
sum2sq = sum(now ** 2 for now in l2)
|
|
|
|
|
sum12 = sum(x * y for x, y in zip(l1, l2))
|
|
|
|
|
|
|
|
|
|
part1 = n * sum12 - sum1 * sum2
|
|
|
|
|
part2 = math.sqrt((n * sum1sq - sum1 ** 2) * (n * sum2sq - sum2 ** 2))
|
|
|
|
|
part2 = math.sqrt((n * sum1sq - sum1 ** 2) * (n * sum2sq - sum2 ** 2))#расчёт числителя и знаменателя формулы корреляции Пирсона
|
|
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
corr = part1 / part2
|
|
|
|
|
return corr
|
|
|
|
|
return corr#коэф корреляции
|
|
|
|
|
except ZeroDivisionError:
|
|
|
|
|
print("Ошибка деления на ноль!")
|
|
|
|
|
return
|
|
|
|
|
|