and 1, respectively). We then calculated the frequency of concordance between the permuted data and the observed data. Finally, we determined the MLE by plotting the smoothed trend line of frequency versus number of risk genes using local polynomial regression fitting (loess in R) and used the predict function in R to estimate the MLE of the number of risk genes (420 genes).