---
title: Multiplicative Selection Model
url: https://www.emergentmind.com/topics/multiplicative-selection-model
type: topic
---

# Multiplicative Selection Model

Searching arXiv for the supplied papers to ground the article and verify the relevant usages of “multiplicative selection model.”
arxiv_search(query="1208.3160", max_results=5, sort_by="relevance")
arxiv_search(query="1208.3160", max_results=5)
Searching for arXiv:1208.3160.
{"query":"1208.3160","max_results":5}
A multiplicative selection model is a family of formal constructions in which selection, attrition, or response acts through multiplicative reweighting of a current state. In population genetics, the canonical instance is the weak-selection haploid model, where allele frequencies evolve by the discrete replicator equation and coincide exactly with multiplicative-weights updates in a coordination game between genes [1208.3160, 1408.6270]. Closely related multiplicative formulations also appear in rank-ordered survival modeling [1703.10548], in missing-not-at-random identification with instrumental variables [2509.22499], and in statistical selection procedures based on multiplicative weights or multiplicative priors [2002.08663, 1903.03531].

## 1. Terminological scope

The term is not reserved for a single universally standardized model. In the evolutionary literature, it refers to fitness-proportional updating of allele frequencies under haploid selection. In sociophysics and electoral modeling, it refers to a multiplicative-decrease process generating a rank-ordered survival function. In missing-data theory, it refers to a multiplicative restriction on the nonresponse mechanism. A plausible unifying description is that the state variable is transformed by proportional factors rather than by additive increments.

| Domain | Core multiplicative form | Representative source |
|---|---|---|
| Weak-selection haploid evolution | Allele-frequency update by multiplicative fitness weighting | [1208.3160], [1408.6270] |
| Rank-ordered survival | \(F(i+1,s+1)=(1-\mu(i,s))F(i,s)\) | [1703.10548] |
| MNAR missing data with IVs | \(P\{R=0\mid Z,U,X\}=\exp\{\alpha_z+\alpha_u\}\) | [2509.22499] |
| Graphical-model learning and sparsity selection | Multiplicative-weights or multiplicative-prior updates | [2002.08663], [1903.03531] |

This multiplicative viewpoint matters because it links disparate areas to common analytical machinery: discrete replicator dynamics, regret bounds, transport equations, influence-function theory, and hierarchical Bayesian sparsity selection.

## 2. Weak selection in haploid population genetics

In "Multiplicative Updates in Coordination Games and the Theory of Evolution" [1208.3160], the model is specified for two genes with allele-frequency vectors
\[
x(t)=(x_1(t),\dots,x_m(t))^T,\qquad y(t)=(y_1(t),\dots,y_n(t))^T,
\]
and a fitness-landscape matrix \(W=(w_{ij})\) with \(w_{ij}\in[1-s,1+s]\) for small \(s>0\). Writing
\[
\Delta_{ij}:=\frac{w_{ij}-1}{s},
\]
the weak-selection dynamics on the Wright manifold are
\[
x_i(t+1)=x_i(t)\cdot \frac{1+s\cdot q_A(i;t)}{Z_x(t)},\qquad
y_j(t+1)=y_j(t)\cdot \frac{1+s\cdot q_B(j;t)}{Z_y(t)},
\]
where
\[
q_A(i;t):=\sum_{j=1}^n y_j(t)\Delta_{ij},\qquad
q_B(j;t):=\sum_{i=1}^m x_i(t)\Delta_{ij},
\]
and
\[
Z_x(t)=\sum_{i=1}^m x_i(t)[1+s\cdot q_A(i;t)],\qquad
Z_y(t)=\sum_{j=1}^n y_j(t)[1+s\cdot q_B(j;t)].
\]
Since \(\sum_i x_i=1\), one can also write
\[
Z_x(t)=1+s\cdot \bar q(t),\qquad
\bar q(t)=\sum_{i,j}x_i(t)y_j(t)\Delta_{ij}.
\]

This is exactly the weak-selection form of the discrete replicator equation. If \(f_i(t)=\sum_j y_j(t)w_{ij}\) and \(\phi(t)=\sum_{i,j}x_i y_j w_{ij}\), then
\[
x_i(t+1)=x_i(t)\cdot \frac{f_i(t)}{\phi(t)}
=x_i(t)\cdot \frac{1+s\cdot q_A(i;t)}{1+s\cdot \bar q(t)}.
\]
The biological interpretation is that alleles are boosted in proportion to their expected fitness against the current allelic distribution at the other locus.

The 2014 analysis of haploid selection uses the same multiplicative structure in the standard single-locus form
\[
p_i(t+1)=\frac{w_i\,p_i(t)}{\sum_j w_j\,p_j(t)},
\]
and then lifts it to a two-locus coordination-game representation. There, natural selection is described as applying the discrete multiplicative-weights update separately to each gene [1408.6270].

## 3. Coordination games, replicator dynamics, and mixability

The game-theoretic reformulation in [1208.3160] treats gene \(A\) and gene \(B\) as the two players of a coordination game. Gene \(A\) chooses an allele \(i\in\{1,\dots,m\}\), gene \(B\) chooses an allele \(j\in\{1,\dots,n\}\), and both players receive payoff \(\Delta_{ij}\). The mixed strategies are precisely the allele-frequency vectors \(x\) and \(y\). The expected payoff to pure strategy \(i\) is
\[
q_A(i)=\sum_j y_j\Delta_{ij},
\]
and symmetrically
\[
q_B(j)=\sum_i x_i\Delta_{ij}.
\]
Under multiplicative weights with step size \(\epsilon=s\),
\[
x_i\leftarrow x_i\cdot \frac{1+\epsilon\cdot q_A(i)}{1+\epsilon\cdot \sum_k x_k q_A(k)},
\]
with an analogous update for \(y\). Thus natural selection under weak selection is tantamount to multiplicative updates in a coordination game.

A central quantity in this formulation is mixability. The mixability of allele \(i\) of gene \(A\) at time \(t\) is
\[
m_x^t(i):=\sum_j y_j(t)\Delta_{ij},
\]
and that of allele \(j\) of gene \(B\) is
\[
m_y^t(j):=\sum_i x_i(t)\Delta_{ij}.
\]
These utilities are exactly the terms that drive the multiplicative updates. In this sense, the model gives a rigorous realization of the idea that natural selection under sex favors alleles with high average fitness against varied genetic backgrounds [1208.3160].

The 2014 coordination-game treatment makes the same identification in the discrete replicator form
\[
x_i' = x_i\cdot \frac{(A y)_i}{x^T A y},\qquad
y_j' = y_j\cdot \frac{(A^T x)_j}{x^T A y},
\]
where \(A\) is the common-payoff matrix of a symmetric coordination game. When \(A\) is identified with the genotype-fitness matrix \(B\), the term \((B y)_i\) is the average fitness, or mixability, of allele \(i\) at gene 1 [1408.6270].

## 4. Fixed points, support size, and the fate of diversity

Two major results define the long-run theory. First, [1208.3160] shows that under weak selection the discrete dynamics remain within \(O(s)\) of the Wright manifold and converge. Equivalently, the multiplicative-weights regret bound implies that in the limit \(T\to\infty\), the empirical payoffs of the evolving \(x(t),y(t)\) do no worse than the best fixed pure pair \((i^*,j^*)\) in hindsight. In particular, the population focuses on alleles of maximal cumulative mixability.

Second, [1408.6270] proves a pointwise convergence theorem for coordination games with all row and column entries distinct: the replicator map is a diffeomorphism on the interior of \(\Delta_n\times\Delta_n\), every orbit converges pointwise to a fixed point, and, except for an initial set of Lebesgue measure zero, every trajectory converges to a pure Nash equilibrium. The proof uses the Losert–Akin convergence theorem, a linearization argument showing that stable fixed points are weakly stable Nash equilibria, and the Center–Stable-Manifold Theorem to show that unstable fixed points have measure-zero basins of attraction.

A common misconception is that the existence of many mixed equilibria by itself implies long-run maintenance of polymorphism. The two analyses separate existence from attraction. At a fixed point, all alleles in the support of \(x\) must have the same expected payoff, and likewise for \(y\). If \(\operatorname{supp}(x)\) and \(\operatorname{supp}(y)\) have size \(k\), then necessarily both supports have size \(k\), and the corresponding \(k\times k\) submatrix \(A\) of \(W\) satisfies
\[
A x = a\cdot 1,\qquad A^T y = a\cdot 1,
\]
for some scalar \(a>1\) and probability vectors \(x,y>0\). Writing \(A=U+sB\), one obtains
\[
B x=(a-1)\cdot 1,\qquad B^T y=(a-1)\cdot 1.
\]
This yields explicit conditions for nontrivial mixed equilibria and shows that equilibria in two-person coordination games can have large supports [1208.3160].

The same paper derives a probabilistic lower bound: if the \(\Delta_{ij}\) are iid symmetric continuous on \([-1,1]\), then a random \(k\times k\) matrix \(B\) has \(B^{-1}\) with all row- and column-sums positive with probability at least \(2^{-(2k-1)}\). Hence any \(k\times k\) submatrix \(A=U+sB\) has probability at least \(2^{-(2k-1)}\) of supporting a nontrivial equilibrium, and the expected number of such equilibria is at least
\[
{m\choose k}{n\choose k}\cdot 2^{-(2k-1)}
\ge 2\cdot \left(\frac{mn}{4k^2}\right)^k.
\]
This suggests a two-level picture: mixed equilibria can be combinatorially abundant, but under the genericity assumptions of [1408.6270] they do not attract a positive-measure set of initial conditions. Biologically, natural selection alone therefore leads almost surely to fixation, so long-term preservation of genetic diversity requires additional mechanisms such as mutation, recombination, or speciation [1408.6270].

## 5. Rank-ordered survival and beta-like selection

In "A multiplicative process for generating a beta-like survival function with application to the UK 2016 EU referendum results" [1703.10548], the multiplicative selection idea is recast as a generative model for rank-ordered vote shares. Let \(F(i,s)\in[0,1]\) be the expected survival in district \(i\) after processing up to stage \(s\), and let \(\mu(i,s)\in[0,1]\) be an attrition or mortality function. The discrete dynamics are
\[
F(0,s)=1,\qquad
F(i+1,s+1)=(1-\mu(i,s))F(i,s),
\]
with boundary condition \(F(i,0)=0\) for \(i>0\). In a continuous approximation,
\[
\frac{\partial f}{\partial s}+\frac{\partial f}{\partial i}+\mu(i,s)f(i,s)=0,
\]
and the solution is
\[
f(i,s)=\exp\!\Bigl(-\!\!\int_0^i \mu(i-t,s-t)\,dt\Bigr).
\]

For the mixture attrition
\[
\mu(i,s)=\frac{\alpha}{i+\kappa}+\frac{\beta}{s},
\qquad \alpha,\beta>0,\ \kappa>0,
\]
the solution becomes
\[
f(i,s)=\Bigl(\frac{\kappa}{i+\kappa}\Bigr)^{\alpha}\Bigl(1-\frac{i}{s}\Bigr)^{\beta}.
\]
Fixing \(\kappa=0.5\) and introducing a scale constant \(C\) yields the beta-like survival function
\[
S(i)=C\;\Bigl(\frac{0.5}{i+0.5}\Bigr)^{\alpha}\Bigl(1-\frac{i}{s}\Bigr)^{\beta}.
\]
For large \(i\), the asymptotic form is
\[
S(i)\approx K\,i^{-\alpha}(s-i)^{\beta},
\qquad K=\kappa^{\alpha}s^{-\beta}.
\]

The parameters have distinct shape roles: \(\alpha\) controls the intermediate-rank power-law decay, \(\beta\) controls the cutoff near \(i\to s\), and \(C\) scales the curve. Estimation is performed by nonlinear least squares or direct maximum likelihood using
\[
\mathrm{SSE}(\alpha,\beta,C)
=
\sum_{i=0}^{n}\bigl[V_i-C(0.5/(i+0.5))^\alpha(1-i/s)^\beta\bigr]^2.
\]
Applied to the UK result with \(n=382\) Local Authorities, the nonlinear fit for rank-ordered Leave shares gave \(\hat\alpha\approx0.0357\), \(\hat\beta\approx0.2094\), \(\hat C\approx0.7801\), and \(R^2\approx0.9913\); for Remain shares it gave \(\hat\alpha\approx0.1286\), \(\hat\beta\approx0.1244\), \(\hat C\approx1.0740\), and \(R^2\approx0.9930\). The paper also rank-orders census covariates and reports that this method outperformed simple linear regression on raw census percentages, especially for covariates with strong geographic clustering [1703.10548].

## 6. Missing-data identification and related selection frameworks

In missing-data theory, the multiplicative selection model appears as a restriction on the nonresponse mechanism rather than on a dynamical update. In "A Multiplicative Instrumental Variable Model for Data Missing Not-at-Random" [2509.22499], each subject contributes outcome \(Y\), missingness indicator \(R\), covariates \(X\), instrument \(Z\), and unobserved factor \(U\). The core multiplicative selection assumption is
\[
P\{R=0\mid Z=z,U=u,X=x\}
=
\exp\{\alpha_z(z,x)+\alpha_u(u,x)\},
\]
combined with \(Y\perp (Z,R)\mid (U,X)\) and \(U\perp Z\mid X\). This excludes interaction between \(Z\) and \(U\) on the log scale and leaves the degree of selection bias on the outcome scale unrestricted. Identification proceeds through a “single-arm Wald ratio”: if
\[
\pi_z(x)=P(R=1\mid Z=z,X=x),\qquad
\mu_z(x)=E[R\,h(Y;\psi_0)\mid Z=z,X=x],
\]
and
\[
\delta^R(x)=\pi_1(x)-\pi_0(x),\qquad
\delta^Y(x)=\mu_1(x)-\mu_0(x),\qquad
\delta(x)=\frac{\delta^Y(x)}{\delta^R(x)},
\]
then
\[
E[h(Y;\psi_0)\mid R=0]=E[\delta(X)\mid R=0].
\]
The paper derives the efficient influence function, develops semiparametric multiply robust IV estimators, extends the framework to polytomous and continuous instruments, and reports that in a Botswana HIV survey application the IF-based MIV estimate of HIV prevalence among nonresponders was \(38.4\%\) \((95\%\ \mathrm{CI}\ 33.1\text{–}43.7)\), versus \(21.4\%\) among responders, leading to an adjusted overall estimate of \(24.7\%\) \((95\%\ \mathrm{CI}\ 18.6\text{–}30.8)\) [2509.22499].

Related multiplicative formulations also appear in statistical selection and structure learning. For Gaussian graphical models, "Learning Gaussian Graphical Models via Multiplicative Weights" uses a Sparsitron update in which a positive weight vector is updated coordinatewise by
\[
v^{t+1}_j=v^t_j\cdot \beta^{\ell^t_j},
\]
with per-node runtime \(O(mp)\), total runtime \(O(mp^2)\), and sample complexity
\[
m=T+M=O\!\left(\lambda^4\nu_{\max}^2\theta_{\max}^2\kappa^{-4}\cdot \ln^3(p/\delta)\right)
\]
under the stated edge-strength, \(\ell_1\)-norm, variance, and diagonal bounds [2002.08663]. For high-dimensional Gaussian DAG models, "Consistent Bayesian Sparsity Selection for High-dimensional Gaussian DAG Models with Multiplicative and Beta-mixture Priors" introduces a multiplicative prior in which node-specific latent weights \(\omega_j\sim \mathrm{Beta}(\alpha_1,\alpha_2)\) determine edge sparsity in the Cholesky factor of the precision matrix, and proves posterior-mode consistency under assumptions A1–A6; in the reported simulation table at \(3\%\) sparsity, SSC–Multiplicative attained PPV \(0.83\), TPR \(0.87\), and MCC \(0.84\) [1903.03531].

These later uses are not identical to the haploid evolutionary model, but they preserve the same structural idea: selection or recovery is implemented through multiplicative reweighting, multiplicative factorization, or multiplicative priors. Across domains, that structure makes the models analytically tractable and connects them to established methods in online learning, dynamical systems, semiparametric inference, and Bayesian high-dimensional statistics.

Source: https://www.emergentmind.com/topics/multiplicative-selection-model