Allelic Association Test

The Allelic Association Operation will be performed along with a report and 2 charts. Only Markers and Samples that have passed the basic Matrix QA and the Hardy-Weinberg test, filtered by specified missingness thresholds, HW equilibrium conformity as well as non-mismatching markers status will be considered in this test.

You may find further details on how to perform basic statistical analysis in a population-based genetic association case-control study in this Nature Protocls article.

Using the recount of Genotypes done at the Genotype Frequency count, a 3×2 contingency table is considered for each marker as in the following example:

Observed AA Aa aa Row Total
Case 110 85 15 210
Control 90 72 48 210
Column Total 200 157 63 420

Expected genotype values must be calculated for every cell, as per the following formulas:

expected case AA = (observed case Row Total × observed AA Column Total) / total genotypes
expected case Aa = (observed case Row Total × observed Aa Column Total) / total genotypes
expected case aa = (observed case Row Total × observed aa Column Total) / total genotypes
expected control AA = (observed control Row Total × observed AA Column Total) / total genotypes
expected control Aa = (observed control Row Total × observed Aa Column Total) / total genotypes
expected control aa = (observed control Row Total × observed aa Column Total) / total genotypes

Resulting in following expected values table:

Expected AA Aa aa
Case 110 85 15
Control 45 36 24

We determine each allele’s frequency by Case/Control:

obsCaseA = 2 × observed case AA count + observed case Aa count
obsCasea = 2 × observed case aa count + observed case Aa count
expCaseA = 2 × expected case AA count + expected case Aa count
expCasea = 2 × expected case aa count + expected case Aa count


obsCtrlA = 2 × observed control AA count + observed control Aa count
obsCtrla = 2 × observed control aa count + observed control Aa count
expCtrlA = 2 × expected control AA count + expected control Aa count
expCtrla = 2 × expected control aa count + expected control Aa count

From this GWASpi calculates the X² value as follows:


X² = ((obsCaseA-expCaseA)/expCaseA)² +
    ((obsCasea-expCasea)/expCasea)² +
    ((obsCtrlA-expCtrlA)/expCtrlA)² +
    ((obsCtrla-expCtrla)/expCtrla)²

From this X² value, a p-value can be calculated for a X² distribution with 1 degree of freedom.

Allelic Association Odds Ratio
The Odds Ratio (OR) is calculated as follows:


1) Check if any value is 0
if YES: OR = NaN


if NOT:

A a
Case 305 115
Control 252 168


OR = Case A × Control a / Case a × Control A,


Then:
If OR is < 1:
OR = 1/OR

See Association Test Reports section for details on table displays and charts.