The GCTA program as described by Yang and colleagues (Yang et al. 2011) was used to estimate the proportion of phenotypic variance for maxdrinks explained by the genome-wide SNPs genotyped on the Illumina 1M chip in the SAGE dataset. This program fits a linear model to the phenotype data and uses a restricted maximum likelihood to estimate the variance explained by the SNPs.
