The significant improvement of MAE to reaction energies for the PBE functional compared to B3LYP-D or M06-2X indicates that the bond-based correction method relying on linear regression is better at capturing systematic errors in reaction energy (from less pa-rameterized DFT functionals) compared to the highly parameterized ones. The significant reduction in number of reactions with higher than 0.25 eV error assures the effectiveness of this method for controlling large errors. The computational cost of the proposed method is negligible and it can be easily implemented in simulation libraries such as the atomistic simulation environment (ASE). Although the correction terms were calculated using the errors in reaction energies from post-SCF DFT, the application to real systems shows that it can be equally applied to correct reaction energies from regular DFT. Lastly, the method can be extended to other datasets to include a more diverse set of bonds compared to the ones considered in this study.