Question: filter IMPUTE2 genotype probabilities on .metrics INFO value?
3
gravatar for Krisr
5.7 years ago by
Krisr460
United States
Krisr460 wrote:

HI,

 I have a dataset from dbGAP that provides imputed genotype probability files (.gprobs) for each human chromosome generated with IMPUTE2.  Furthermore, a ".metrics" file is also included.  The documentation says the INFO column within the .metrics file can be used to QC/filter low quality SNPs (INFO < .3) from the .gprob files.   

I am having trouble figuring out which program can handle these 2 files for the purposes of filtering out those SNPs which have a INFO score < .3.    I've looked at QCTOOL, fcgene, and GTOOLs.  However, .gprob seems to be a beagle format, and I can not use it with the .metrics file (which is IMPUTE2) without each program throwing an error(s).  I can probably write a perl script to do this, but given  the apparent frequency of this task, I figured there must be a method for it already.  

So. my question, is there a utility that can perform this task?   I found a paper which did so, but they do not describe how,      "We performed downstream analysis of the complete, imputed merged dataset to take into account the uncertainty of the imputed genotypes. We filtered data based on info score of 0.7 after looking at the distribution of markers at all possible info scores. "         http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4263197/

Any help would be greatly appreciated!    

gprob dbgap imputation impute2 gwas • 4.5k views
ADD COMMENTlink modified 5.7 years ago • written 5.7 years ago by Krisr460
1
gravatar for Krisr
5.7 years ago by
Krisr460
United States
Krisr460 wrote:

ah, figured this out.    

The data to calculate the INFO metric is contained within the gprobs file.  Executing the below command determines the INFO metric for each SNP, then applies the filter.      

./qctool -g chr22.gprobs  -og subsetted.gen   -info .3 1 

ADD COMMENTlink modified 5.7 years ago • written 5.7 years ago by Krisr460
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1874 users visited in the last hour