I collected 10 NBA teams of the 2017 season  and the following variables:

W – the number of games won
L – the number of games lost
P – the number of points (or runs, goals, etc) scored by the team
PA – the number of points allowed by the team

Data displayed in the table below:

The Plot of log(W/L) against log(P/PA) with the best fitted line of the form k log(P/PA) and Residuals against log(P/PA) is presented below: The formula of best fit is given by “lm” method as: log(W/L)=13.76*log(PA/P) which implies the corresponding Pythagorean formula is W/L=(PA/P)^13.76.

Lucky or unlucky, there is no apparent unusual observation in the residual plot, all the observations are bounded by 0.2 log(W/L) away from the fitted line, and the high leverage observation which is to the very left also fits well with the formula.