On the Sampling-based Computation of Nash Equilibria under Uncertainty via the Nikaido-Isoda Function
Luke Marrinan, Farzad Yousefian, Uday V. Shanbhag
TL;DR
The paper addresses computing Nash equilibria in stochastic NEPs with convex, Lipschitz objectives by reformulating the problem with the Nikaido-Isoda (NI) function and its regularized variant $V_\alpha$. It develops a sampling-enabled, inexact projected gradient method that handles inexact best-responses and uses stochastic gradient estimates to minimize $V_\alpha$, proving almost-sure convergence to stationary points of $V_\alpha$ and deriving explicit rate and complexity guarantees. Crucially, under a suitable regularity condition, stationary points of $V_\alpha$ correspond to Nash equilibria, enabling equilibrium computation without requiring gradient map monotonicity or game potentiality. The results yield concrete projection and sampling complexities: with a fixed stepsize, $\mathcal{O}(\varepsilon^{-2})$ projections and $\mathcal{O}(\varepsilon^{-4})$ samples; with diminishing stepsizes, $\mathcal{O}(\varepsilon^{-4})$ projections and $\mathcal{O}(\varepsilon^{-6})$ samples, offering practical guarantees for stochastic equilibrium computation in uncertain multi-agent settings.
Abstract
We consider the computation of an equilibrium of a stochastic Nash equilibrium problem, where the player objectives are assumed to be $L_0$-Lipschitz continuous and convex given rival decisions with convex and closed player-specific feasibility sets. To address this problem, we consider minimizing a suitably defined value function associated with the Nikaido-Isoda function. Such an avenue does not necessitate either monotonicity properties of the concatenated gradient map or potentiality requirements on the game but does require a suitable regularity requirement under which a stationary point is a Nash equilibrium. We design and analyze a sampling-enabled projected gradient descent-type method, reliant on inexact resolution of a player-level best-response subproblem. By deriving suitable Lipschitzian guarantees on the value function, we derive both asymptotic guarantees for the sequence of iterates as well as rate and complexity guarantees for computing a stationary point by appropriate choices of the sampling rate and inexactness sequence.
