Estimating and testing selection: the two-alleles, genic selection diffusion model

1979 ◽  
Vol 11 (1) ◽  
pp. 14-30 ◽  
Author(s):  
G. A. Watterson

The estimation, and testing for the presence, of a selective advantage of one allele over another is considered. It is assumed that a population's allele frequencies are known from some initial time until fixation of one or other allele occurs. The statistics needed to perform the estimation and testing are the heterozygosity of the population summed over all generations, and the observation of which allele fixes. It is shown that certain asymptotic probability distributions arise which are similar to those proved by Brown and Hewitt for statistical inference from diffusion processes, but their results assumed that the diffusion had a stationary density which is not the case for alleles which fix.The genetic diffusion may be transformed to Brownian motion with constant drift, and the inference questions concerning selection can be transformed to questions about the first exit of a Brownian motion from an interval. It is thus possible to construct significance tests, and to calculate the power of those tests, for detecting selection.

1979 ◽  
Vol 11 (01) ◽  
pp. 14-30 ◽  
Author(s):  
G. A. Watterson

The estimation, and testing for the presence, of a selective advantage of one allele over another is considered. It is assumed that a population's allele frequencies are known from some initial time until fixation of one or other allele occurs. The statistics needed to perform the estimation and testing are the heterozygosity of the population summed over all generations, and the observation of which allele fixes. It is shown that certain asymptotic probability distributions arise which are similar to those proved by Brown and Hewitt for statistical inference from diffusion processes, but their results assumed that the diffusion had a stationary density which is not the case for alleles which fix. The genetic diffusion may be transformed to Brownian motion with constant drift, and the inference questions concerning selection can be transformed to questions about the first exit of a Brownian motion from an interval. It is thus possible to construct significance tests, and to calculate the power of those tests, for detecting selection.


Author(s):  
Jeffrey J. Gory ◽  
Radu Herbei ◽  
Laura S. Kubatko

Abstract The increasing availability of population-level allele frequency data across one or more related populations necessitates the development of methods that can efficiently estimate population genetics parameters, such as the strength of selection acting on the population(s), from such data. Existing methods for this problem in the setting of the Wright-Fisher diffusion model are primarily likelihood-based, and rely on numerical approximation for likelihood computation and on bootstrapping for assessment of variability in the resulting estimates, requiring extensive computation. Recent work has provided a method for obtaining exact samples from general Wright-Fisher diffusion processes, enabling the development of methods for Bayesian estimation in this setting. We develop and implement a Bayesian method for estimating the strength of selection based on the Wright-Fisher diffusion for data sampled at a single time point. The method utilizes the latest algorithms for exact sampling to devise a Markov chain Monte Carlo procedure to draw samples from the joint posterior distribution of the selection coefficient and the allele frequencies. We demonstrate that when assumptions about the initial allele frequencies are accurate the method performs well for both simulated data and for an empirical data set on hypoxia in flies, where we find evidence for strong positive selection in a region of chromosome 2L previously identified. We discuss possible extensions of our method to the more general settings commonly encountered in practice, highlighting the advantages of Bayesian approaches to inference in this setting.


Author(s):  
ILYA V. TELYATNIKOV

We consider surface measures on the set of trajectories in a smooth compact Riemannian submanifold of Euclidean space generated by diffusion processes in the ambient space. A construction of surface measures on the path space of a smooth compact Riemannian submanifold of Euclidean space was introduced by Smolyanov and Weizsäcker for the case of the standard Brownian motion. The result presented in this paper extends the result of Smolyanov and Weizsäcker to the case when we consider measures generated by diffusion processes in the ambient space with nonidentical correlation operators. For every partition of the time interval, we consider the marginal distribution of the diffusion process in the ambient space under the condition that it visits the manifold at all times of the partition, when the mesh of the partition tends to zero. We prove the existence of some limit surface measures and the equivalence of the above measures to the distribution of some diffusion process on the manifold.


Author(s):  
Tianyu Ma ◽  
Vladimir S. Matveev ◽  
Ilya Pavlyukevich

AbstractWe show that geodesic random walks on a complete Finsler manifold of bounded geometry converge to a diffusion process which is, up to a drift, the Brownian motion corresponding to a Riemannian metric.


Entropy ◽  
2018 ◽  
Vol 20 (11) ◽  
pp. 813 ◽  
Author(s):  
José Amigó ◽  
Sámuel Balogh ◽  
Sergio Hernández

Entropy appears in many contexts (thermodynamics, statistical mechanics, information theory, measure-preserving dynamical systems, topological dynamics, etc.) as a measure of different properties (energy that cannot produce work, disorder, uncertainty, randomness, complexity, etc.). In this review, we focus on the so-called generalized entropies, which from a mathematical point of view are nonnegative functions defined on probability distributions that satisfy the first three Shannon–Khinchin axioms: continuity, maximality and expansibility. While these three axioms are expected to be satisfied by all macroscopic physical systems, the fourth axiom (separability or strong additivity) is in general violated by non-ergodic systems with long range forces, this having been the main reason for exploring weaker axiomatic settings. Currently, non-additive generalized entropies are being used also to study new phenomena in complex dynamics (multifractality), quantum systems (entanglement), soft sciences, and more. Besides going through the axiomatic framework, we review the characterization of generalized entropies via two scaling exponents introduced by Hanel and Thurner. In turn, the first of these exponents is related to the diffusion scaling exponent of diffusion processes, as we also discuss. Applications are addressed as the description of the main generalized entropies advances.


2014 ◽  
Vol 13 (04) ◽  
pp. 1430001 ◽  
Author(s):  
Jaume Masoliver

We review the level-crossing problem which includes the first-passage and escape problems as well as the theory of extreme values (the maximum, the minimum, the maximum absolute value and the range or span). We set the definitions and general results and apply them to one-dimensional diffusion processes with explicit results for the Brownian motion and the Ornstein–Uhlenbeck (OU) process.


Sign in / Sign up

Export Citation Format

Share Document