在无限位突变模型下,用于时间均的凝聚过程的矩阵分析采样公式
Asger Hobolth1, Simon Boitard2, Andreas Futschik3
1Department of Mathematics, Aarhus University, Aarhus, Denmark.
Theoretical population biology
|April 3, 2025
概括
我们开发了一个计算框架,用凝聚过程和突变模型计算遗传样本概率. 这种方法是高效和稳定的,提供了对人口遗传学数据的洞察力.
科学领域:
- 人口遗传学 人口遗传学
- 计算生物学 计算生物学
- 进化生物学 进化生物学
背景情况:
- 凝聚理论是从DNA序列推断人口遗传史的基础.
- 无限位突变模型是分析遗传变异的标准假设.
- 在这些模型下计算样本概率对于统计推理至关重要.
研究的目的:
- 为计算遗传样本概率开发一个通用和计算效率高的框架.
- 提供适用于各种凝聚和突变模型的矩阵分析方法.
- 评估该框架对分析人口遗传数据和比较人口模型的有用性.
主要方法:
- 利用多变量相型理论来定义凝聚过程.
- 开发了一个基于速率矩阵,初始状态向量和奖励矩阵的概率生成函数.
- 实现了一个计算稳定的算法,涉及用于概率计算的矩阵运算.
主要成果:
- 推导出一种计算人口遗传数据集概率的一般方法.
- 对于少量突变来说,已经证明了计算吸引力.
- 展示了该方法对不同样本表示和人口模型 (例如结构化和Beta-coalescents) 的适用性.
结论:
- 开发的框架提供了一个计算稳定和高效的方式来计算样本概率.
- 塔吉马的D-统计是频谱概率的糟糕预测器,强调了需要更复杂的方法.
- 这项研究提供了对人口参数如何影响遗传谱概率的更深入的理解.
相关概念视频
Hardy-Weinberg Principle
71.4K
Diploid organisms have two alleles of each gene, one from each parent, in their somatic cells. Therefore, each individual contributes two alleles to the gene pool of the population. The gene pool of a population is the sum of every allele of all genes within that population and has some degree of variation. Genetic variation is typically expressed as a relative frequency, which is the percentage of the total population that has a given allele, genotype or phenotype.
71.4K
Mutation, Gene Flow, and Genetic Drift
57.7K
In a population that is not at Hardy-Weinberg equilibrium, the frequency of alleles changes over time. Therefore, any deviations from the five conditions of Hardy-Weinberg equilibrium can alter the genetic variation of a given population. Conditions that change the genetic variability of a population include mutations, natural selection, non-random mating, gene flow, and genetic drift (small population size).
57.7K
Mechanistic Models: Compartment Models in Individual and Population Analysis
21
Mechanistic models are utilized in individual analysis using single-source data, but imperfections arise due to data collection errors, preventing perfect prediction of observed data. The mathematical equation involves known values (Xi), observed concentrations (Ci), measurement errors (εi), model parameters (ϕj), and the related function (ƒi) for i number of values. Different least-squares metrics quantify differences between predicted and observed values. The ordinary least...
21
Noncompartmental Analysis: Statistical Moment Theory
65
Noncompartmental analyses leverage statistical moment theory to examine time-related changes in macroscopic events, encapsulating the collective outcomes stemming from the constituent elements in play. Statistical moment theory is a mathematical approach used to describe the time course of drug concentration in the body without assuming a specific compartmental model. SMT provides insights into drug absorption, distribution, metabolism, and elimination by treating drug concentration versus time...
65
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
35
Mechanistic models play a crucial role in algorithms for numerical problem-solving, particularly in nonlinear mixed effects modeling (NMEM). These models aim to minimize specific objective functions by evaluating various parameter estimates, leading to the development of systematic algorithms. In some cases, linearization techniques approximate the model using linear equations.
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
35


