Random selective advantages of a gene in a finite population

1969 ◽  
Vol 6 (1) ◽  
pp. 19-37 ◽  
Author(s):  
Louis Jensen ◽  
Edward Pollak

A problem of interest to many population geneticists is the process of change in a gene frequency. A popular model used to describe the change in a gene frequency involves the assumption that the gene frequency is Markovian. The probabilities in a Markov process can be approximated by the solution of a partial differential equation known as the Fokker-Planck equation or the forward Kolmogorov equation. Mathematically this equation is where subscripts indicate partial differentiation. In this equation, f(p, x; t) is the probability density that the frequency of a gene is x at time t, given that the frequency was p at time t = o. The expressions MΔX and VΔx are, respectively, the first and second moments of the change in the gene frequency during one generation. A rigorous derivation of this equation was given by Kolmogorov (1931).

1969 ◽  
Vol 6 (01) ◽  
pp. 19-37 ◽  
Author(s):  
Louis Jensen ◽  
Edward Pollak

A problem of interest to many population geneticists is the process of change in a gene frequency. A popular model used to describe the change in a gene frequency involves the assumption that the gene frequency is Markovian. The probabilities in a Markov process can be approximated by the solution of a partial differential equation known as the Fokker-Planck equation or the forward Kolmogorov equation. Mathematically this equation is where subscripts indicate partial differentiation. In this equation, f(p, x; t) is the probability density that the frequency of a gene is x at time t, given that the frequency was p at time t = o. The expressions MΔX and VΔx are, respectively, the first and second moments of the change in the gene frequency during one generation. A rigorous derivation of this equation was given by Kolmogorov (1931).


Author(s):  
J. R. Cannon ◽  
Yanping Lin ◽  
Shingmin Wang

AbstractThe authors consider in this paper the inverse problem of finding a pair of functions (u, p) such thatwhere F, f, E, s, αi, βi, γi, gi, i = 1, 2, are given functions.The existence and uniqueness of a smooth global solution pair (u, p) which depends continuously upon the data are demonstrated under certain assumptions on the data.


2020 ◽  
Vol 62 (1) ◽  
pp. 3-17 ◽  
Author(s):  
MATTHEW J. SIMPSON

We consider an extension of the classical Fisher–Kolmogorov equation, called the “Fisher–Stefan” model, which is a moving boundary problem on $0<x<L(t)$. A key property of the Fisher–Stefan model is the “spreading–vanishing dichotomy”, where solutions with $L(t)>L_{\text{c}}$ will eventually spread as $t\rightarrow \infty$, whereas solutions where $L(t)\ngtr L_{\text{c}}$ will vanish as $t\rightarrow \infty$. In one dimension it is well known that the critical length is $L_{\text{c}}=\unicode[STIX]{x1D70B}/2$. In this work, we re-formulate the Fisher–Stefan model in higher dimensions and calculate $L_{\text{c}}$ as a function of spatial dimensions in a radially symmetric coordinate system. Our results show how $L_{\text{c}}$ depends upon the dimension of the problem, and numerical solutions of the governing partial differential equation are consistent with our calculations.


1980 ◽  
Vol 87 (3) ◽  
pp. 515-521
Author(s):  
Albert E. Heins

In a recent paper, hereafter referred to as I (1) we derived two alternate forms for the fundamental solution of the axially symmetric wave equation. We demonstrated that for α > 0, the fundamental solution (the so-called free space Green's function) of the partial differential equationcould be written asif b > rorif r > b.


1966 ◽  
Vol 7 (3) ◽  
pp. 136-144 ◽  
Author(s):  
Ian N. Sneddon

The statical Reissner-Sagoci problem [1, 2, 3] is that of determining the components of stress and displacement in the interior of the semi-infinite homogeneous isotropic elastic solid z ≧ 0 when a circular area (0 ≦ p ≦ a, z = 0) of the boundary surface is forced to rotate through an angle a about an axis which is normal to the undeformed plane surface of the medium. It is easily shown that, if we use cylindrical coordinates (p, φ, z), the displacement vector has only one non-vanishing component uφ (p, z), and the stress tensor has only two non-vanishing components, σρπ(p, z) and σπz(p, z). The stress-strain relations reduce to the two simple equationswhere μ is the shear modulus of the material. From these equations, it follows immediately that the equilibrium equationis satisfied provided that the function uπ(ρ, z) is a solution of the partial differential equationThe boundary conditions can be written in the formwhere, in the case in which we are most interested, f(p) = αρ. We also assume that, as r → ∞, uπ, σρπ and σπz all tend to zero.


1922 ◽  
Vol 41 ◽  
pp. 76-81
Author(s):  
E. T. Copson

Riemann's method of solution of a linear second order partial differential equation of hyperbolic type was introduced in his memoir on sound waves. It has been used by Darboux in discussing the equationwhere α, β, γ are functions of x and y.


1959 ◽  
Vol 11 ◽  
pp. 148-155 ◽  
Author(s):  
Louis Weisner

On replacing the parameter n in Bessel's differential equation1.1by the operator y(∂/∂y), the partial differential equation Lu = 0 is constructed, where1.2This operator annuls u(x, y) = v(x)yn if, and only if, v(x) satisfies (1.1) and hence is a cylindrical function of order n. Thus every generating function of a set of cylindrical functions is a solution of Lu = 0.It is shown in § 2 that the partial differential equation Lu = 0 is invariant under a three-parameter Lie group. This group is then applied to the systematic determination of generating functions for Bessel functions, following the methods employed in two previous papers (4; 5).


1972 ◽  
Vol 15 (2) ◽  
pp. 229-234
Author(s):  
Julius A. Krantzberg

We consider the initial-boundary value problem for the parabolic partial differential equation1.1in the bounded domain D, contained in the upper half of the xy-plane, where a part of the x-axis lies on the boundary B(see Fig.1).


1974 ◽  
Vol 11 (1) ◽  
pp. 1-15 ◽  
Author(s):  
Louis Jensen

A technique for solving a partial differential equation is presented. The technique is based upon the known solution of a similar equation. The method is used to attempt to solve the equation describing the change in the frequency of an allele in the presence of selection and random drift in a finite population. Two cases can be solved within a reasonable degree of approximation: (a) the viabilities are additive and (b) heterozygotes are symmetrically overdominant to the homozygotes. The solutions in both cases are compared with the exact discrete solutions found by powering the re evant transition matrix.


1973 ◽  
Vol 73 (1) ◽  
pp. 157-165 ◽  
Author(s):  
R. E. Scraton

AbstractAn approximation is found to the solution of the partial differential equationin the region −1 ≤ x ≤ 1, t > 0, where u satisfies a general linear boundary condition on x = ± 1. This approximation is a polynomial in x, and is an exact solution of a perturbed form of the differential equation. By choosing the perturbation appropriately, this approach is mathematically equivalent to some recent methods for solving the differential equation in the form of a Chebyshev series. Better approximations to the required solution (and particularly to the eigenvalues) are obtained by choosing the perturbation to satisfy a least squares criterion.


Sign in / Sign up

Export Citation Format

Share Document