 Technical advance
 Open Access
 Open Peer Review
 Published:
The importance of distinguishing between the odds ratio and the incidence rate ratio in GWAS
BMC Medical Genetics volume 16, Article number: 71 (2015)
Abstract
Background
In recent years, genome wide association studies have identified many genetic variants that are consistently associated with common complex diseases, but the amount of heritability explained by these risk alleles is still low. Part of the missing heritability may be due to genetic heterogeneity and small sample sizes, but nonoptimal study designs in many genome wide association studies may also have contributed to the failure of identifying gene variants causing a predisposition to disease. The normally used odds ratio from a classical casecontrol study measures the association between genotype and being diseased. In comparison, under incidence density sampling, the incidence rate ratio measures the association between genotype and becoming diseased. We estimate the differences between the odds ratio and the incidence rate ratio under the presence of events precluding the disease of interest. Such events may arise due to pleiotropy and are known as competing events. In addition, we investigate how these differences impact the association test.
Methods
We simulate life spans of individuals whose gene variants are subject to competing events. To estimate the association between genotype and disease, we applied classical casecontrol studies and incidence density sampling.
Results
We find significant numerical differences between the odds ratio and the incidence rate ratio when the fact that gene variant may be associated with competing events, e.g. lifetime, is ignored. The only scenario showing little or no difference is an association with a rare disease and no other present associations. Furthermore, we find that pvalues for association tests differed between the two study designs.
Conclusions
If the interest is on the aetiology of the disease, a design based on incidence density sampling provides the preferred interpretation of the estimate. Under a classical casecontrol design and in the presence of competing events, the change in pvalues in the association test may lead to false positive findings and, more importantly, false negative findings. The ranking of the SNPs according to pvalues may differ between the two study designs.
Background
The genome wide association study (GWAS) is a powerful tool to associate genetic variation with disease, and thousands of associations have been established and replicated [1]. Through the development of new methodology, the results from the GWAS's have been useful in other settings. One example is in healthrelated epidemiology studies (e.g. Ligthart et. al. [2] and Simonson et. al [3]). Here, the estimated odds ratio from the GWAS was used as a proxy for genetic variation through the polygenic score [4, 5]. The estimated value of the odds ratio is most often based on a number of unrelated cases and controls, and estimates the association between genetic variation and the prevalence of disease.
However, events precluding the disease of interest may play an undesired and crucial role when sampling cases and controls for a GWAS. If some of the genotyped SNPs are correlated with a disease but also with other events precluding this e.g. early death, a carefully chosen study design is needed to identify true associations with the aetiology of the disease. Precluding events may arise due to pleiotropy and are often referred to as competing events [6, 7]. These possible (and most often hidden) associations can be taken into account using incidence density sampling matching on age. This design measures the association between SNPs and acquisition of the disease and the effect size is expressed as the incidence rate ratio. Whether or not to take competing events into account in a GWAS depends on the scientific question. The classical casecontrol study provides a measure of the association between genotype and being diseased, whereas the incidence density sampled study associates genotype with becoming diseased.
In this paper, we explain two widely used measures of association: the odds ratio (OR) and the incidence rate ratio (IRR) [8]. We explain, using an analogy with leaves falling from a tree, the different interpretations of the prevalence and the incidence rate, and how competing events can affect the interpretation of disease associations. In order to understand and quantify in detail the difference between the OR and the IRR of disease we consider two sampling methods for controls: incidence density sampling (or matching) and unrelated casecontrol sampling. We use simulations to quantify the difference between the association estimates from the two study designs under different effects of competing events.
We find considerable numerical differences between the OR and the IRR. The differences are reflected in a change in the pvalues of the association tests. Our findings have significant implications for GWASs, since they point out the risk of serious false negative associations caused by competing events.
Methods
Prevalence vs. Incidence rate
To illustrate the prevalence, the incidence rate, and the effect of competing events, one may imagine a large tree in a park. The tree has two kinds of leaves: heavy leaves and light leaves. Consider a dominant model, corresponding to a population of individuals with two levels of one exposure (presence or absence of risk allele). When leaves fall from the tree they can either fall onto the sticky soil or onto the nonsticky concrete. Leaves falling on the concrete will eventually move to the soil. Falling onto the soil will decompose the leaf, corresponding to dying, and falling onto the concrete corresponds to becoming sick. The healthy leaves are still hanging on the tree. The heavy leaves will not be carried by the wind as easily as the light leaves, and therefore the heavy leaves have a tendency to fall on the soil under the tree more often than the light leaves, and the light leaves have a tendency to fall onto the concrete. Among leaves falling on the concrete, the heavy leaves have a tendency to stay longer on the concrete, before they are blown onto the soil, than the light leaves.
The number of leaves falling from the tree onto the concrete per time unit corresponds to the incidence rate of disease i.e., the light leaves have a higher incidence rate than the heavy leaves. The number of leaves falling onto the soil per time unit corresponds to the incidence rate of death i.e. the light leaves have a smaller incidence rate of death than the heavy leaves. The prevalence of disease is the fraction of leaves on the concrete at a given time period compared to the total number of leaves on the tree and on the concrete. The odds of disease is the number of leaves on the concrete at a given time point compared to the number of leaves on the tree at the same time, i.e. the prevalence of disease divided by one minus the prevalence of disease.
The OR, normally calculated in a GWAS, is a function of the prevalence. It is the ratio between the odds of disease for light leaves compared to the odds of disease for heavy leaves. It measures the association between the exposure and being diseased. The IRR is the fraction of the incidence rate for light leaves vs. the incidence rate for heavy leaves, and it measures the association between the exposure and becoming diseased. The OR depends on the incidence rate of leaves falling onto the concrete and the incidence rate of leaves moving onto the soil both from the tree and from the concrete. Because of these dependencies SNPs may be associated with a disease through the OR, but not through the IRR and vice versa.
Measures of association
When the exposure or the outcome is not available for all individuals in a population, it is impossible to set up a full cohort study, and different types of sampling are used. Most often the measure of disease occurrence is calculated as the number of cases divided by the number of controls, maybe for different exposure groups. The measure of association between exposure and disease is then calculated as the ratio between the measures of disease occurrence in two different groups of the exposure [9].
When using a casecontrol sample to estimate the association between an exposure and a disease of interest, the interpretation of the association depends on the procedure used to select controls. Two different methods are normally used to select controls according to the onset of disease: the unrelated casecontrol sampling and the incidence density sampling.
The two sampling methods are displayed in Fig. 1. Each horizontal line in the figure corresponds to a lifespan of an individual. The individual is born when the line starts and dies when the line ends. A black diamond marks time of onset of the disease of interest. The green hatched area corresponds to a sampling period for capturing cases and controls for a classical casecontrol study. The green and red dots mark individuals, who are eligible controls for the two different sampling methods described below.

1)
In the classical casecontrol study all cases are chosen within a certain time period. The sampling time is specified by the investigators. The controls are chosen uniformly at random among those who have not experienced the event after the sampling period (green dots in Fig. 1). The number of cases proportional to the number of controls estimates the odds of disease, and the OR is estimated as the fraction between two odds. The measure of disease occurrence in this setup is the odds, which is a function of the prevalence [10].

2)
Using incidence density sampling controls are sampled longitudinally in time. The controls are chosen uniformly at random among persons at risk of developing the disease at the time when a case appears (red dots in Fig. 1), i.e., each person who develops the disease of interest during the sampling period has their own set of potential controls. This is also called timematching [11]. The measure of disease can be interpreted as the probability of becoming diseased in the next time increment, i.e., the incidence rate of disease. The IRR is the association measure [12]. If the data support an assumption of constant instantaneous probability of disease in small intervals of time, then the controls may instead be chosen within a small time window (and not at specific time points). Contrary to the OR, the IRR is not dependent on the prevalence of the disease.
Simulation study
To quantify the differences between the OR and the IRR, if any, we simulate a longitudinal study with one exposure and two events; disease and death. The longitudinal study includes 10,000 individuals. From this study we generate two different casecontrol studies using the sampling designs described above to collect controls: the classical casecontrol study and the matched casecontrol study using incidence density sampling. In detail we do the following:
For all 10,000 individuals in the simulated data we generate the number of minor alleles under the assumption of HardyWeinberg equilibrium and a minor allele frequency (MAF) of 0.1, 0.25 and 0.5.
In order to generate a longitudinal data source with one exposure, we use the algorithm discussed by Beyersmann et al. [13] using two events: disease and death. An individual’s waiting time to any event is generated as a function of the total incidence rate of the two events over time. At the given time point of an event, the probability of disease is the proportion between the incidence rate of disease and the total incidence rate. Individuals still alive at age 90 are censored.
We use a Cox proportional hazard model where the baseline hazards follow a Gompertz distribution [14] for both events. The baseline for the incidence of dying follows a Gompertz distribution with mode 85 and shape parameter 0.0004. The densities for different effects of the number of minor alleles on the IRR of death are shown in Fig. 2. Two different parameter settings are used for the baseline of disease: one with mode 25 and shape parameter 0.95, and one with mode 50 and shape parameter 0.1. Figure 2 shows the densities for the two parameter settings for different effects of the number of minor alleles on the IRR. The black lines in Fig. 2 correspond to the overall baseline density. We multiply the incidence rate of disease by a constant (less than 1) in order to scale the lifetime risk of disease [15]. Different values of the constant are used to consider common versus rare diseases; the smaller the constant the rarer the disease. We assume a log linear proportional effect of the number of minor alleles on the incidence rate of disease and death, with different incidence rate ratios of 0.5, 1.0, 1.1, 1.2, 1.5, 1.7, 2.0 and 3.0 for disease and 0.5, 0.75, 0.9, 1.0, 1.1, 1.2, 1.5, 1.7, 2.0, and 3.0 for death.
From the simulated longitudinal data we generate two epidemiological study designs to estimate the association between the number of minor alleles and disease. For the incidence density sampling, we choose 5 individuals uniformly at random among the persons not diseased and not dead at the age when an individual becomes diseased. This result in a 1: 5 individually time matched nested casecontrol study, where cases and controls are matched on age. The association between the number of minor alleles and disease is estimated using conditional logistic regression (see chapter 7 in [16]).
The cases in the unrelated casecontrol sampling are individuals diseased before the age of 90. The controls are chosen among all persons not diseased before the age of 90, that is, the sampling period is from birth to age 90 (see Fig. 1). If the number of possible controls is larger than five times the number of cases, the controls are chosen uniformly at random at a 1:5 ratio between cases and controls. We estimate the association between the number of minor alleles and disease in the casecontrol sample using logistic regression (see chapter 6 in [16]). Additionally, we simulate an age covariate for each individual in the classical casecontrol study. For controls, the age is sampled from a uniform distribution on the interval form birth to end of followup. For cases the covariate is the age at which the individual is diagnosed. We then estimate the association between the number of minor alleles and disease adjusted for a linear age trend in the classical casecontrol study.
In total we generate 1000 longitudinal data sources for each configuration of MAF, IRR of disease, IRR of death, the rarity of the disease and the two parameter settings of the baseline rate of disease. The measure of association in a given analysis is estimated as the average of the 1000 replications. The null hypothesis is no association between the number of minor alleles and disease prevalence in the classical casecontrol study, and is no association between minor alleles and disease incidence in the incidence density sampled study. The pvalue of this hypothesis is estimated as 2 times the proportion of estimates less than 0 if the measure of association is larger than 0, and it is estimated as 2 times 1 minus the proportion of estimates less than 0 if the measure of association is less than or equal to 0. The proportion less than 0 are calculated based on the 1000 replicates (see Additional file 1: R code for programming details).
Results
The estimate from the incidence density sampling of controls is not numerically different from the true IRR determined by the simulation study irrespective of the baseline of disease, the IRR of death and the MAF (results not shown).
For the classical casecontrol study, Fig. 3 shows the estimated associations between the disease and the number of minor alleles for each parameter setting of the baseline rate of disease and for a common and a rare disease (the IRR for death equals 0.5, 1.0 and 3.0. The remaining results are given in the online Additional file 2: Figure S2, including corresponding figures for the ageadjusted model in Figure S1 and Figure S3). In each subfigure one estimate is shown for different settings of the true IRR of disease, the IRR of death and the MAF. The probability of disease at age 90, i.e. the cumulative incidence at age 90, is approximately 22 % for the common disease and 2.5 % for the rare disease (see the legend of Fig. 3). The black solid diagonal line indicates no numerical differences between the estimated OR and the true IRR, whereas the black solid horizontal line at 1 indicates no association according to the estimated OR.
Figure 4 shows the pvalues for the null hypothesis of no association for the two different sampling methods and for different values of the true IRR of disease and the IRR of death. The three identical coloursymbol combinations correspond to the different MAF. The scenarios and symbols in Fig. 4 are similar to those in Fig. 3. The pvalues from the matched casecontrol study are plotted on the horizontal axis and the pvalues from the classical casecontrol study are plotted on the vertical axis. For each subfigure we see false positive findings below the diagonal and false negative findings above the diagonal. Scenarios where both pvalues equal zero are not plotted.
If the number of minor alleles is associated with longevity, i.e. the IRR for death is less than 1, then the estimated OR is larger than the true IRR of disease, irrespective of the rarity of the disease, the baseline rate of disease and the MAF. The differences become larger when the value of the true IRR for disease is large, and the differences are larger for common diseases than for rare diseases (see blue symbols in Fig. 3). The larger estimates are reflected in false positive findings. In the case of no true association (the true IRR for disease equals 1) the estimated OR is significantly different from 1 for common diseases (see blue symbols in Fig. 4). Generally, the pvalues from the classical casecontrol study are smaller than the pvalues from the matched casecontrol study.
When there is no association between the number of minor alleles and lifetime, i.e. the IRR of death is 1, the estimated OR is larger than the true IRR for common diseases when the true IRR for disease is larger than 1 (green symbols in Fig. 3a and c). For the rare diseases, there is no difference between the estimated OR and the true IRR of disease (green symbols in Fig. 3b and d). The pvalues from the classical casecontrol study and from the matched casecontrol study are similar, thus all green symbols in Fig. 4 follow the solid black line.
When the true IRR for death is larger than 1, the estimated OR is smaller than the true IRR for disease (red symbols in Fig. 3). In Fig. 3a the estimated OR is larger than the true IRR for death, for large true values. The smaller estimated OR gives rise to false positive associations (see + in Fig. 4a and d). False negative associations occurred when the estimated OR is close to 1 and the true IRR is larger than 1, e.g. when the IRR of the true association with disease is 1.7 (see in Fig. 3 and 4c or d). If the true IRR equals 1.1 (x) Fig. 4c reveals a positive test for association. The same symbols, x, on Fig. 3c shows the estimated OR is significantly smaller than 1, whereas the true association is an increase in the incidence rate of disease (the true IRR equals 1.1).
Including age in the model introduces slightly smaller differences for the estimates in scenario A and B in Fig. 3 (see Additional file 2: Figure S1). For scenario C there are smaller differences when the IRR for death is larger than 1, but a slightly larger difference if the IRR for death is smaller than or equal to 1. The largest effect of the adjustment is seen in scenario D, where the introduction of age as a covariate considerably reduce the differences irrespective of the value of the IRR of death.The differences between the estimated OR and the true IRR for disease, for different values of the MAF and the IRR for death, are displayed in Additional figure 2: Figure S4 and Figure S5 (age adjusted). The relative differences between the two are displayed in Additional figure 2: Figure S6 and Figure S7 (age adjusted).
For scenario D in Figs. 3 and 4, Fig. 5 displays the comparison of the pvalues for all simulated results where the IRR for disease is 1.0, 1.1 and 1.7 (scenarios A, B and C are displayed in the online Additional file 2: Figure S8, S9 and S10). The –log transformed pvalues above the black horizontal line represent the classical casecontrol study, whereas those below the horizontal line originate from the matched casecontrol study. The dashed lines indicate a genomewide significance level of 5 · 10^{−8}. For clarity, all pvalues lower than 10^{−16} are truncated. For large and small values of the IRR for death and for the IRR for death equal to 1.0 and 1.1, the pvalues from the classical casecontrol study are lower than those from the incidence density sampled study, to the extent of genomewide significance. When the IRR for disease equals 1.7, the pvalues from the casecontrol study decrease more rapidly than those from the incidence density sampled study. It falls below the genomewide significance level.
Discussion
For each configuration of the rarity of the disease, the parameter setting of the baseline rate of disease, the MAF, and the IRR of disease and death, we simulated 1000 estimates of the odds ratio of disease. For each configuration, the standardized estimates of the OR for disease are, according to the central limit theorem (e.g. [17]), standard normal distributed.
There is no difference between the estimated and true association of disease and the number of minor alleles using incidence density sampling. In the classical casecontrol study we observe differences between the estimated OR and the true IRR in both directions, depending on the IRR for death. We find a larger estimated OR if there is an association between the number of minor alleles and longevity, i.e., the IRR for death is less than 1. If the association between the number of minor alleles and death is larger than 1, the estimated OR is smaller. However, these differences are small when the effect of the competing association is small. Small effects of competing events (below 1.3) are expected in a GWAS, since no evidence of large effects is found so far. When there is no association between the number of minor alleles and competing events we find a larger estimated OR if the disease is common, whereas there are no differences if the disease is rare. Others find similar results, e.g. Fradin et al. [18] and Wang et al. [19].
The fact that we do not find any differences if the SNP is not associated with death and the disease is rare, is known as the “rare disease assumption”, i.e. the OR is asymptotically equal to the IRR if the disease is rare [9, 20]. Notice that the “rare disease assumption” is compromised if there is an association between the SNP and competing events. An inference on the probability of disease may be biased due to competing events, even in a longitudinal study. Fine and Gray [21] developed a method allowing association testing of the effect of a SNP on the probability of disease.
From Fig. 5 we identify differences in the pvalues both towards false positive hypothesis and false negative hypothesis. These differences are due partly to the larger or smaller estimated OR and to the effect of the total number of cases on the statistical power of the studies. Pvalue differences may introduce different rankings of SNPs between the two study designs. The ranking of the SNPs in a GWAS using a classical casecontrol study gives a list of SNPs most associated with the prevalence of disease, whereas the incidence density sampled study ranks SNPs according to the impact of the SNP on the aetiology of the disease. If the aim of the association study is to identify genetic markers involved in the genetic aetiology of the disease, then incidence density sampling allows for the correct interpretation of the estimates.
When a summary of the OR estimates from a classical casecontrol study is used as a proxy for genetic variation, e.g. in the polygenic score, the differences between the estimated OR and the IRR are aggregated into one measure. It is unknown if and to what extent this aggregation affects the validity of the approximation. The possible change in rankings can lead to summary statistics based on false positive SNPs. The approximation may be affected more, though, by the exclusion of false negative SNPs. The impact of ignoring the false negative SNPs is unknown and request further research.
In studies using a welldocumented dichotomous confounder, the incidence density sampled study has higher power than the classical casecontrol study [22, 23]. However, in practice this may not always be the case. If a predefined incidence density sampled study cannot collect all needed information on both cases and controls the matching may be broken. In a 1:1 matched design this means a loss of both the case and the control in the analysis, reducing the power due to loss of information [24]. In the classical casecontrol study, loss of information for a single subject will not imply a loss of information for other subjects in the study.
The presence of competing events will decrease the number of cases observed, and will decrease the number of available matched sets in the incidence density sampling. The power of a matched study depends on the number of sets, and therefore the power is dependent on the effect of the competing event. The decreasing number of cases will also affect the power of a classical casecontrol study. This effect is, however, assumed to be smaller than under incidence density sampling [25]. The literature in this field of research is sparse, and many questions are still to be answered, especially when including quantitative measures and more complex models.
Planning an incidence density sampled study requires information on established confounders, and these are not always present. Under incidence density sampling, matching on noninformative confounders may in some situations introduce data overfitting, and thereby bias on both the estimated IRR and on the standard error [25]. Even matching on an established confounder introduces bias, if this confounder is correlated with one of the unmatched confounders in the study.
When matching on age, the birthdays of all individuals, and the date or period of which the case become diagnosed are particularly important. This information may be hard to retrieve using surveys or other selfreporting sampling techniques. Hence, the incidence density sampled study is difficult to implement in countries without access to nationwide populationbased registers.
The simulations in this study are a simplified example of a GWAS only including one time scale and one genetic marker. Since risk of death increases with increasing age, age is an important factor to consider for matching. In other diseases, calendar time has a major effect due to changes in diagnostic systems, changes in politics etc. In such cases, one should consider which of the two timescales is the most important if the possibility of using both is limited. One should also consider other factors for matching that have a large impact on the risk of disease, e.g. gender, place of birth or work environment.
When selecting controls in the incidence density sampling, the controls should be chosen at random among persons still at risk at the time of disease onset of the case. The controls are a random sample of the population at that specific time point. Leaving out controls that later become cases or never experience other diseases, will cause bias to the estimated IRR [11, 26].
In comparison with the analogy of the tree, this study does not take into account the association between leaf type and the likelihood of falling onto the soil, among the leaves on the concrete. Including this association will make the interpretation of the estimated OR and the effect on the pvalues more complex, since this association will be reflected in the estimated OR. This association will however most likely be present in any population. For simplicity, we have also chosen to consider death as the only competing event. One should keep in mind that other events exist, such as emigration, which, for some age groups, are more likely than death [27]. The difference between the OR and the IRR is affected by an aggregation of all possible competing events. We urge the genetic community to carefully consider the specific research question and hence the sampling for any GWAS.
Conclusions
Competing events are universally related to human life precluding individuals from specific outcomes. For the identification of aetiological factors related to disease, interest is in the risk of becoming diseased. Using a classical casecontrol study to estimate the association between risk of becoming diseased and genetic variants may bias the association estimate and the pvalue of the association test. The incidence density sampled study will provide the correct interpretation of the estimated association and the proper ranking of SNPs for further analysis.
References
 1.
Stranger BE, Stahl EA, Raj T. Progress and promise of genomewide association studies for human complex trait genetics. Genetics. 2011;187(2):367–83.
 2.
Ligthart L, Hottenga JJ, Lewis CM, Farmer AE, Craig IW, Breen G, et al. Genetic risk score analysis indicates migraine with and without comorbid depression are genetically different disorders. Hum Genet. 2014;133(2):173–86.
 3.
Simonson MA, Wills AG, Keller MC, McQueen MB. Recent methods for polygenic analysis of genomewide data implicate an important effect of common variants on cardiovascular disease risk. BMC Med Genet. 2011;12:146.
 4.
Wray NR, Lee SH, Mehta D, Vinkhuyzen AA, Dudbridge F, Middeldorp CM. Research review: polygenic methods and their application to psychiatric traits. J Child Psychol Psychiatry. 2014;55(10):1068–87.
 5.
Dudbridge F. Power and predictive accuracy of polygenic risk scores. PLoS Genet. 2013;9(3), e1003348.
 6.
Andersen PK, Geskus RB, de Witte T, Putter H. Competing risks in epidemiology: possibilities and pitfalls. Int J Epidemiol. 2012;41(3):861–70.
 7.
Allignol A, Schumacher M, Wanner C, Drechsler C, Beyersmann J. Understanding competing risks: a simulation point of veiw. BMC Medical Research Methodology 2011, 11(86). doi:10.1186/147122881186
 8.
Pearce N. What does the odds ratio estimate in a casecontrol study. Int J Epidemiol. 1993;22(6):1189–92.
 9.
Clayton D, Hills M. Statistical models in epidemiology. New York: Oxford University Press Inc.; 1998.
 10.
Sistrom CL, Garvan CW. Proportions, odds, and risk. Radiology. 2004;230:12–9.
 11.
Wacholder S, Silverman DT, McLaughlin JK, Mandel JS. Seletion of controls in casecontrol studies. III. Design options. Am J Epidemiol. 1992;135(9):1042–50.
 12.
Prentice RL, Breslow NE. Retrospective studies and failure time models. Biometrika. 1978;65(1):153–8.
 13.
Beyersmann J, Latouche A, Buchholz A, Schumacher M. Simulating competing risks data in survival analysis. Stat Med. 2009;28(6):956–71.
 14.
Bender R, Augustin T, Blettner M. Generating survival times to simulate Cox proportional hazards models. Stat Med. 2005;24(11):1713–23.
 15.
Rosthøj S, Andersen PK, Abildstrom SZ. SAS macros for estimation of the cumulative incidence functions based on a Cox regression model for competing risks survival data. Comput Methods Programs Biomed. 2004;74(1):69–75.
 16.
Breslow NE, Day NE. Statistical Methods in Cancer Research. Volume I—The Analysis of CaseControl Studies. Lyon: International Agency for Research on Cancer (IARC Scientific Publications No. 32); 1980.
 17.
HoffmannJørgensen J. Probability With a View Towards Statistics, Volume 1, vol. 1. New York: Chapmann & Hall; 1994.
 18.
Fradin DD, Fallin MD. Influence of control selection in genomewide association studies: the example of diabetes in the Framingham Heart Study. BMC Preceedings. 2009;3(7):S113.
 19.
Wang MH, Shugart YY, Cole SR, Platz EA. A simulation study of control sampling methods for nested casecontrol studies of genetic and molecular biomarkers and prostate cancer progression. Cancer Epidemiol Biomarkers Prev. 2009;18(3):706–11.
 20.
Greenland S, Thomas DC. On the need for the rare disease assumption in casecontrol studies. Am J Epidemiol. 1982;116(3):547–53.
 21.
Fine JP, Gray RJ. A proportional hazards model for the subdistribution of a competing risk. J Am Stat Assoc. 1999;94(446):496–509.
 22.
Karon JM, Kupper LL. In defense of matching. Am J Epidemiol. 1982;116(5):852–66.
 23.
Kupper LL, Karon JM, Kleinbaum DG, Morgenstern H, Lewis DK. Matching in epidemiologic studies: validity and efficiency considerations. Biometrics. 1981;37(2):271–91.
 24.
Rose S, Laan MJ. Why match? Investigating matched casecontrol study designs with causal effect estimation. Int J Biostat. 2009;5(1):1.
 25.
Thomas DC, Greenland S. The relative efficiencies of matched and independent sample designs for casecontrol studies. J Chronic Dis. 1983;36(10):685–97.
 26.
Schwartz S, Susser E. Genomewide association studies: does only size matter? Am J Epidemiol. 2010;167(7):741–4.
 27.
Pedersen CB, Mortensen PB, CantorGraae E. Do risk factors for schizophrenia predispose to emigration? Schizophr Res. 2011;127(1–3):229–34.
Acknowledgement
This study is funded by The Lundbeck Foundation, though The Lundbeck Foundation Initiative for Integrative Psychiatric Research, iPSYCH as well as the Stanley Medical Research Institute, Bethesda, Maryland. None of the authors have any conflicts of interest.
Author information
Additional information
Authors’ contributions
BLW designed the study, developed the method, conducted the analysesand wrote the manuscript. AH, CBP and MN provided critical advice on methodology, result interpretation and manuscript preparation. All authors read and approved the final manuscript.
Additional files
Additional file 1:
The R code for simulations. (ZIP 2.06 kb)
Additional file 2:
Supplementary Figures S1–S10 with captions. (PDF 3.29 mb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Waltoft, B.L., Pedersen, C.B., Nyegaard, M. et al. The importance of distinguishing between the odds ratio and the incidence rate ratio in GWAS. BMC Med Genet 16, 71 (2015) doi:10.1186/s1288101502101
Received
Accepted
Published
DOI
Keywords
 Genome wide association study
 Study design
 Matched casecontrol study
 Competing risk
 Logistic regression
 Conditional logistic regression
 Rare disease assumption