Procedure 1
SNP data set: 1000 Genomes Pilot 1 r2 threshold: 0.7 Population panel: CEU Distance limit: 100
plink --bfile bed_DB --epistasis --epi1 1 --epi2 1 --set-test --set sets.txt --allow-no-sex --out bed_DB_sets
$$ p = \phi \left( \frac{1}{\sqrt{k}} \sum_{i=1}^{k} \phi^{-1}(p_i) \right) $$ where $ \phi $ is the Normal distribution cumulative function.