How to handle "no call" on SNP imputation results
0
0
Entering edit mode
6 weeks ago
Babbala • 0

Hi, i just imputed some SNPs using IMPUTE2, i pruned them and then kept all the haplotypes with probabilities higher than 90% (thresholding was done using Plink1.90)

but what i get is that some SNPs that do not pass the threshold are undefined with the symbol ./. (instead of the normal 0/0, 0/1 or 1/1)

I don't know if what i am doing is right and i need some advice on how to handle them because i would like to keep as much SNPs as possible and not removing samples.

I am new about imputation so any feedback is appreciated.

Thank you.

GWAS SNP genotype Imputation • 166 views
ADD COMMENT
0
Entering edit mode

Can you please show the code that you have used for pre-phasing, imputation, and then post-conversion of the imputed data to VCF? ./. means that, while a variant at this position was called in one or more other samples in the dataset, in the sample with ./. there was not enough information to make any call, usually no aligned reads.

ADD REPLY

Login before adding your answer.

Traffic: 2194 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