In order to identify the bus drivers who are about to cause accidents, the data set was obtained by combining the Bus Safety Management System database, the Baidu application programming interface(API) and web crawling technology. K-Nearest Neighbor algorithm was used to supplement the missing values and data from 1893 drivers in 42 lines was obtained. The basic characteristic variables included driver, vehicle, route characteristics, violations, accidents, management, and further construct derived characteristics on this basis. An integrated method, including recursive feature elimination, logistic regression with penalty terms, random forest and others, was designed and used for feature selection. The model was built using 6 machine methods like XGBoost and optimized for the hyper-parameters using Bayesian methods. The results indicate that among the six classification models constructed, the model constructed by XGBoost method has the best area under receiver operating characteristic(ROC) area under curve(AUC) evaluation results. Bayesian optimization can improve the AUC of ROC to a certain extent. For the accident driver's prediction accuracy rate reaches 98.66%, the operating unit can also weigh the false positive rate and true positive rate according to its own situation. Moreover, the nonlinear influence effect of features is found in the model results. The characteristics of vehicle service time, driving age, violations, punishment and other characteristics have a very obvious role in the accident risk.
