Genomic Correction
2
2
Entering edit mode
12.0 years ago
makaonte ▴ 70

What does it mean "Genomic Correction" when you have a data set with SNPs and p-values from a GWAS (Genome Wide Association Study)? What you do with the data ?

Thanks in advance for any help help.

gwas genomics genetics • 8.0k views
ADD COMMENT
0
Entering edit mode

More details please, and provide a reference to the term to put it into context.

ADD REPLY
0
Entering edit mode

It is odd but the term 'Genomic Correction' came from a colaborator when I send to him a file with a long list of SNPs with their p-values from a GWAS.

ADD REPLY
0
Entering edit mode

Maybe you are referring to correction of p-values for the number of genome wide association tests performed.

ADD REPLY
0
Entering edit mode

Without context it could mean almost anything. Perhaps as Michael suggests you mean genomic control?

ADD REPLY
0
Entering edit mode

You mean "genomic control".

ADD REPLY
4
Entering edit mode
12.0 years ago
Genotepes ▴ 950

Hi

basically you are correction for an inflation of your test statistic (here in Genome wide association studies). If you assume that most of your SNPs are NOT assciated with the trait of interest, then the overwhelming majority of statistical tests'p-values (one per SNP) should follow a uniform distribution (H0). Similarly, the vast majority of your tests (this is mainly in this context that we use GC) should follow of chi-square with a mean of 1. Even in the presence of, say 100 true positives, the mean of all your tests should not be very far from 1. If they are, this can be due to population stratification, diffrential DNA quality (in case controls), or even to polygenic component (infinity of very low effects).

Whatever the reason, this can give an excess of false positive results and on way to correct is to create a new statistic based on the mean of your chi-square distribution (or the median or other robust moments). The simplest way is to divide all your test statistics by the mean of the test statistics.

for SNP i, if your statistic is 2 and and the mean of all your stats is 1.3 (typical of an inflation in case-control studies), then your new corrected statistics will be

2/1.3 = 1.58. This is one of all the possible corrections.

Your new statistic is weaker than the original one (usually if the mean is below 1, we are not correcting - ie not creating a higher chi-square).

Hope this helps

Christian

ADD COMMENT
0
Entering edit mode

Thanks, where could I read more on this?

ADD REPLY
1
Entering edit mode

Usually papers by Bacanu and Devlin are good (advanced) introductions to this : http://www.ncbi.nlm.nih.gov.gate2.inist.fr/pubmed/11754464 http://www.ncbi.nlm.nih.gov.gate2.inist.fr/pubmed/15514657

ADD REPLY
0
Entering edit mode

Thanks for answer!

ADD REPLY
1
Entering edit mode
12.0 years ago
makaonte ▴ 70

From another source I got the following answer: "genomic correction. it's a trick to account for inflation (meaning too low p-values) due to disproportionate relatedness in the population." Does it sound correct ?

ADD COMMENT
1
Entering edit mode
ADD REPLY

Login before adding your answer.

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