Rxivist logo

Statistical properties of simple random-effects models for genetic heritability

By David Steinsaltz, Andrew Dahl, Kenneth W Wachter

Posted 14 Nov 2016
bioRxiv DOI: 10.1101/087304 (published DOI: 10.1214/17-EJS1386)

Random-effects models are a popular tool for analysing total narrow-sense heritability for simple quantitative phenotypes on the basis of large-scale SNP data. Recently, there have been disputes over the validity of conclusions that may be drawn from such analysis. We derive some of the fundamental statistical properties of heritability estimates arising from these models, showing that the bias will generally be small. We show that that the score function may be manipulated into a form that facilitates intelligible interpretations of the results. We use this score function to explore the behavior of the model when certain key assumptions of the model are not satisfied -- shared environment, measurement error, and genetic effects that are confined to a small subset of sites -- as well as to elucidate the meaning of negative heritability estimates that may arise. The variance and bias depend crucially on the variance of certain functionals of the singular values of the genotype matrix. A useful baseline is the singular value distribution associated with genotypes that are completely independent --- that is, with no linkage and no relatedness --- for a given number of individuals and sites. We calculate the corresponding variance and bias for this setting.

Download data

  • Downloaded 809 times
  • Download rankings, all-time:
    • Site-wide: 17,101 out of 94,912
    • In genetics: 1,084 out of 4,824
  • Year to date:
    • Site-wide: 54,250 out of 94,912
  • Since beginning of last month:
    • Site-wide: 45,857 out of 94,912

Altmetric data

Downloads over time

Distribution of downloads per paper, site-wide


Sign up for the Rxivist weekly newsletter! (Click here for more details.)