Introduction: Accurate assessment of the risk for bronchopulmonary dysplasia (BPD) is critical to determine the prognosis and identify infants who will benefit from preventive therapies. Clinical prediction models can support the identification of high-risk patients. In this study, we investigated the potential risk factors for BPD and compared machine learning models for predicting the outcome of BPD/death on days 1, 7, 14, and 28 in preterm infants. We also developed a local BPD estimator. Methods: This study involved 124 infants. We evaluated the composite outcome of BPD/death at a postmenstrual age of 36 weeks and identified risk factors that would improve BPD/death prediction. SPSS for Windows Version 11.5 and Weka 3.9 software were used for the data analysis. Results: To evaluate the combined effect of all variables, all risk factors were taken into consideration. Gestational age, birth weight, mode of respiratory support, intraventricular hemorrhage, necrotizing enterocolitis, surfactant requirement, and late-onset sepsis were risk factors on postnatal days 7, 14, and 28. In a comparison of four different time points (postnatal days 1, 7, 14, and 28), the day 7 model provided the best prediction. According to this model, when a patient was diagnosed with BPD/death, the accuracy rate was 89.5%. Conclusion: The postnatal day 7 model was the best predictor of BPD or death. Future validation studies will help identify infants who may benefit from preventive therapies and develop individualized care.