This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Selecting SNPs to impute in Plink

Hi all,

I want to construct a workflow that will allow me to impute small batches of SNPs that were untyped in my GWAS data (in .bed .bim .fam fomat).

Is it possible to meaningfully subset my reference panel (hapmap release 23) so that I only impute the ~8000 missing SNPs? Would it be erroneous to start by subsetting the reference by tag SNPs to those 8000 SNPs, then merging the reference and GWAS data, then imputing genotypes?

Thanks for the help in advance!

plink impute
• 3,561 views • link
updated by Ram Ram • written by hessjl • hessjl

1 answer

Is there a reason you still want to use Plink here? Its imputation algorithm was already considered obsolete 5 years ago: see http://www.ncbi.nlm.nih.gov/pubmed/19089453.

I suggest choosing a more effective imputation tool and applying it to your data. BEAGLE, IMPUTE2, and MaCH continue to be reasonable choices, and it should be straightforward to convert your data to their formats. (1000 Genomes phase 3 reference panels are also now available for all three programs; they should yield more accurate results than HapMap release 23.)

• 0 views • link
updated by Ram Ram • written by chrchang523 chrchang523

Thanks for this insight! As subtext, are you indicating that whole GWAS imputation is the way to go; or can I still pre-filter my reference set to expedite the imputation turnaround?

• 0 views • link
updated by Ram Ram • written by hessjl • hessjl

Log in to answer this question.

More posts like this