how to do QC on data imputed on Michigan Imputation Server?
0
3
Entering edit mode
4.6 years ago
anamaria ▴ 220

Hello,

I had .CEL files and than I converted those to plink format and did QC steps following https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6001694/ Then imputed those SNPs on Michigan imputation server and than I did QC steps on those file like this:

plink --vcf  chr1.dose.vcf.gz --biallelic-only --make-bed --double-id --out s1
plink --bfile s1 --bmerge s1 --merge-mode 6
plink --bfile s1 --exclude plink.missnp --make-bed --out s2
plink --bfile s2 --list-duplicate-vars
plink --bfile s2 --exclude plink.dupvar --make-bed --out s3
plink --bfile s3 --qual-scores chr1.info 7 1 1 --qual-threshold 0.5 --make-bed --out s4
plink --bfile s4 --maf 0.01 --hwe 0.000001 --snps-only --make-bed --out s5
plink --bfile s5  --geno 0.02 --mind 0.02 --make-bed --out FINAL_DATA_QC1

Is there is anything you would change in these parameters or do you have any GWAS QC workflow available which would be more relevant?

Thanks!

plink gwas imputation • 1.9k views
ADD COMMENT

Login before adding your answer.

Traffic: 2000 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6