Info Score > 1 Plink?
1
0
Entering edit mode
2.9 years ago
bg1116 • 0

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.

Plink GWAS Score • 1.8k views
ADD COMMENT
0
Entering edit mode

did you check the plink manual pages:

The fields INFO and NPRX refer to how well PLINK managed, if at all, to impute the SNP. If NPRX is zero, then it could not be even poorly imputed. If INFO ranges from between 0 and 1, although it can be greater than 1 occasionally. A higher value general means a better imputed SNP; roughly speaking, only looking at imputed SNPs with a INFO value greater than 0.8 or so is probably good practice. More specific details on these metrics will be posted soon.

...I put in bold the important line.

Kevin

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
1
Entering edit mode
2.9 years ago

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.

ADD COMMENT
0
Entering edit mode

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?

ADD REPLY

Login before adding your answer.

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