The hole expansion ratio (HER) observed in a standardized hole expansion test (HET) is commonly used to determine the edge fracture of steel sheets. A large variation of the measured HER restricts the practical application of the method. This study presents a systematic investigation on uncertainties in the HER of DP800 sheet material, including the hole-edge quality, pre-strain due to the hole-punching process, the friction coefficient, and the determination of fracture. An artificial neural network was trained to develop a surrogate model using a database gained from a thousand finite element simulations of the HET. Monte-Carlo simulations were performed using the trained surrogate model to characterize the distribution of the HER. Sensitivity analysis via Sobol's indices is calculated to determine the influence of the input variables on the output. It is found that the pre-strain and pre-damage generated during the hole punching process in the shear-affected zone dominate the variation of the HER. Discussions on reducing the output's variation are detailed. In general, these findings provide valuable insights for the determination of HER as well as the edge crack behavior of the investigated DP800 steel sheet.