On random quadratic forms: supports of potential local maxima

Abstract

Abstract The selection model in population genetics is a dynamic system on the set of the probability distributions 𝒑=(p1,…,pn) of the alleles A1…,An, with pi(t+1) proportional to pi(t) multiplied by ∑jfi,jpj(t), and fi,j=fj,i interpreted as a fitness of the gene pair (Ai,Aj). It is known that 𝒑̂ is a locally stable equilibrium if and only if 𝒑̂ is a strict local maximum of the quadratic form 𝒑T𝒇𝒑. Usually, there are multiple local maxima and lim𝒑(t) depends on 𝒑(0). To address the question of a typical behavior of {𝒑(t)}, John Kingman considered the case when the fi,j are independent and [0,1]-uniform. He proved that with high probability (w.h.p.) no local maximum may have more than 2.49n1∕2 positive components, and reduced 2.49 to 2.14 for a nonbiological case of exponentials on [0,∞). We show that the constant 2.14 serves a broad class of smooth densities on [0,1] with the increasing hazard rate. As for a lower bound, we prove that w.h.p. for all k≤2n1∕3, there are many k-element subsets of [n] that pass a partial test to be a support of a local maximum. Still, it may well be that w.h.p. the actual supports are much smaller. In that direction, we prove that w.h.p. a support of a local maximum, which does not contain a support of a local equilibrium, is very unlikely to have size exceeding ⅔log2n and, for the uniform fitnesses, there are super-polynomially many potential supports free of local equilibriums of size close to ½log2n.