The evaluation of the significance of the association of the individual SNPs with theta EROs is determined by the selection of individual p-values based on control of the false discovery rate. Given the fact that this is a candidate gene study, in which at least a small proportion of the tested variables are expected to be significant, the control of the false discovery rate by the characteristics of data is appropriate. The method used to determine the false discovery rate for any particular level of p-value must be capable of dealing with correlated data, as p-values are correlated as a result of both the non-parametric regression method and the phenotypic correlations described in Chorlian et al. (2015) and summarized in section 3.1.2. The method used is derived from that in Storey and Tibshirani (2003), which uses the distribution of the calculated p-values to determine the prevalence of results which reflect the absence of association, and correspondingly, the prevalence of results which reflect the presence of association. (The part of the original algorithm which estimates the proportion of the non-associated results
