Fixation Maximization in the Positional Moran Process

Abstract

The Moran process is a classic stochastic process that models spread dynamics on graphs. A single "mutant" (e.g., a new opinion, strain, social trait etc.) invades a population of "residents" scattered over the nodes of a graph. The mutant fitness advantage δ ≥ 0 determines how aggressively mutants propagate to their neighbors. The quantity of interest is the fixation probability, i.e., the probability that the initial mutant eventually takes over the whole population. However, in realistic settings, the invading mutant has an advantage only in certain locations. E.g., the ability to metabolize a certain sugar is an advantageous trait to bacteria only when the sugar is actually present in their surroundings. In this paper, we introduce the positional Moran process, a natural generalization in which the mutant fitness advantage is only realized on specific nodes called active nodes, and study the problem of fixation maximization: given a budget k, choose a set of k active nodes that maximize the fixation probability of the invading mutant. We show that the problem is NP-hard, while the optimization function is not submodular, indicating strong computational hardness. We then focus on two natural limits: at the limit of δ → ∞ (strong selection), although the problem remains NP-hard, the optimization function becomes submodular and thus admits a constant-factor greedy approximation; at the limit of δ → 0 (weak selection), we show that we can obtain a tight approximation in O(m ω ) time, where m is the number of edges and ω is the matrix-multiplication exponent. An experimental evaluation of the new algorithms along with some proposed heuristics corroborates our results.

Fixation Maximization in the Positional Moran Process

“…In the rest of this section, we state Lemma 7 and give intuition about the function α(•). For the formal proof of Lemma 7, see the full version (Brendborg et al 2022).…”