import re, numpy as np
from scipy import stats
s = """сюда вставь весь список"""
pairs = re.findall(r'\(([^;]+);([^)]+)\)', s)
xy = [(float(a), float(b)) for a, b in pairs if a != 'NA' and b != 'NA']
x = np.array([i[0] for i in xy])
y = np.array([i[1] for i in xy])
n = len(x)
r = np.corrcoef(x, y)[0, 1]
t = r * ((n - 2) / (1 - r**2))**0.5
p = 2 * (1 - stats.t.cdf(abs(t), n - 2))
print(n)
print(f"{r:.8f}")
print(f"{t:.8f}")
print(f"{p:.8f}")
print("Да" if p < 0.04 else "Нет")