<- See Stata 18's new features
Additive models of relative risk (RR)
Relative excess risk due to interaction (RERI)
Attributable proportion (AP)
Synergy index (SI)
Binomial generalized linear
Cox and other survival models
How do exposures interact to increase risk? Do you suspect that interaction is additive? Use reri to find out. Three measures of two-way interactions are provided: RERI, AP, and SI. Many models that estimate RR are supported, including logistic, binomial generalized linear, and survival.
The new reri command estimates additive interactions in binomial generalized linear models; logistic, Poisson, and negative binomial regressions; and Cox and other survival models.
Two-way interactions in models that estimate RR are typically expressed as multiplicative interactions. However, for many biological processes, an additive model for the interaction may model the process better than a multiplicative one (see, for example, Andersson et al. ).
For binary exposures A and B, we often want to model the risk of a positive outcome for subjects having both A and B compared with subjects having only A and those having only B. That is, we want to examine the interaction of A and B.
In a multiplicative model, the risk for subjects with A and B is hypothesized to equal the risk for only A times the risk for only B. If the risk for A and B equals the multiplicative risk, we say there is no interaction of A and B. If the risk is greater than the multiplicative risk, we say there is a positive interaction. Multiplicative models are convenient because they are simple to specify, and it is easy to test whether the interaction is greater (or less) than multiplicative.
In an additive model, the risk for subjects with A and B is hypothesized to equal the risk for only A plus the risk for only B. For this hypothesis, we want to look at how much the risk for A and B is greater (or less) than the additive risk.
The RERI statistic is this risk difference formulated using RRs (Rothman, Greenland, and Lash 2008). A RERI statistic of 0 means the risk is additive. A positive RERI statistic means the risk is superadditive, and a negative RERI statistic means the risk is subadditive. The test of statistical significance of the RERI statistic is straightforward, but it is not as simple as the test for the multiplicative model.
Besides the RERI statistic, the reri command reports two other related statistics: AP and SI. AP is the proportion of risk of A and B due to the superadditivity of exposures A and B. SI recasts the RERI statistic as a ratio rather than a difference.
We illustrate reri with an example that estimates the risk of birth defects from exposures to nitrates. We use a simulated dataset (nitrates.dta) that mimics the results from Brender et al. (2013).
. webuse nitrates (Simulated prenatal nitrate intake and birth defects data)
Our outcome measure is tube, representing neural tube defects. The exposures are two measures of nitrate intake: drug, with values 0/1, indicating nitrosatable drug exposure; and nitrate, also with values 0/1, representing either low or high daily consumption of nitrates from drinking water.
We want to fit an additive model for the interaction of drug and nitrate. Because this is an observational study, a binomial model of risk is appropriate, and we fit it using binreg:
. reri binreg tube drug nitrate Fitting binreg ... Interaction of drug and nitrate on an additive scale Model: Binomial generalized linear Number of observations = 1,284 drug#nitrate - + No#High + - Yes#Low + + Yes#High
|ERR std. err. z P>|z| [95% conf. interval]|
|- +||.1161419 .1571024 0.78 0.435 -.1529507 .4707206|
|+ -||.2385481 .2392316 1.11 0.268 -.1517982 .8085336|
|+ +||.7769424 .407042 2.51 0.012 .1342 1.783922|
|RERI||.4222523 .453286 0.93 0.352 -.466172 1.310677|
|Attr. prop.||.2376286 .2187062 1.09 0.277 -.1910276 .6662848|
|Synergy ind.||2.190483 2.034328 0.84 0.398 .3548285 13.52263|
The interactions created by reri represent disjoint groups in the data. The -+ interaction is for drug = No and nitrate = High. The +- and ++ interactions are defined similarly, as shown in the legend. The interaction -- is the reference category.
The estimates next to -+, +-, and ++ are measures of excess relative risk (ERR), defined as RR minus one. ERR greater than 0 means that RR is greater than 1.
If the model is additive, we would have
ERR(++) = ERR(-+) + ERR(+-)
The RERI statistic is the difference between the two sides of this equation:
RERI = ERR(++) - ERR(-+) - ERR(+-)
The estimate of ERR(++) is 0.777, which is greater than ERR(-+) + ERR(+-) = 0.116 + 0.239 = 0.355. So the model is superadditive. The RERI statistic is 0.422, the difference between 0.777 and 0.355.
The AP statistic is the proportion of risk for the ++ interaction that is due to the risk that is above additive. In our example, the estimate of AP is 0.24. The SI statistic is the excess risk expressed as a ratio rather than a difference. Its estimate of 2.19 is greater than 1, which implies a positive interaction or superadditivity.
As well as binreg models, reri can fit logistic, poisson, nbreg, stcox, stintcox, streg, and stintreg models.
Andersson, T., L. Alfredsson, H. Källberg, S. Zdravkovic, and A. Ahlbom. 2005. Calculating measures of biological interaction. European Journal of Epidemiology 20: 575–579.
Brender, J. D., P. J. Weyer, P. A. Romitti, B. P. Mohanty, M. U. Shinde, A. M. Vuong, J. R. Sharkey, et al. 2013. Prenatal nitrate intake from drinking water and selected birth defects in offspring of participants in the National Birth Defects Prevention Study. Environmental Health Perspectives 121: 1083–1089.
Rothman, K. J., S. Greenland, and T. L. Lash. 2008. Modern Epidemiology. 3rd ed. Philadelphia: Lippincott Williams & Wilkins.
Read more about the RERI, AP, and SI statistics estimated by the reri command in Stata Base Reference Manual; see [R] reri ->
View all the new features in Stata 18 ->