Owing to the continuous deterioration in the quality of iron ore and scrap, there is an increasing focus on improving the Basic Oxygen Furnace (BOF) process to utilize lower grade input materials. The present paper discusses dephosphorization in BOF steelmaking from a data science perspective, which thus enables steelmakers to produce medium and low phosphorus steel grades. In the present study, data from two steel mills (Plant I and Plant II) were collected and various statistical methods were employed to analyze the data. While most operators in steel plants use spreadsheet-based techniques and linear regression to analyze data, this paper discusses on the suitability of selecting various statistical methods, and benchmarking tests to analyze such dephosphorization data sets. The data contains a wide range of operating conditions, both low and high phosphorus input loads, different slag basicity’s, different slag chemistries, and different end point temperatures, etc. The predicted phosphorus partition from various statistical models is compared against plant data and verified against previously published research.