Billiards are flat cavities where a particle is free to move between elastic collisions with the boundary. In chaos theory, these systems are simple prototypes. The conservative dynamics of a billiard may vary from regular to chaotic, depending only on the shape of the boundary. This work aims to shed light into the quantization of classically chaotic systems. We present numerical results on classical and quantum properties in two bi-parametric families of billiards, namely Elliptical Stadium Billiard (ESB) and Elliptical-C3 Billiards (E-C3B). Both families entail elliptical perturbations of chaotic billiards with originally circular sectors on their borders. Our numerical calculations provide evidence that these elliptical families may exhibit a mixed classical dynamics, identified by the chaotic fraction of the phase space, the parameter χc<1. We use this quantity to guide our analysis of quantum spectra. We explore the short-range correlations through nearest neighbor spacing distribution p(s), revealing that in the mixed region of the classical phase space, p(s) is well described by the Berry–Robnik–Brody (BRB) distributions for the ESB. In agreement with the expectation from the so-called ergodic parameter α=tH/tT, the ratio between the Heisenberg time and the classical diffusive-like transport time. Our findings indicate the possibility of quantum dynamical localization when α<1. For the E-C3B family, the eigenstates can be split into singlets and doublets. BRB successfully describes p(s) for singlets as the previous family in the mixed region. However, for doublets, new distributions recently introduced in the literature come into play, providing descriptions for p(s) with a focus on cases where χc=1. We observed that as χc decreases, the p(s) distributions simultaneously deviate from the Gaussian Orthogonal Ensemble (GOE) for singlets, and Gaussian Unitary Ensemble (GUE) for doublets.