Monte Carlo estimation of mixed models for large complex pedigrees

Area of Science:

• Human quantitative genetics
• Statistical genetics
• Computational biology

Background:

• Computational complexity limits mixed model estimation in human quantitative genetics to small pedigrees.
• Large pedigrees offer more genetic transmission information and homogeneity than pooled nuclear families.
• Existing methods struggle with the practical reality of large, complex pedigrees.

Purpose of the Study:

• To develop a computationally feasible method for estimating mixed models with major gene effects in large pedigrees.
• To provide a robust approach for analyzing genetic transmission in complex family structures.
• To enable more accurate genetic parameter estimation using extensive pedigree data.

Main Methods:

• A novel Monte Carlo method combining the EM algorithm and Gibbs sampler.
• Joint application of Expectation-Maximization (EM) algorithm and Gibbs sampling for parameter estimation.
• Monte Carlo estimation of the asymptotic variance-covariance matrix for model parameters.

Main Results:

• The proposed method successfully estimates mixed models in large pedigrees.
• The approach provides a Monte Carlo estimate of the variance-covariance matrix.
• The methods are conceptually simple, easy to implement, and versatile.

Conclusions:

• The Monte Carlo EM-Gibbs sampler method overcomes computational barriers in human quantitative genetics.
• This approach enhances the analysis of genetic models in large, complex pedigrees.
• The method supports multiple heritable and nonheritable random components, offering broad applicability.