I am trying to do a GWAS study, where my pedigree info and marker data is coming from axiom array. I have .ped and .map plink compatible files. I converted them to binary .bed, .fam and .bim files for plink usage.
I have to perform QC on these files in plink before imputation. I want to separately perform QC on sample first and then marker QC. And most importantly, at each step of QC, I want to have a separate file. Not overwritten. I want to automate this pipeline as I have huge data. Is there any way I can do this and get separate files independent of one another for every step? In plink, there is an option "--script". For that as per I understood, we need to provide all filters together with specified thresholds and at the end it will give us one filtered file. Is there any other option other than this? Or one has to write python script for automating? Pardon me, I am not good at python.
Any help is highly appreciated.