对非常大的数据集和复杂的单步基因组模型的 (协同) 差异组件的估计
Matias Bermann1, Andres Legarra2,3, Ignacio Aguilar4
1Department of Animal and Dairy Science, University of Georgia, Athens, GA, 30602, USA. mbermann@uga.edu.
Genetics, selection, evolution : GSE
|October 31, 2025
概括
大型基因组数据集的估计方差组件的计算局限性被蒙特卡洛单步基因组REML (MC-ss-GREML) 克服. 这种新方法使用显著减少计算时间和内存准确估计差异组件,使复杂的遗传分析可行.
科学领域:
- 动物育种与遗传学
- 统计遗传学 统计遗传学
- 计算生物学 计算生物学
背景情况:
- 准确估计线性混合模型中的方差成分对于遗传分析至关重要.
- 计算约束通常需要数据子集或模型简化,可能引入偏差.
- 现有的蒙特卡罗REML (MC-REML) 方法缺乏用于单步基因组分析的扩展.
研究的目的:
- 扩展蒙特卡罗REML (MC-REML),将大型基因组数据集纳入单步基因组最佳线性无偏预测 (ssGBLUP) 模型中.
- 开发一种计算效率高的方法,用于大规模遗传评估中估计差异组件.
主要方法:
- 开发了蒙特卡洛单步基因组REML (MC-ss-GREML),通过模拟繁殖值和解决混合模型方程.
- 使用预期最大化和平均信息来优化REML.
- 使用三特征肉牛生长模型和大出生体重模型验证了该方法.
主要成果:
- 与精确ss-GREML相比,MC-ss-GREML在差异成分估计中没有差异.
- 与精确方法相比,该方法减少了86%的计算时间和99%的内存使用量.
- 证明了可扩展性,使用700万只动物的大数据集,在11轮汇聚.
结论:
- MC-ss-GREML有效地估计了大型的基因型种群的变异成分.
- 该方法可显著降低计算成本 (时间和内存).
- 通过复杂的模型和广泛的基因组数据,能够进行准确的遗传评估.
相关概念视频
One-Compartment Open Model: Wagner-Nelson and Loo Riegelman Method for ka Estimation
1.1K
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...
1.1K
Genome-wide Association Studies-GWAS
15.3K
Genome-wide association studies or GWAS are used to identify whether common SNPs are associated with certain diseases. Suppose specific SNPs are more frequently observed in individuals with a particular disease than those without the disease. In that case, those SNPs are said to be associated with the disease. Chi-square analysis is performed to check the probability of the allele likely to be associated with the disease.
GWAS does not require the identification of the target gene involved in...
GWAS does not require the identification of the target gene involved in...
15.3K
Estimating Population Standard Deviation
3.3K
When the population standard deviation is unknown and the sample size is large, the sample standard deviation s is commonly used as a point estimate of σ. However, it can sometimes under or overestimate the population standard deviation. To overcome this drawback, confidence intervals are determined to estimate population parameters and eliminate any calculation bias accurately. However, this only applies to random samples from normally distributed populations. Knowing the sample mean and...
3.3K
Estimating Population Mean with Unknown Standard Deviation
8.8K
In practice, we rarely know the population standard deviation. In the past, when the sample size was large, this did not present a problem to statisticians. They used the sample standard deviation s as an estimate for σ and proceeded as before to calculate a confidence interval with close enough results. However, statisticians ran into problems when the sample size was small. A small sample size caused inaccuracies in the confidence interval.
William S. Gosset (1876–1937) of the...
William S. Gosset (1876–1937) of the...
8.8K
Evolutionary Relationships through Genome Comparisons
6.8K
Genome comparison is one of the excellent ways to interpret the evolutionary relationships between organisms. The basic principle of genome comparison is that if two species share a common feature, it is likely encoded by the DNA sequence conserved between both species. The advent of genome sequencing technologies in the late 20th century enabled scientists to understand the concept of conservation of domains between species and helped them to deduce evolutionary relationships across diverse...
6.8K
Estimating Population Mean with Known Standard Deviation
9.6K
To construct a confidence interval for a single unknown population mean μ, where the population standard deviation is known, we need sample mean as an estimate for μ and we need the margin of error. Here, the margin of error (EBM) is called the error bound for a population mean (abbreviated EBM). The sample mean is the point estimate of the unknown population mean μ.
The confidence interval estimate will have the form as follows:
(point estimate - error bound, point estimate +...
The confidence interval estimate will have the form as follows:
(point estimate - error bound, point estimate +...
9.6K


