r_squared = r2_score(y_test, predictions)

1 - (1 - r_squared) * (len(y_test) - 1) / (len(y_test) - X_test.shape[1] - 1)