In this paper, we discuss, both analytically and numerically, the paraxial propagation of the radially polarized Laguerre-Gaussian-correlated Schell-model (LGCSM) beams orthogonal to the optical axis in uniaxial crystals. The analytical expression for the cross-spectral density function and the second-order moments of the radially polarized LGCSM beams are derived, and the evolution properties of the normalized intensity distribution, the spectral degree of the coherence (SDOC), and the spectral degree of the polarization (SDOP) in uniaxial crystals are elucidated by numerical examples. It is found that the intensity distribution of the radially polarized LGCSM beams evolves from a doughnut shape into a solid shape and finally converts into an elliptical symmetric hollow-ring profile in uniaxial crystals due to the combined effect of special correlation functions and the anisotropy effect of the uniaxial crystals. The evolution of the SDOC and SDOP for the radially polarized LGCSM beams is quite different from that of the radially polarized Gaussian-Schell-model beams. In addition, the propagation properties of the radially polarized LGCSM beams are closely related to the spatial coherence length, the mode order, and the ratio of extraordinary and ordinary reflective indices. The results show that the uniaxial crystals could modulate the evolution properties of the radially polarized LGCSM beams.