Multivariate Genomewide Association Analysis with IHT
Chu, B. B.; Ko, S.; Zhou, J. J.; Zhou, H.; Sinsheimer, J. S.; Lange, K. L.
Show abstract
1In genome-wide association studies (GWAS), analyzing multiple correlated traits is potentially superior to conducting multiple univariate analyses. Standard methods for multivariate GWAS operate marker-by-marker and are computationally intensive. We present a penalized regression algorithm for multivariate GWAS based on iterative hard thresholding (IHT) and implement it in a convenient Julia package MendelIHT.jl (https://github.com/OpenMendel/MendelIHT.jl). In simulation studies with up to 100 traits, IHT exhibits similar true positive rates, smaller false positive rates, and faster execution times than GEMMAs linear mixed models and mv-PLINKs canonical correlation analysis. On UK Biobank data, our IHT software completed a 3-trait joint analysis in 20 hours and an 18-trait joint analysis in 53 hours, requiring up to 80GB of computer memory. In short, our software enables geneticists to fit a single regression model that simultaneously considers the effect of all SNPs and dozens of traits.
Matching journals
The top 5 journals account for 50% of the predicted probability mass.
Similar papers in this journal
Similar papers in this journal
Similar papers in this journal
"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.