Tunings for leapfrog integration of Hamiltonian Monte Carlo for estimating genetic parameters
Arakawa, A.; Hayashi, T.; Taniguchi, M.; Mikawa, S.; Nishio, M.
Show abstract
A Hamiltonian Monte Carlo algorithm is a Markov Chain Monte Carlo method that is considered more effective than the conventional Gibbs sampling method. Hamiltonian Monte Carlo is based on Hamiltonian dynamics, and it follows Hamiltons equations, which are expressed as two differential equations. In the sampling process of Hamiltonian Monte Carlo, a numerical integration method called leapfrog integration is used to approximately solve Hamiltons equations, and the integration is required to set the number of discrete time steps and the integration stepsize. These two parameters require some amount of tuning and calibration for effective sampling. In this study, we applied the Hamiltonian Monte Carlo method to animal breeding data and identified the optimal tunings of leapfrog integration for normal and inverse chi-square distributions. Then, using real pig data, we revealed the properties of the Hamiltonian Monte Carlo method with the optimal tuning by applying models including variance explained by pedigree information or genomic information. Compared with the Gibbs sampling method, the Hamiltonian Monte Carlo method had superior performance in both models. We have provided the source codes of this method written in the R language.
Matching journals
The top 5 journals account for 50% of the predicted probability mass.
Similar papers in this journal
- Interpretable Artificial Neural Networks incorporating Bayesian Alphabet Models for Genome-wide Prediction and Association Studies 96%
- A Multiple-trait Bayesian Variable Selection Regression Method for Integrating Phenotypic Causal Networks in Genome-Wide Association Studies 96%
- Performing Parentage Analysis for Polysomic Inheritances Based on Allelic Phenotypes 95%
Similar papers in this journal
Similar papers in this journal
- Multifactorial Methods Integrating Haplotype and Epistasis Effects for Genomic Estimation and Prediction of Quantitative Traits 95%
- Validation of the linear regression method to evaluate population accuracy and bias of predictions for non-linear models 95%
- Accelerated matrix-vector multiplications for matrices involving genotype covariates with applications in genomic prediction 93%
Similar papers in this journal
- Using encrypted genotypes and phenotypes for collaborative genomic analyses to maintain data confidentiality 94%
- Estimating SNP heritability in presence of population substructure in large biobank-scale data 94%
- Robust Genomic Modelling of Additive and Non-Additive Variation with Intuitive Inclusion of Expert Knowledge 94%
Similar papers in this journal
- A logratio approach to the analysis of autosomal genotype frequencies across multiple samples 93%
- Estimating contemporary migration numbers of adults based on kinship relationships in iteroparous species 93%
- freqpcr: estimation of population allele frequency using qPCR ΔΔCq measures from bulk samples 93%
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.