The purpose of this paper is to demonstrate how modern statistical techniques of non-stationary time-series analysis can be used to characterize the mutual interaction among three coexisting rhythms in nephron pressure and flow regulation. Besides a relatively fast vasomotoric rhythm with a period of 5-8 s and a somewhat slower mode arising from an instability in the tubuloglomerular feedback mechanism, we also observe a very slow mode with a period of 100-200 s. Double-wavelet techniques are used to study how the very slow rhythm influences the two faster modes. In a broader perspective, the paper emphasizes the significance of complex dynamic phenomena in the normal and pathological function of physiological systems and discusses how simulation methods can help to understand the underlying biological mechanisms. At the present there is no causal explanation of the very slow mode. However, vascular oscillations with similar frequencies have been observed in other tissues.