Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
Comment
. 2015 Oct 6;112(40):E5449-51.
doi: 10.1073/pnas.1502868112. Epub 2015 Sep 28.

Implications of simplified linkage equilibrium SNP simulation

Affiliations
Comment

Implications of simplified linkage equilibrium SNP simulation

Sang Hong Lee. Proc Natl Acad Sci U S A. .
No abstract available

PubMed Disclaimer

Conflict of interest statement

The author declares no conflict of interest.

Figures

Fig. 1.
Fig. 1.
Artifactual correlation between the eigenvectors of the GRM and disease status generated from the SLES simulation. (A) With the SLES simulation, the association between the eigenvectors and case–control status is unrealistically inflated when the value for N (# individuals)/M (# SNPs) increases. The correlation between the first principal component and disease status is 0.14, 0.70, and 0.63 with the value for N/M = 1.3, 13, and 8, respectively. Population disease risk of K = 0.01 and proportion of cases in the sample of P = 0.5 were used. Red represents cases, and blue represents controls. (B) With the GLDS simulation, the association between the eigenvectors and case–control status is negligible, regardless of the values for N/M, i.e., more realistic compared with the SLES. The correlation between the first principal component and disease status is 0.04, 0.03, and 0.06, with the value for N/M = 1.3, 13, and 130, respectively (GLDS could not simulate n = 8,000, because a GWAS data set of ∼400,000 individuals would be needed; instead, we tested it with an extreme with M = 10, i.e., N/M = 130). Population disease risk of K = 0.01 and proportion of cases in the sample of P = 0.5 were used. Red represents cases, and blue represents controls.
Fig. 2.
Fig. 2.
REML bias is negligible with a more realistic simulation and real data. (A) The average of estimated SNP heritability and empirical SE bar of the mean estimate from REML with GLDS simulation (blue line) and SLES simulation (red line) over 50 replicates. The true simulated SNP heritability is 0.5. (B) Estimated SNP heritability from REML and PCGC with real data analyses [to be compared with figure 2B in Golan et al. (1)]. We excluded two diseases that had highly confounded population structure [figure S1 in Gusev et al. (3)]. HE, Haseman–Elston regression. (C) Estimation error assuming that PCGC estimate is true value [to be compared with figure S4 in Golan et al. (1)].

Comment in

Comment on

References

    1. Golan D, Lander ES, Rosset S. Measuring missing heritability: Inferring the contribution of common variants. Proc Natl Acad Sci USA. 2014;111(49):E5272–E5281. - PMC - PubMed
    1. Lee SH, Wray NR, Goddard ME, Visscher PM. Estimating missing heritability for disease from genome-wide association studies. Am J Hum Genet. 2011;88(3):294–305. - PMC - PubMed
    1. Gusev A, et al. Schizophrenia Working Group of the Psychiatric Genomics Consortium SWE-SCZ Consortium Schizophrenia Working Group of the Psychiatric Genomics Consortium SWE-SCZ Consortium Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases. Am J Hum Genet. 2014;95(5):535–552. - PMC - PubMed
    1. Thompson EA. Identity by descent: Variation in meiosis, across genomes, and in populations. Genetics. 2013;194(2):301–326. - PMC - PubMed
    1. Chen G-B. Estimating heritability of complex traits from genome-wide association studies using IBS-based Haseman–Elston regression. Front Genet. 2014;5:107. - PMC - PubMed

LinkOut - more resources