Heritability Calculator 
=======================

This package acompanies the paper: 

"The Mystery of Missing Heritability: Genetic Interactions Create Phantom Heritability" 
O. Zuk, E. Hechter, S. Sunyaev and E.S. Lander 
Proceedings of the National Academy of Sciences (2012) 



The package enables heritability calculations using different mehotds and for many differnet disease models (genetic architectures)


enables simulations of different sequencing pooling experiments, including pooling design, modeling sequencing errors, reconstruction of genotype vectors and evaluating reconstruction accuracy. 
It may serve as a guideline for designing pooling sequencing experiments.
One may also enter real allele counts collected in pooled experiments and reconstruct the genotype.

Please look at the following files for the different applications:

(1) test_estimate_heritability_using_IBD.m - A script testing heritability estimation by local-regression on IBD. The script simulates IBD structure, genotypes 
and phenotypes, with phenotypes being a non-linear noisy function of genotypes. Then, we estimate narrow-sense heritability by regressing phenotypic correlation 
on IBD sharing and compare to the true narrow-sense heritability.

(2) compute_k_of_N_liabilities_statistics.m - computing all relevant genetic epidemiological parameters for a generalized Limiting-Pathway genetic architecture.

(3) estimate_IBD_sharing_from_snps.m - Estimate IBD sharing between pairs of individuals (used for heritability calculations).

(4) compute_association_power.m - Power calculations for different tests (e.g. marginal effect/epistasis/pathway-epistasis tests) and different genetic architectures. 


The package was tested on Matlab 2011a, but should work on most Matlab versions. 


Please send any comments/bug reports to: Or Zuk: orzuk@broadinstitute.org, orzuk@ttic.edu


Ver 2.0: Apr. 29th, 2013