myvars = c("FID1", "IID1", "FID2", "IID2", "PI_HAT")
problem_pairs[myvars]
plink --file GWAS_clean3 --pheno pheno.txt --pheno-name Aff --hardy
q()
plink --file GWAS_clean2 --remove IBS_excluded.txt --recode --out GWAS_clean3
plink --file GWAS_clean3 --pheno pheno.txt --pheno-name Aff --hardy