Archive for PariSanté campus

mostly Monte Carlo [20/02]

Posted in Statistics with tags , , , , , , , , , , , , , , , , on February 16, 2026 by xi'an

A new episode of our mostly Monte Carlo seminar, very soon coming near you (if in Paris):

On Friday 20/02/26, from 3-5pm at PariSanté Campus

15h: Paul Mangold (École Polytechnique, Palaiseau)

Convergence and Linear Speed-Up in Stochastic Federated Learning

In federated learning, multiple users collaboratively train a machine learning model without sharing local data. To reduce communication, users perform multiple local stochastic gradient steps that are then aggregated by a central server. However, due to data heterogeneity, local training introduces bias. In this talk, I will present a novel interpretation of the Federated Averaging algorithm, establishing its convergence to a stationary distribution. By analyzing this distribution, we show that the bias consists of two components: one due to heterogeneity and another due to gradient stochasticity. I will then extend this analysis to the Scaffold algorithm, demonstrating that it effectively mitigates heterogeneity bias but not stochasticity bias. Finally, we show that both algorithms achieve linear speed-up in the number of agents, a key property in federated stochastic optimization.

16h: Alain Durmus (École Polytechnique, Palaiseau)

A Mixture-based Framework for Guiding Diffusion Models

Inverse problems—such as image restoration from noisy or incomplete measurements and musical source separation—are ill-posed, making Bayesian approaches with learned generative priors especially appealing. Diffusion models provide powerful priors, but existing posterior sampling methods often rely on crude likelihood-gradient approximations and heavy task-specific tuning. In this talk, I will introduce a novel principled approach specifically designed to overcome these limitations. The core contribution of this approach is the construction of a mixture approximation of intermediate posterior distributions defined by the diffusion model. The sampling is carried out sequentially via Gibbs sampling, a Markov Chain Monte Carlo method, using a careful data augmentation scheme. Gibbs sampling is employed here due to its simplicity and theoretical guarantees, allowing for exact conditional updates at each iteration, thus ensuring stability and efficiency. One key advantage of the presented algorithm is its flexibility: it adapts to varying levels of computational resources by adjusting the number of Gibbs iterations. Consequently, substantial performance gains can be achieved by increasing inference-time computational effort. I will present extensive experimental results demonstrating empirical performance across diverse image restoration tasks, involving both pixel-space and latent-space diffusion models, and showcase its successful application in musical source separation.

 

January session of the mostly Monte Carlo seminar (16/01, 3pm)

Posted in Statistics, University life with tags , , , , , , , , , , , on January 9, 2026 by xi'an

corrida de Noël, Issy-les-Moulineaux [10K, 41’00, 701/3335, 1/38 M6M, 4´06/km]

Posted in pictures, Running, Travel with tags , , , , , , , , , , on December 16, 2025 by xi'an

mostly Monte Carlo [last session of 2025]

Posted in Books, pictures, Statistics, Travel, University life with tags , , , , , , , , , , on December 10, 2025 by xi'an

Rather hurriedly, here is the announcement for the last mostly MC seminar this year, to take place at 3-5pm this very Friday, Dec 12, 2025. It will take place in Salle 03, PariSanté Campus


3pm Least squares variational inference

Yvann Le Fay CREST, ENSAE

Variational inference seeks the best approximation of a target distribution within a chosen family, where “best” means minimizing Kullback-Leibler divergence. When the approximation family is exponential, the optimal approximation satisfies a fixed-point equation. We introduce LSVI (Least Squares Variational Inference), a gradient-free, Monte Carlo-based scheme for the fixed-point recursion, where each iteration boils down to performing ordinary least squares regression on tempered log-target evaluations under the variational approximation. We show that LSVI is equivalent to biased stochastic natural gradient descent and use this to derive convergence rates with respect to the numbers of samples and iterations. When the approximation family is Gaussian, LSVI involves inverting the Fisher information matrix, whose size grows quadratically with dimension d. We exploit the regression formulation to eliminate the need for this inversion, yielding O(d³) complexity in the full-covariance case and O(d) in the mean-field case. Finally, we numerically demonstrate LSVI’s performance on various tasks, including logistic regression, discrete variable selection, and Bayesian synthetic likelihood, showing competitive results with state-of-the-art methods, even when gradients are unavailable.

4pm Beyond the Unified Skew-Normal: Extended Models for Bayesian Classification

Paolo Onorati CEREMADE, Université Paris Dauphine – PSL

Binary classification models typically lose the conjugacy and computational simplicity enjoyed by Gaussian models. While the Unified Skew-Normal (SUN) family has recently been shown to be conjugated under the probit model, two new developments are presented that extend this idea to a broader class of link functions, including both logit and probit. In the parametric setting, the Perturbed Unified Skew-Normal (pSUN) distribution is introduced; it is conjugate to any binary regression model whose link admits a scale-mixture representation of Gaussian random variables, enabling tractable posterior summaries, efficient sampling schemes, and strong performance in high-dimensional covariate settings. The discussion then moves to the nonparametric domain, where the Quasi SUN family and the associated stochastic process provide conjugacy for nonparametric logit and probit models while preserving key closure properties. A stochastic representation of this process yields practical computational improvements over existing Gaussian-based approximations. Together, these SUN-type extensions offer promising tools for Bayesian classification with accurate posterior inference.

mostly Monte Carlo, November

Posted in pictures, Statistics, Travel, University life with tags , , , , , , , , , , , on November 7, 2025 by xi'an

The November session of the Mostly (and monthly) Monte Carlo seminar will take place next week on Thursday, November 13, 2025, at 3PM in Salle 08, PariSanté Campus.  With two exciting speakers:

Abstracts are available on the seminar’s website