---
title: Performative Stability in Adaptive Prediction
url: https://www.emergentmind.com/topics/performative-stability
type: topic
---

# Performative Stability in Adaptive Prediction

Performative stability is the central equilibrium concept in performative prediction and its extensions. It formalizes settings in which a deployed model, policy, or decision changes the very distribution or environment on which it is later evaluated. In its canonical form, one has a parameter space $\Theta$, a distribution map $D:\Theta\to\Delta(\mathcal Z)$, and loss $\ell(z;\theta)$. The deployed parameter $\theta$ induces data from $D(\theta)$, and a parameter is performatively stable if retraining on its own induced distribution returns the same parameter. In this sense, performative stability is a fixed point of a retraining or best-response map, and it differs from ordinary empirical risk minimization precisely because the data-generating process is endogenous to deployment [2002.06673].

## 1. Foundations and solution concepts

The original single-agent formulation distinguishes between the performative risk
$$
PR(\theta)=\mathbb E_{Z\sim D(\theta)}[\ell(Z;\theta)]
$$
and the decoupled or cross-risk
$$
R(\theta_{\mathrm{induce}},\theta_{\mathrm{eval}})
=\mathbb E_{Z\sim D(\theta_{\mathrm{induce}})}[\ell(Z;\theta_{\mathrm{eval}})].
$$
A parameter $\theta^\dagger$ is performatively stable when
$$
\theta^\dagger\in\arg\min_{\theta'}R(\theta^\dagger,\theta')
\quad\Longleftrightarrow\quad
\theta^\dagger=\arg\min_{\theta'}\mathbb E_{Z\sim D(\theta^\dagger)}[\ell(Z;\theta')],
$$
equivalently, a fixed point of
$$
G(\theta)=\arg\min_{\theta'}R(\theta,\theta').
$$
The concept was introduced as the equilibrium notion appropriate for learning systems whose predictions reshape future outcomes [2002.06673].

Performative stability is distinct from performative optimality. A performative optimum minimizes the self-referential objective
$$
\theta^\star\in\arg\min_{\theta}PR(\theta)=\arg\min_{\theta}R(\theta,\theta),
$$
whereas a performatively stable point only requires self-consistency under retraining. The two notions need not coincide. The survey literature emphasizes that stability captures “stability under induced shifts”: once the environment has reacted to $\theta_{PS}$, the best counterfactual model for the resulting distribution is $\theta_{PS}$ itself [2602.10176].

The same fixed-point structure recurs in later generalizations. In multi-agent performative prediction, with joint profile $\theta=(\theta^1,\dots,\theta^n)$ and joint-induced distribution $\mathcal D(\cdot\mid\theta)$, a profile $\theta^\ast\in\Theta^n$ is performatively stable if for every agent $i$,
$$
\theta^{\ast i}\in\arg\min_{\phi\in\Theta}\ell(\theta^\ast;\phi),
$$
so no single agent can unilaterally improve its own loss on the distribution induced by the full profile [2201.10483]. In reinforcement learning, a performatively stable policy is defined analogously as a policy that is optimal in the MDP induced by itself [2207.00046]. In control, a performatively stable control (PSC) solution is a fixed point of a distribution-shifted cost-minimizer mapping [2410.23251]. These formulations differ in state space and objective, but preserve the same equilibrium template.

## 2. Existence, uniqueness, and contraction regimes

The classical positive theory is based on contraction. In the foundational performative-prediction model, if $\ell$ is $\beta$-jointly smooth and $\gamma$-strongly convex in $\theta$, and if the distribution map is $\epsilon$-sensitive in Wasserstein-1,
$$
W_1(D(\theta),D(\theta'))\le \epsilon\|\theta-\theta'\|_2,
$$
then the retraining map satisfies
$$
\|G(\theta)-G(\theta')\|_2\le \frac{\epsilon\beta}{\gamma}\|\theta-\theta'\|_2.
$$
Hence, when $\epsilon<\gamma/\beta$, $G$ is a contraction, the performatively stable point is unique, and repeated risk minimization converges linearly to it [2002.06673].

A related stateful generalization replaces the stateless map $D(\theta)$ by a transition operator
$$
d_t=Tr(d_{t-1},\theta_t),
$$
so the population state itself carries memory. There, a stable pair $(d_S,\theta_S)$ satisfies both $Tr(d_S,\theta_S)=d_S$ and $\theta_S=G(d_S)$. If the transition map is $\epsilon$-jointly sensitive and the loss is $\gamma$-strongly convex and $\beta$-smooth, then the combined update is a contraction with coefficient
$$
\kappa=\epsilon\Bigl(1+\frac{\beta}{\gamma}\Bigr),
$$
and uniqueness plus linear convergence follow when $\kappa<1$ [2011.03885].

Later work relaxes the object with respect to which sensitivity is measured. For neural networks, the key assumption is Lipschitz sensitivity of the induced distribution with respect to the model’s predictions rather than its parameters:
$$
\chi^2(\mathcal D(f'),\mathcal D(f))\le L_p\|f-f'\|_{L^2(p)}^2.
$$
Combined with a bounded density-ratio condition and strong convexity of the loss in its prediction argument, this yields contraction of repeated risk minimization in function space. The contraction factor is
$$
\kappa=\frac{\sqrt{C\,L_p}\,M}{\gamma},
$$
and if $\kappa<1$ there is linear convergence to a unique performatively stable classifier. For nonconvex function classes, an approximation version yields convergence to an $O(\kappa)$-neighborhood under squared-error loss [2304.06879]. This shifts the regularity burden from parameterization to predictive behavior.

A different line of work shows that contraction is sufficient but not necessary for stability-like guarantees. In online performative prediction, if an online algorithm achieves sublinear external regret on the adaptive loss sequence $\ell_t(\theta)=\ell(z_t;\theta)$ with $z_t\sim\mathcal D(\theta_t)$, then the uniform distribution over its iterates is $(Reg(T)/T)$-performatively stable:
$$
\mathbb E_{\theta\sim\mu}\mathbb E_{z\sim\mathcal D(\theta)}[\ell(z;\theta)]
\le
\min_{\theta'}\mathbb E_{\theta\sim\mu}\mathbb E_{z\sim\mathcal D(\theta)}[\ell(z;\theta')]
+\frac{Reg(T)}{T}.
$$
This reduction requires no convexity, no smoothness, and no continuity assumptions on $\mathcal D(\cdot)$; it replaces deterministic fixed points by mixed equilibria [2602.24207].

## 3. Algorithms for reaching stable points

The basic computational procedures are repeated risk minimization (RRM) and repeated gradient descent (RGD). RRM iterates
$$
\theta_{t+1}\leftarrow \arg\min_{\theta\in\Theta}\mathbb E_{Z\sim D(\theta_t)}[\ell(Z;\theta)],
$$
whereas RGD performs one gradient step on the risk induced by the currently deployed model,
$$
\theta_{t+1}
=\theta_t-\eta\,\mathbb E_{Z\sim D(\theta_t)}[\nabla_\theta\ell(Z;\theta_t)].
$$
Under the standard contraction assumptions, both methods converge linearly to the same unique stable point [2002.06673].

For nonconvex losses, exact best-response maps are generally unavailable or multivalued, so the relevant target becomes a stationary performative stable (SPS) solution. In one formulation, $\theta^\star$ is a $\delta$-SPS solution if
$$
\|\nabla J(\theta^\star;\theta^\star)\|^2\le \delta,
$$
where
$$
J(\theta_1;\theta_2)=\mathbb E_{Z\sim\mathcal D(\theta_2)}[\ell(\theta_1;Z)].
$$
Greedy-deployment SGD,
$$
\theta_{t+1}=\theta_t-\gamma_{t+1}\nabla\ell(\theta_t;Z_{t+1}),\qquad Z_{t+1}\sim\mathcal D(\theta_t),
$$
converges in expectation to a biased SPS solution under smoothness, controlled gradient noise, and either Wasserstein-1 or total-variation sensitivity conditions. With constant-scale tuning, the asymptotic bias is
$$
\mathcal O(\sigma_0\epsilon+(1+\sigma_1^2)\epsilon^2),
$$
and a lazy deployment scheme reduces the dominant sensitivity-driven bias to $\mathcal O(\epsilon^2)$ [2405.17922].

Variance reduction improves this picture. The SPRINT algorithm adapts the SVRG idea to performative prediction by combining a snapshot gradient
$$
g^s=\nabla_\theta \mathbb E_{z\sim D(\tilde\theta^s)}[\ell(\tilde\theta^s;z)]
$$
with inner-loop control variates
$$
v_k^{s+1}
=
\nabla\ell(\theta_k^{s+1};z_{i_k})
-
\nabla\ell(\tilde\theta^s;z_{i_k})
+
g^s.
$$
Under $L$-smoothness in $\theta$, $\epsilon$-sensitivity of the distribution map, and pointwise Lipschitzness in $z$, SPRINT attains an $\mathcal O(1/T)$ rate to an SPS solution,
$$
\frac1T\sum_{t=0}^{T-1}\mathbb E[\|\nabla R(\theta_t)\|^2]
\le \frac{\Delta_0}{T\Gamma}+\Delta_1,
$$
with
$$
\Delta_1=\mathcal O(L_0^2\epsilon^2+L_0^4\epsilon^4),
$$
and this residual neighborhood is independent of stochastic-gradient variance [2509.17304].

The online-learning perspective yields a conceptually different algorithmic guarantee. Any no-regret method induces mixed performative stability through time-averaging. This includes standard online gradient descent with $\eta_t=O(1/\sqrt t)$, whose iterate mixture is $O(GD/\sqrt T)$-stable for convex losses, and stronger $O((1+\log T)/T)$ stability under strong convexity [2602.24207]. A plausible implication is that time-averaging is not merely an analysis device but a structural way to stabilize systems in which deterministic retraining maps may be discontinuous.

## 4. Multi-agent, dynamical, and structured extensions

In multi-agent performative prediction, multiple decision makers simultaneously influence a common outcome. The decoupled loss of agent $i$ is
$$
\ell(\theta;\theta^i)=\mathbb E_{(x,y)\sim\mathcal D(\cdot\mid\theta)}[\ell_{\mathrm{data}}((x,y),\theta^i)],
$$
and performative stability becomes a Nash-like fixed point of unilateral retraining [2201.10483]. In a location-scale regression model,
$$
y=\langle \theta^0-\sum_{i=1}^n\lambda_i\theta^i,x\rangle+\text{noise},
$$
small exponentiated-gradient learning rates satisfying thresholds $\eta_\ast$ and $R_\eta$ imply global convergence from any interior start. The stable profile is also the unique global minimizer of the total loss $\sum_i\ell(\theta;\theta^i)$, so in that regime stability implies optimality [2201.10483].

Performative control extends the same idea to policy-dependent dynamical systems. With disturbance-action policy parameter $M$, cost
$$
C_T(M;M')
=
\mathbb E\Big[\sum_{t=0}^T c_t(x_t^{(M)},u_t^{(M)})\Big],
$$
and policy-dependent disturbance laws $\mathcal D_t(M')$, a PSC solution is a fixed point
$$
M^{PS}=\Phi(M^{PS}),\qquad
\Phi(M')=\arg\min_{M\in\mathbb M}C_T(M;M').
$$
Existence and uniqueness follow if a weighted sensitivity-propagation condition holds:
$$
\sum_{t=0}^{T-1}\Bigl(\varepsilon_t\sum_{i=t+1}^{T}\nu_i\Bigr)<\widetilde\mu.
$$
For almost-sure strongly stable dynamics this reduces to a small-sensitivity sum condition, whereas for almost-sure unstable dynamics it requires temporally backward decay of the sensitivities. A repeated stochastic gradient descent scheme then converges to the PSC solution with a non-asymptotic bound whose second term is $O(1/N)$ [2410.23251].

Performative reinforcement learning replaces decision-dependent supervised data with policy-dependent rewards and transition kernels. In one formulation, a performatively stable policy $\pi_S$ satisfies
$$
\pi_S\in\arg\max_{\pi'}V^{\pi'}_{\pi_S}(\rho),
$$
meaning that once the environment freezes at $M(\pi_S)$, no alternative policy performs better there. Regularized repeated optimization of the occupancy-measure program converges under Lipschitz dependence of rewards and transitions on the occupancy measure and sufficiently large regularization $\lambda$; projected-gradient and finite-sample versions also converge [2207.00046]. Subsequent work in performative RL separates stability from optimality and develops a performative policy gradient method for the latter, which underscores that the two solution concepts remain distinct in sequential decision problems [2512.20576].

Several later frameworks preserve the fixed-point idea while changing the ambient optimization problem. In decision-dependent chance-constrained optimization, a performative solution is a self-feasible optimizer $x^\ast\in\Phi(x^\ast)$, existence follows from Kakutani’s fixed-point theorem, and a stochastic scenario iteration with
$$
N_t=C\,(d+\log t)
$$
samples converges almost surely to the unique fixed point under a deterministic contraction condition $K<1$ [2603.29982]. In federated bilevel performative prediction, the federated bilevel performatively stable (FBPS) point is the fixed point of a decoupled upper-level response map, and a contraction condition of the form
$$
\frac{\bar L_{F_c}^x(\varepsilon_c,\varepsilon_d)}{\gamma_f}\le 1
$$
implies uniqueness and linear convergence of federated repeated-risk minimization, while a stochastic federated hypergradient method converges at rate $\mathcal O(1/r)$ in mean square under diminishing step sizes and sufficiently small sensitivities [2606.19734].

## 5. Instability, phase transitions, and computational hardness

Performative stability is not universally benign. The multi-agent theory exhibits a phase transition from convergence to instability and then to formal chaos. In the two-dimensional location-scale setting with common learning rate $\eta>0$ and total influence $L_n=\sum_i\lambda_i$ sufficiently large, the exponentiated-gradient dynamics become Li–Yorke chaotic. Under symmetric initialization, the system reduces to the one-dimensional map
$$
x_{t+1}=f_{u,v}(x_t)
=
\frac{x_t}{x_t+(1-x_t)e^{u(x_t-v)}},
$$
and existence of a period-3 orbit implies Li–Yorke chaos and periodic orbits of all periods. Simulations show monotone convergence for settings such as $(L,\eta)=(14,0.001)$ or $(1.4,0.05)$, but aperiodic sensitive oscillations for $(14,0.05)$ [2201.10483].

A separate line of work establishes a sharp complexity transition at the classical contraction threshold
$$
\rho=L\beta/\alpha.
$$
When $\rho<1$, repeated risk minimization is a contraction and converges linearly to the unique stable point. At $\rho=1$, even simple one-dimensional quadratic examples can cycle. For any sufficiently small $\epsilon>0$, computing an $\epsilon$-performatively stable point is PPAD-complete even when
$$
\rho\le 1+O(\epsilon),
$$
and this hardness persists for quadratic loss with linear distribution shifts,
$$
\ell(\theta;z)=\tfrac12\|\theta-z\|_2^2,\qquad D(\theta)=\delta_{A\theta+b}.
$$
The hardness also extends from hypercubes to general well-bounded convex domains. In strategic classification, computing a strategic local optimum is PLS-hard [2601.20180].

These negative results clarify a common misunderstanding. Stability is often introduced through contraction-based sufficient conditions, but those conditions do not characterize the full landscape. Outside the weak-performativity regime, deterministic stable points may fail to exist, may be hard to find, or may be dynamically unstable. The mixed-equilibrium reduction via no-regret learning provides one route around this barrier by replacing deterministic fixed points with distributions over iterates [2602.24207].

## 6. Fairness, nonstationarity, and statistical inference

Performative stability can conflict with social desiderata. In group-structured performative prediction with demographic groups $s\in\mathcal S$, the literature defines loss disparity
$$
\Delta_{\mathcal L}^{(t)}
=
\Bigl|
\mathbb E_{z\sim\mathcal D_a^{(t)}}[\ell(\theta^{(t)};z)]
-
\mathbb E_{z\sim\mathcal D_b^{(t)}}[\ell(\theta^{(t)};z)]
\Bigr|
$$
and participation disparity
$$
\Delta_p^{(t)}=|p_a^{(t)}-p_b^{(t)}|.
$$
Examples show that repeated risk minimization can converge to a performatively stable point whose induced distribution is highly polarized, including cases where one group vanishes and group-wise loss disparity becomes large [2406.16756].

Standard static-data fairness interventions may also break stability. Group-loss-variance regularization can render the per-round objective non-convex, and repeated DRO with $\chi^2$ constraints can generate oscillation with no fixed point. To address this, Fair-RRM replaces the standard per-round objective with one of three modified objectives: a group-level penalty, a sample-level penalty, or loss-guided re-weighting. Under the same strong convexity and joint-smoothness assumptions as standard performative prediction, but with modified smoothness parameter $\tilde\beta$, uniqueness and convergence are recovered when
$$
\epsilon(1+\tilde\beta/\gamma)<1.
$$
In the retention setting, larger fairness weight $\rho$ yields smaller group-loss disparity at the Fair-PS point [2406.16756].

Performative stability has also been extended to environments with both endogenous and exogenous drift. In partially performative prediction, the round-$t$ data law is
$$
D_t(\theta_t)=(1-\alpha_t)\,D(\theta_t)+\alpha_t\,P_t,
$$
and the relevant comparator is the time-varying stable point $\theta_t^{PS}$. The corresponding online criterion is stability regret,
$$
Reg_T^{PS}(A)=\sum_{t=1}^T\bigl(\mathbb E[PR_t(\theta_t)]-PR_t(\theta_t^{PS})\bigr).
$$
Under strong convexity, smoothness, sensitivity, and a path-length measure
$$
\Delta_t^{PS}=\|\theta_t^{PS}-\theta_{t+1}^{PS}\|,
$$
repeated risk minimization satisfies
$$
Reg_T^{PS}(RRM)
\le
\frac{L}{1-\gamma}\Bigl[\|\theta_1-\theta_1^{PS}\|+\sum_{t=1}^{T-1}\Delta_t^{PS}\Bigr],
$$
while repeated gradient descent achieves
$$
Reg_T^{PS}(RGD)=O\Bigl(1+\sum_{t=1}^{T-1}\Delta_t^{PS}\Bigr).
$$
This recasts stability as a tracking problem rather than convergence to a single fixed point [2606.07890].

Finally, recent work develops an inferential theory for performative stability itself. In a unified single- and multi-player setting, repeated risk minimization produces estimators $\hat\theta_t$ satisfying
$$
\sqrt N\,(\hat\theta_t-\theta_t)\Rightarrow N(0,\Sigma_t),
$$
with a recursive closed-form covariance formula involving the Jacobian of the solution map and the gradient covariance under the induced distribution. The same work shows that the estimator attains the semiparametric efficiency bound under the stated regularity conditions [2602.03049]. This shifts performative stability from a purely optimization-theoretic target to an estimable statistical object with uncertainty quantification.

Performative stability therefore occupies a broad middle ground between optimization, equilibrium analysis, and dynamical systems. In weakly performative regimes it is a contractive fixed point reachable by straightforward retraining; in richer settings it becomes a mixed equilibrium, a moving target, or a stability notion for control, reinforcement learning, bilevel optimization, and chance-constrained decision-making. The modern literature shows both its utility and its limits: stability can imply exact convergence and, in some regimes, optimality, but it can also coexist with polarization, undergo transitions to chaos, or become computationally intractable once feedback strength crosses the contraction threshold [2002.06673].

Source: https://www.emergentmind.com/topics/performative-stability