Variations in scaling behavior in the flux and emissions of distant astronomical sources with respect to their cosmic time are important l phenomena that can provide valuable information about the dynamics within the sources and their cosmological evolution with time. Different studies have been applying linear analysis to understand and model quasars' light curves. Here, we study the multifractal behavior of selected quasars' radio emissions in their observed frame (at 22 and 37 GHz bands) and and their rest frame. To this end, we apply the wavelet transform-based multifractal analysis formalism called wavelet transform modulus maxima. In addition, we verify whether the autoregressive integrated moving average (ARIMA) models fit our data or not. In our work, we observe strong multifractal behavior for all the sources. Additionally, we find that the degree of multifractality is strongly similar for each source and significantly different between sources at 22 and 37 GHz. This similarity implies that the two frequencies have the same radiation region and mechanism, whereas the difference indicates that the sources have intrinsically different dynamics. Furthermore, we show that the degree of multifractality is the same in the observed and rest frames of the quasars, i.e., multifractality is an intrinsic property of radio quasars. Finally, we show that the ARIMA models fit the 3C 345 quasar at 22 GHz and partially fit most of the time series with the exception of the 3C 273 and 3C 279 quasars at 37 GHz, for which the models are found to be inadequate.