The prediction of Jominy hardness curves and the effect of alloying elements on the hardenability of boron steels (19 different steels) are investigated using multiple regression analysis. To evaluate the hardenability of boron steels, Jominy end quenching tests are performed. Regardless of the alloy type, lath martensite structure is observed at the quenching end, and ferrite and pearlite structures are detected in the core. Some bainite microstructure also appears in areas where hardness is sharply reduced. Through multiple regression analysis method, the average multiplying factor (regression coefficient) for each alloying element is derived. As a result, B is found to be 6308.6, C is 71.5, Si is 59.4, Mn is 25.5, Ti is 13.8, and Cr is 24.5. The valid concentration ranges of the main alloying elements are 19 ppm < B < 28 ppm, 0.17 < C < 0.27 wt%, 0.19 < Si < 0.30 wt%, 0.75 < Mn < 1.15 wt%, 0.15 < Cr < 0.82 wt%, and 3 < N < 7 ppm. It is possible to predict changes of hardenability and hardness curves based on the above method. In the validation results of the multiple regression analysis, it is confirmed that the measured hardness values are within the error range of the predicted curves, regardless of alloy type.