I would like to know the association between the snps and the clinical variables available. For this, I have to do association for all possible combinations of one snp- one variable. I have come across an option in Plink for this. But to follow this process for 100+ snps and 150+ variables is cumbersome. Would you suggest any other option in plink/R/any other software ?