My guess would be that it does correlate, even though it is not a 1:1 correlation. I would use a linear regression, including a ratio variable for each base stat, plus an additional one for total stats. I would probably also add one for the sum of each Pokemon's higher attack stat and speed, plus one for the sum of both defenses and HP. Then, I would add a binary dummy variable for every type and one for every ability. This might create some sample size problems, particularly in cases of extremely rare abilities, like Hydration, or those mostly given to unusually powerful Pokemon (I'm thinking specifically of Pressure). This list would of course not contain the distorting stats of the Ubers, since it would be based on OU. NFEs would also be omitted (except the totally different ones, like Vigoroth). The variable of interest would be usage in OU.
Of course, this would not consider movepool. Including a dummy variable for every single move, in addition to being very time-consuming, would not work due to massive multicollinearity. It might be possible to include binary variables just for a few specific coverage combinations normally considered to be powerful, such as Ice/Electric, Water/Ice/Grass, and Fire/Dragon. There could also be one more variable for "outclassed," which would be determined via consensus among high level players; it would be "1" for anything considered to have a flat-out superior counterpart, such as Blaziken, Rotom basic, etc, and "0" for everything else. If I were running the data, I would try it with and without this variable to see if it distorts any of the others, and my guess would be that it does have a statistically significant relationship with usage.