The discrepancy arises because the user is calculating R-squared for univariate regressions incorrectly by subtracting the mean of Y, which is not appropriate for models without an intercept. To align with the geometric interpretation, R-squared should be computed as the squared ratio of the norm of the predicted values to the norm of Y. Here's how to fix it:\n\n1. Replace costheta1 = np.linalg.norm(Y1 - np.mean(Y)) / np.linalg.norm(Y - np.mean(Y)) with costheta1 = np.linalg.norm(Y1) / np.linalg.norm(Y).\n2. Similarly, replace costheta2 = ... with costheta2 = np.linalg.norm(Y2) / np.linalg.norm(Y).\n\nThis adjustment ensures that the calculated R-squared values (cos²(theta)) for univariate regressions match the angles (45° and 120°) set in the script.