We investigate the charging energy level statistics of disordered interacting electrons in quantum dots by numerical calculations using the Hartree approximation. The aim is to obtain a global picture of the statistics as a function of disorder and interaction strengths. We find Poisson statistics at very strong disorder, Wigner-Dyson statistics for weak disorder and interactions, and a Gaussian intermediate regime. These regimes are as expected from previous studies and fundamental considerations, but we also find interesting and rather broad crossover regimes. In particular, intermediate between the Gaussian and Poisson regimes we find a two-sided exponential distribution for the energy level spacings. In comparing with experiment, we find that this distribution may be realized in some quantum dots.