Statistical analyses of the spectral distributions of rotational bands in 51 deformed prolate even–even nuclei in the 152 ≤ A ≤ 250 mass region $$R_{{4_{1}^{ + } /2_{1}^{ + } }} \ge 3.00$$
R
4
1
+
/
2
1
+
≥
3.00
are examined in terms of nearest neighbor spacing distributions. Specifically, the focus is on data for 0+, 2+, and 4+ energy levels of the ground, gamma, and beta bands. The chaotic behavior of the gamma band, especially the position of the $$2_{\gamma }^{ + }$$
2
γ
+
band-head compared to other levels and bands, is clear. The levels are analyzed within the framework of two models, namely, a SU(3)-partial dynamical symmetry Hamiltonian and a SU(3) two-coupled quasi-dynamical symmetry Hamiltonian, with results that are further analyzed using random matrix theory. The partial and quasi dynamics both yield outcomes that are in reasonable agreement with the known experimental results. However, due to the degeneracy of the beta and gamma bands within the simplest SU(3) picture, the theory cannot be used to describe the fluctuation properties of excited bands. By changing relative weights of the different terms in the partial and quasi dynamical Hamiltonians, results are obtained that show more GOE-like statistics in the partial dynamical formalism as the strength of the pairing term is increased. Also, in the quasi-dynamical symmetry limit, more correlations are found because of the stronger couplings.