graphld.blup

June 5, 2026 · View on GitHub

Best Linear Unbiased Prediction (BLUP) for effect size estimation.

Under the infinitesimal model with per-s.d. effect sizes βN(0,D)\beta \sim N(0, D), the BLUP effect sizes are:

E(β)=nD(nD+R1)1R1zE(\beta) = \sqrt{n} D (nD + R^{-1})^{-1} R^{-1}z

where R1R^{-1} is approximated with the LDGM precision matrix. Public BLUP workflows take heritability for the analyzed variant scope and set trace(D)\mathrm{trace}(D) equal to that value across matched LDGM effect indices.

For usage examples, see the BLUP guide.

::: graphld.blup options: show_root_heading: true members_order: source