theses.hal.science
Divergence-minimization for variational inference, black-box global optimization, and importance sampling
Many methods across applied mathematics aim at constructing specific parametric probability distributions. Examples of these tasks include evolution strategies or simulated annealing for black-box global optimization, Monte Carlo methods based on adaptive importance sampling, and variational inference algorithms in machine learning. The construction of such distributions can often be formulated as the minimization of a statistical divergence. However, these divergence-minimization problems are challenging because of the following reasons. First, the specific geometry of the considered set of parametric probability distributions needs to be taken into account. Second, efficient evaluations of statistical divergences come with important noise. Third, divergence-minimization problems are generally non-convex. Because of these difficulties, standard methods may fail to converge to good solutions. We tackle these challenges in this thesis.First, we show that evolution strategies, which are sampling-based algorithms for black-box global optimization problems, can be analysed through the lens of divergence-minimization problems. Our approach allows to establish and quantify the improvement brought at each iteration of the algorithms. We show that existing methods fit within our framework, yielding a new approach for their analysis. We also establish improvement results for two novel algorithms, one related with mixture models, and another one using heavy-tailed parametric probability distributions.Second, we consider the minimization of a regularized Rényi divergence over an exponential family. We propose to solve this problem with a stochastic Bregman proximal-gradient algorithm, with biased gradient estimator. By leveraging the geometry of the exponential family, we prove strong convergence guarantees for our algorithm, with proof techniques that are of interest beyond the considered problem. We then extend this algorithm to propose an adaptive simulated annealing algorithm with solid theoretical understanding. We show through a rigorous benchmarking that our algorithm outperforms similar non-adaptive algorithms.Finally, we go beyond exponential families and look at variational inference problems over lambda-exponential families. Using generalized convexity tools, we give new sufficient optimality conditions for these problems, which generalize existing similar results for the exponential family. For the resolution of these problems, we propose novel proximal-like algorithms that exploit the geometry underlying the lambda-exponential family. These results are especially useful for heavy-tailed distributions. We then leverage our results to propose an adaptive importance sampling algorithm to cover these cases. We show that our algorithm is able to learn Student distributions that capture the location, scale, and tail behaviour of target distributions, both in heavy-tailed and light-tailed cases.