The exact or Wilson renormalization group equations can be formulated as a functional Fokker–Planck equation in the infinite-dimensional configuration space of a field theory, suggesting a stochastic process in the space of couplings. Indeed, the ordinary renormalization group differential equations can be supplemented with noise, making them stochastic Langevin equations. Furthermore, if the renormalization group is a gradient flow, the space of couplings can be endowed with a supersymmetric structure a la Parisi–Sourlas. The formulation of the renormalization group as supersymmetric quantum mechanics is useful for analysing the topology of the space of couplings by means of Morse theory. We present simple examples with one or two couplings.