Sensibility analysis
The final estimations, with their confidence intervals, of YLD and YLL, were made by performing 10 000 iterations in a Monte Carlo simulation of each parameter using the bootstrapped technique of DALYs calculator package in R program (18). All estimates are performed using a discount rate of 3% and weighting by age (C = 0.1658, β = 0.04). Subgroup analyzes were performed according to age group and sex. The probabilistic sensitivity analysis was performed using the entire range of initial values ​​of each of the parameters used to calculate the DALYs using the standardized regression coefficient method. The simulated overall DALY estimates were regressed against the simulated values for the stochastic input parameters. To facilitate comparison, the independent terms are standardized such that they were normally distributed with mean zero and standard deviation one. The resulting regression coefficients were therefore referred to as standardized regression coefficients. The regression coefficients correspond to the number of standard deviations change in overall DALY given one standard deviation change in the corresponding input parameter. The internal consistency of each parameter was evaluated using the DISMOD II program (19) following the recommendations and methodology described in the manual for national studies of the WHO disease burden (8).