n = 10
x = np.arange(n)
y = 3 * x + 5 + np.random.random(n)
A = np.vstack((x, np.ones(n))).T
model, resid = np.linalg.lstsq(A, y)[:2]
r2 = 1 - resid / (y.size * y.var())
print(r2)
Преобразование остаточного значения Numpy Lstsq в R^2