When performing association tests on dosage data, it looks like my poorly imputed variants (Rsq < .3ish) get NA pvalues. I would usually filter those out anyways, but out of curiosity is there an internal cutoff in PLINK that excludes sites with low Rsq from the test?
Info Score > 1 Plink?
What does an info score > 1 mean? I thought all scores should be between 0 and 1. Info score comes from assoc.dosage file in plink.
• 5,244 views
•
link
1 answer
The INFO field in an .assoc.dosage file contains the MaCH r-squared imputation quality metric, which is [observed dosage variance] / [expected variance, given observed allele frequencies and Hardy-Weinberg equilibrium]. Despite the r-squared name, this is not guaranteed to be less than or equal to 1, especially when there's significant deviation from Hardy-Weinberg equilibrium; but low values can still be trusted to correspond to poorly-imputed SNPs.
• 0 views
•
link
• 0 views
•
link
Log in to answer this question.
did you check the plink manual pages:
...I put in bold the important line.
Kevin
INFO has the same meaning in both files, but that documentation snippet is actually for a slightly different (and defunct) command, --proxy-assoc instead of --dosage.