MCMC-CE:用于估计正方形形状的小右尾概率的新高效算法,在基因组学中具有应用
bioRxiv : the preprint server for biology
|April 1, 2025
概括
我们开发了MCMC-CE,这是一个新算法,用于准确计算大规模二次形式的小p值,对基因组研究至关重要. 这种方法提高了复杂的统计分析的计算效率和可靠性.
科学领域:
- 统计 统计 统计 统计
- 生物信息学是一种生物信息学.
- 基因组学就是基因组学.
背景情况:
- 多变量正常变量的二次形式在基因组学和生物信息学中是必不可少的.
- 计算大规模二次方形的小右尾概率 (p值) 是由于难以处理的分布和数值限制而具有计算挑战.
研究的目的:
- 提出一个创新的算法,MCMC-CE,以高效准确地估计大规模二次方形的小p值.
- 在基因组研究中解决与这些计算相关的计算挑战和数值约束.
主要方法:
- 马尔科夫链蒙特卡洛 (MCMC) 采样与交叉 (CE) 方法的整合.
- 纳入领先的自身价值提取和Satterthwaite类型的近似技术.
- 估计等级超过10,000的二次方形的p值.
主要成果:
- 与现有方法 (戴维斯,Imhof,Farebrother,-Tang-Zhang,点近似) 相比,MCMC-CE表现出优势的数值准确性和计算可靠性.
- 该算法高效地估计大等级二次数形式的小p值.
- 通过广泛的模拟和现实世界的基因组应用,包括GWAS和途径丰富分析来验证.
结论:
- MCMC-CE提供了一个强大的和可扩展的解决方案,用于在大型二次形式中准确的小p值计算.
- 在大型基因组研究中促进更精确的统计推断.
- 克服了生物信息学中复杂的统计问题的现有计算方法的局限性.
相关概念视频
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
34
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...
34
Distributions to Estimate Population Parameter
4.0K
The accurate values of population parameters such as population proportion, population mean, and population standard deviation (or variance) are usually unknown. These are fixed values that can only be estimated from the data collected from the samples. The estimates of each of these parameters are sample proportion, the sample mean, and sample standard deviation (or variance). To obtain the values of these sample statistics, data are required that have particular distribution and central...
4.0K
Parametric Survival Analysis: Weibull and Exponential Methods
295
Parametric survival analysis models survival data by assuming a specific probability distribution for the time until an event occurs. The Weibull and exponential distributions are two of the most commonly used methods in this context, due to their versatility and relatively straightforward application.
Weibull Distribution
The Weibull distribution is a flexible model used in parametric survival analysis. It can handle both increasing and decreasing hazard rates, depending on its shape parameter...
Weibull Distribution
The Weibull distribution is a flexible model used in parametric survival analysis. It can handle both increasing and decreasing hazard rates, depending on its shape parameter...
295
One-Compartment Open Model: Wagner-Nelson and Loo Riegelman Method for ka Estimation
257
This lesson introduces two critical methods in pharmacokinetics, the Wagner-Nelson and Loo-Riegelman methods, used for estimating the absorption rate constant (ka) for drugs administered via non-intravenous routes. The Wagner-Nelson method relates ka to the plasma concentration derived from the slope of a semilog percent unabsorbed time plot. However, it is limited to drugs with one-compartment kinetics and can be impacted by factors like gastrointestinal motility or enzymatic degradation.
On...
On...
257
Binomial Probability Distribution
10.1K
A binomial distribution is a probability distribution for a procedure with a fixed number of trials, where each trial can have only two outcomes.
The outcomes of a binomial experiment fit a binomial probability distribution. A statistical experiment can be classified as a binomial experiment if the following conditions are met:
There are a fixed number of trials. Think of trials as repetitions of an experiment. The letter n denotes the number of trials.
There are only two possible outcomes,...
The outcomes of a binomial experiment fit a binomial probability distribution. A statistical experiment can be classified as a binomial experiment if the following conditions are met:
There are a fixed number of trials. Think of trials as repetitions of an experiment. The letter n denotes the number of trials.
There are only two possible outcomes,...
10.1K
Hardy-Weinberg Principle
71.3K
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.3K


