Voltage sags are power quality disturbances mainly caused by faults with a highly random nature. For this reason, the severity and frequency of occurrence of voltage sags is typically analyzed by means of probabilistic methods, such as Monte Carlo simulations, where different probability distributions are selected to model the uncertain inputs (e.g., fault impedance, fault location, fault rate, etc.). The selection and parameterization of the appropriate probability functions involves many difficulties to ensure that they represent the distribution of input variables in a realistic way and so, many different approaches can be found in literature. This paper studies the sensitivity of voltage sags severity to the uncertainty in the input models in order to focus modeling effort on the input variables that are more critical and effectively have an influence on voltage sags levels, whereas non-influential parameters can be neglected or modeled by means of simplified hypothesis. The method proposed is based on Design of Experiments and Analysis of Variance (ANOVA) and is able to discriminate with a statistic confidence level whether the uncertainty in the input significantly affects voltage sags severity indices or not. The effect of uncertainty on voltage sags indices SARFI90 and SARFI70 has been evaluated in realistic scenarios for IEEE test networks of 24 and 118 buses and in the power system of Ecuador with 357 nodes. General trends have been also established that help to understand the effect of the modeling of input parameters on the estimated number of sags.INDEX TERMS Analysis of variance, design of experiments, sensitivity analysis, power quality, uncertainty, voltage sag.