The discrepancy is because you need to use the running masses of the bottom quark and muon at the 125 GeV energy scale, rather than the pole masses of those particles in the formula that you cite (see generally here).
The running mass of the bottom quark at the 125 GeV energy scale is roughly 3.0 GeV (v. 4.18 GeV pole mass) and the running mass of the muon at that energy scale is about .104 GeV (v. a pole mass of .1057 GeV). See this 1997 paper, and then adjust for the differences between the current pole mass measurements and the ones used in the 1997 paper and for the closest scale they calculate (the Z boson mass energy scale) and the 125 GeV target energy scale, realizing that the running of masses in the Standard Model with energy scale is roughly logarithmic.
The 125 GeV energy scale is special because the experimentally measured value of the Higgs boson mass is currently 125.18 GeV +/ 0.16 GeV.
I get a ratio of 2543 calculating from the table at 125.0 GeV, and I get a ratio of 2488 using the running values of these masses, which is the same to within the uncertainty of the inputs used in the respective calculations (my calculations were made using Microsoft Xcel, the formula provided in the question with a factor of three adjustment for the three b quark colors, and data noted in this answer).
All currently measured masses are from the live edition of the Particle Data Group's website.