An Iterative Bayesian Robbins--Monro Sequence

This study introduces an iterative Bayesian Robbins--Monro (IBRM) sequence, which unites the classical Robbins--Monro sequence with statistical estimation for faster root-finding under noisy observations. Although the standard Robbins--Monro method iteratively approaches solutions, its convergence speed is limited by noisy measurements and naivety to any prior information about the objective function. The proposed Bayesian sequence dynamically updates the prior distribution with newly obtained observations to accelerate convergence rates and robustness. The paper demonstrates almost sure convergence of the sequence and analyses its convergence rates for both one-dimensional and multi-dimensional problems. We evaluate the method in a practical application that suffers from large variability and allows only a few function evaluations, specifically estimating thresholds in noninvasive brain stimulation, where the method is more robust and accurate than conventional alternatives. Simulations involving 25,000 virtual subjects illustrate reduced error margins and decreased outlier frequency with direct impact on clinical use.