Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA hierarchical extension scheme for solutions of the Wright-Fisher model
87
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We develop a global and hierarchical scheme for the forward Kolmogorov (Fokker-Planck) equation of the diffusion approximation of the Wright-Fisher model of population genetics. That model describes the random genetic drift of several alleles at the same locus in a population. The key of our scheme is to connect the solutions before and after the loss of an allele. Whereas in an approach via stochastic processes or partial differential equations, such a loss of an allele leads to a boundary singularity, from a biological or geometric perspective, this is a natural process that can be analyzed in detail. Our method depends on evolution equations for the moments of the process and a careful analysis of the boundary flux.
rate research
Read More
We develop an iterative global solution scheme for the backward Kolmogorov equation of the diffusion approximation of the Wright-Fisher model of population genetics. That model describes the random genetic drift of several alleles at the same locus in a population from a backward perspective. The key of our scheme is to connect the solutions before and after the loss of an allele. Whereas in an approach via stochastic processes or partial differential equations, such a loss of an allele leads to a boundary singularity, from a biological or geometric perspective, this is a natural process that can be analyzed in detail. A clarification of the role of the boundary resolves certain uniqueness issues and enlucidates the construction of hierarchical solutions.
The diffusion approximation of the Wright-Fisher model of population genetics leads to partial differentiable equations, the so-called Kolmogorov equations, with an operator that degenerates at the boundary. Standard tools do not apply, and in fact, solutions lack regularity properties. In this paper, we develop a regularising blow-up scheme for a certain class of solutions of the backward Kolmogorov equation, the iteratively extended global solutions presented in cite{THJ5}, and establish their uniqueness. As the model describes the random genetic drift of several alleles at the same locus from a backward perspective, the singularities result from the loss of an allele. While in an analytical approach, this causes substantial difficulties, from a biological or geometric perspective, this is a natural process that can be analyzed in detail. The presented scheme regularises the solution via a tailored successive transformation of the domain.
We develop a general solution for the Fokker-Planck (Kolomogorov) equation representing the diffusion limit of the Wright-Fisher model of random genetic drift for an arbitrary number of alleles at a single locus. From this solution, we can readily deduce information about the evolution of a Wright-Fisher population.
We investigate the behaviour of the genealogy of a Wright-Fisher population model under the influence of a strong seed-bank effect. More precisely, we consider a simple seed-bank age distribution with two atoms, leading to either classical or long genealogical jumps (the latter modeling the effect of seed-dormancy). We assume that the length of these long jumps scales like a power $N^beta$ of the original population size $N$, thus giving rise to a `strong seed-bank effect. For a certain range of $beta$, we prove that the ancestral process of a sample of $n$ individuals converges under a non-classical time-scaling to Kingmans $n-$coalescent. Further, for a wider range of parameters, we analyze the time to the most recent common ancestor of two individuals analytically and by simulation.
We derive and apply a partial differential equation for the moment generating function of the Wright-Fisher model of population genetics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
