1 / 62

I nteraction and E ffect- M easure M odification

I nteraction and E ffect- M easure M odification. Lydia B. Zablotska, MD, PhD Associate Professor Department of Epidemiology and Biostatistics. Learning Objectives. Statistical interaction Multiplicative and additive interaction Biologic interaction

dale-schneider
Download Presentation

I nteraction and E ffect- M easure M odification

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Interaction and Effect-Measure Modification Lydia B. Zablotska, MD, PhD Associate Professor Department of Epidemiology and Biostatistics

  2. Learning Objectives • Statistical interaction • Multiplicative and additive interaction • Biologic interaction • Evaluation of interaction, presentation of results • Attributable fraction estimation

  3. Review of measures of association • Effect measures vs. measures of association: • Can never achieve counterfactual ideal • Logically impossible to observe the population under both conditions and to estimate true effect measures • Measures of association • Compares what happens in two distinct populations • Constructed to equal the effect measure of interest • Absolute: differences in occurrence measures (rate or risk difference) • Relative: ratios of occurrence measures (rate or risk ratio, relative risk, odds ratio)

  4. Comparison of absolute and relative effect measures Rothman 2002

  5. Concepts of interaction • Terms: • statistical interaction • effect modification or effect-measure modification • synergy (joint action of causal partners) • heterogeneity of effect • departure from additivity of effects on the chosen outcome scale • Definition: • heterogeneity of effect measures across strata of a third variable • Problems: • Scale-dependence, i.e. can be measured on an additive or multiplicative scale • Ambiguity of terms • Types: • Statistical • Biological • Public health interaction (public health costs or benefits from altering one factor must take into account the prevalence of other factors and effects of their reduction) RG Ch 5

  6. Types of interaction:Statistical interaction • If statistical interaction is being described on an additive scale then the measure of effect is the risk difference • R11 - R00 = (R10 - R00) + (R01 - R00). If the 2 sides of the equation are equal the relationship is perfectly additive • If statistical interaction is being described on a multiplicative scale then the measure of effect is the odds ratio or relative risk • R11 / R00 = (R10/R00 )(R01/R00). If the 2 sides of the equation are equal the relationship is perfectly multiplicative RG Ch 5

  7. Types of statistical interaction • Effect modification of the risk difference (absolute effect) corresponds with additive interaction • Effect modification on the risk ratio or odds ratio (relative effect) corresponds with multiplicative interaction • If there is no evidence of interaction on the multiplicative scale (i.e, heterogeneity of RR or OR if OR is a good approximation of RR) there will be evidence of interaction on the additive scale (i.e., heterogeneity of RD) https://www.ctspedia.org/do/view/CTSpedia/InterMultipl RG Ch 5

  8. Statistical interaction • Heterogeneity of effects always refers to a specific type of effect: risk ratios, odds ratios, risk differences • Absence of interaction for one measure does not imply absence of interaction for the other measures of association: • Homogeneity of risk differences implies heterogeneity of risk ratios and vice-versa • Most estimates of effect are based on multiplicative models; specify measures of effect when describing effect modification RG Ch 5

  9. Additive interaction RD = Riskexposed – Riskunexposed A and B are risk factors with risks Ra,- and R-,b and individual risk differences: RDa,- = Ra,- – R-,- RD-,b = R-,b – R-,- RDa,b is a RD for those exposed to both A and B and those exposed to neither • RDa,b = RDa,- + RD-,b – A and B are non-interacting risk factors • RDa,b  RDa,- + RD-,b – Additive interaction between A and B • RDa,b > RDa,- + RD-,b – Additive synergy (positive additive interaction) • RDa,b < RDa,- + RD-,b – Additive antagonism (negative additive interaction)

  10. Multiplicative interaction RR = Riskexposed / RiskunexposedRiskexposed = Riskunexposed x RR A and B are risk factors with risks Ra,- and R-,b and individual risk ratios: RRa,- = Ra,- / R-,- RR-,b = R-,b / R-,- RRa,b is a RR for those exposed to both A and B over those exposed to neither • RDa,b = RDa,- x RD-,b – A and B are non-interacting risk factors • RDa,b  RDa,- x RD-,b – Multiplicative interaction between A and B • RDa,b > RDa,- + RD-,b – Multiplicative synergy (positive multiplicative interaction) • RDa,b < RDa,- + RD-,b – Multiplicative antagonism (negative multiplicative interaction)

  11. Assessment of interaction for binary data

  12. Assessment of interaction for binary data Risk of past-year depression at age 26 according to genotype and stressful life events Dunedin Child-Development Study, Caspi et al. 2002, 2003

  13. Ways to assess interaction • By stratification • By comparing expected and observed joint effects • By graphical assessment • By estimating interaction magnitude • By using statistical tests of interaction: • Test of homogeneity • LRT comparing test and nested models

  14. Assessing interaction by stratification • Effect modification by presence of short allele G on the association between stressful life events E and risk of depression RDE/G is absent = 0.17-0.10=0.07; RRE/G is absent = 0.17/0.10=1.7 RDE/G is present = 0.33-0.10=0.23; RRE/G is present = 0.33/0.10=3.3 • Both RD and RR are heterogeneous: for both measures of association, the effect of E when G is present is larger than the effect of E when G is absent • Quantitative (ordinal) interaction – stratum-specific effects differ by magnitude but not by direction • Qualitative (disordinal) interaction – the same exposure could be protective in one stratum of the other factor but deleterious in the other

  15. Assessing interaction by comparing expected and observed joint effects • What is the individual effect of cause A in the absence of exposure to cause B? • What is the individual effect of cause B in the absence of exposure to cause A? • What is the observed joint effect of A and B? • What is the expected joint effect of A and B in the absence of interaction? • Is the observed joint effect similar to the expected joint effect in the absence of interaction?

  16. What is the individual effect of cause A in the absence of exposure to cause B? What is the individual effect of cause B in the absence of exposure to cause A? What is the observed joint effect of A and B? What is the expected joint effect of A and B in the absence of interaction? Is the observed joint effect similar to the expected joint effect in the absence of interaction? RDE,-=0.17-0.10=0.07 RD-,G=0.10-0.10=0 RDOBSERVED E,G=0.33-0.10=0.23 RDEXPECTED E,G=0.07+0=0.07 RDOBSERVED E,G > RDEXPECTED E,G, additive interaction (superadditivity, additive synergism) Comparing expected and observed joint effects

  17. What is the individual effect of cause A in the absence of exposure to cause B? What is the individual effect of cause A in the absence of exposure to cause A? What is the observed joint effect of A and B? What is the expected joint effect of A and B in the absence of interaction? Is the observed joint effect similar to the expected joint effect in the absence of interaction? What is the interaction magnitude RDE,-=0.17-0.10=0.07 RD-,G=0.10-0.10=0 RDOBSERVED E,G=0.33-0.10=0.23 RDEXPECTED E,G=0.07+0=0.07 RDOBSERVED E,G > RDEXPECTED E,G, additive interaction RDE/ G IS PRESENT – RDE/ G IS ABSENT = 0.23 - 0.07 =0.16 interaction contrast Comparing expected and observed joint effects RRE,-=0.17/0.10=1.7 RR-,G=010/0.10=1.0 RROBSERVED E,G=0.33/0.10=3.3 RREXPECTED E,G=1.7x1.0=1.7 RROBSERVED E,G > RREXPECTED E,G, multiplicative interaction RRE/G IS PRESENT / RRE/G IS ABSENT = 3.3 / 1.7 =1.9

  18. Assessing interaction by graphical representation Zablotska et al. 2013

  19. Assessing interaction by graphical representation Preston et al. 2007

  20. Trouble with assessment of synergy Interaction of vulnerability factors (e.g., fear of intimacy) and stressful life events in causing depression • Analysis on the additive scale: • Analysis on the multiplicative scale: Brown and Harris 1978

  21. Trouble with assessment of synergy Interaction of vulnerability factors (e.g., fear of intimacy) and stressful life events in causing depression • Analysis on the additive scale: 1. RDE,-=0.10-0.01=0.09; 2. RD-,M=0.03-0.01=0.02; 3. RDOBSERVED E,M=0.32-0.01=0.31; 4. RDEXPECTED E,M=0.09+0.02=0.11; 5. RDOBSERVED E,M>RDEXPECTED E,M; • additive interaction • Analysis on the multiplicative scale: RRE,-=0.10/0.01=10; 2. RR-,M=0.03/0.01=3; 3. RROBSERVED E,M=0.32/0.01=32; 4. RREXPECTED E,M=10x3=30; 5. RROBSERVED E,M~RREXPECTED E,M • no multiplicative interaction Brown and Harris 1978 Tenent and Bebbington 1978

  22. The conundrum • Each of these alternative interpretations is consistent with the premises of the mathematical models that were used: • Brown and Harris assumed that, absent interaction, risk factors add in their effects • Tennet and Bebbington assumed that, absent interaction, risk factors multiply in their effects • What is the answer and what could be done to elucidate one correct answer?

  23. Biological interaction • Terms: • Biological interaction • Causal interaction • Definition: • Modification of potential-response types • A process that explains potential mechanisms that can account for observed cases of disease • Exchangeability (i.e., the same data pattern would result if exposure status was switched or the rate in E would be equal to not E if E were not exposed) is required to test for interaction

  24. Biological interaction • Biological interaction can be defined under the counterfactual approach and the sufficient cause approach • Sufficient cause approach • 2 exposures are 2 component causes in a sufficient cause for the disease where the presence of both exposures is required to complete the sufficient cause ie., they are insufficient but necessary component causes of a unnecessary but sufficient cause (INUS partners) • interaction between component causes is implicit in the sufficient cause model • each component cause requires the presence of the others to act, their action is interdependent • Two possible types of interaction between component causes: • Synergism (biological interaction) • the two component causes work together as INUS partners in the same sufficient cause, they act in synergy. The individual would get disease if they are exposed to both A and B but not get disease if exposed to only A or B. • Parallelism • the two component causes compete to be INUS partners in the same sufficient cause, they act in parallel. The individual would get disease if they are exposed to either A and/or B but not get disease if exposed to neither. • synergy and parallelism have different component causes i.e, A and B, A or B.

  25. 4 possible types: parallelism A or B or AB U

  26. Biological vs. statistical interaction • When two factors have effects but risk ratios within the strata of the second factor are homogeneous, there is no interaction on the multiplicative scale • This implies that there is heterogeneity of the corresponding risk differences • The non-additivity of risk differences implies the presence of some type of biologic interaction RG Ch 5

  27. Interaction vs. effect modification Effect modification by Q of the effect of E on D without interaction between the effects of E and Q on D Potential interaction between the effects of E and Q on D without effect modification by Q of the effect of E on D E - drug in an RCT, D- hypertension, X - person’s genotype, and Q - person’s hair color. E - drug in an RCT for weight loss for obese children, D– weight at 6 mos, X – sugar intake, and Q – some measure of exercise. VanderWeele 2009

  28. Biological interaction, cont’ed • Biological interaction can be defined under the counterfactual approach and the sufficient cause approach • Counterfactual approach (potential outcome) • 4 exposure categories for 2 binary variables=16 possible patterns of response types (given disease or no disease) • 10 categories can be considered interaction (interdependence) of some type (i.e., both of the 2 exposure types have an effect) and interaction contrast not equal 0 • If it is assumed the effect is causal, Type 8 in the counterfactual approach is equivalent to causal or biological synergy. Each exposure only causes disease if the other is present.

  29. Possible response types for binary exposure

  30. Interaction contrast • Causal additivity = no causal interaction R11– R00= (R10 – R00) + (R01 – R00)=(p6+p13-p11-p13) + (p4+p11-p11-p13) =(0+0-0-0) + (0+0-0-0)=0 • Interaction contrast=difference in risk differences IC = RDX,-– RD-,Z = (R11 – R01)-(R10 – R00) = (R11 – R10)-(R01 – R00) = R11 – R10 – R01 + R00 = (p3+p5+2p7+p8+p15) – (p2+p9+2p10+p12+p14) RG Ch 5, p. 77

  31. Necessary conditions for interaction • Departures from additivity can only occur when interaction causal types are present in the cohort • Absence of interaction does not imply absence of interaction types because sometimes different interaction types counterbalance each other’s effect on the average risk • Definitions of response types depend on the definition of the outcome under study (if it changes, then response type can change too) RG Ch 5

  32. Departures from additivity • Superadditivity: RD11>RD10+RD01 – type 8 MUST be present • Subadditivity: RD11<RD10+RD01 – type 2 MUST be present • However, presence of synergistic responders (type 8) or competitive responders (type 2) does not imply departures from additivity • If neither factor is ever preventive: IC = p8 –p2, • i.e. synergism – parallelism = additive interaction

  33. This is all good, but how do we know the response types? 16 1 6 8 R R R R

  34. Simplified assessment of synergy based on 5 response types • p8 = (R11 – R01) – (R10 – R00) • Effect of Z (effect modifier) when X=1 – Effect of Z when X=0 • Assumptions when only 5 types are used • Effect measure is the Risk Difference, biologic interaction is then interaction for risk differences • p5 > 0, biologic interaction must be positive (although one can reparameterise the exposures X and Z to get a negative interaction) • Huge reduction of person types, from 16 to 5! • Keep in mind that this is a "biologic“ model

  35. Summary of R&G scheme under counterfactual theory • The reduction from 16 person types to 5 makes it possible to get the p’s for the 5 types, by using the 4 observed probabilities, and the fact that the 4 R’s sum to 1. • By solving the equations we get that the person type “synergy” is equal to additive interaction, with risk differences as measure of effect

  36. Critique of R&G scheme • Rothman and Greenland's model is simplistic. • One reasonable person type is missing! p2 -Parallelism • If A and B are both causal, then it is reasonable to think that some individuals in the population will develop the disease when exposed to only A, only B or both A and B.

  37. Darroch, J. “Biologic Synergism and Parallelism”, AmJEpi 1997; 145:7 page 661-668 • John Darroch discusses an expansion of the ideas by Rothman and Greenland. He assumes 6 person types, including "parallelism". • By using 6 person types he covers all the possible person types if A and B are directly causal in their effect on disease.

  38. 16 1 6 8 2 R R R R

  39. Simplified assessment of synergy based on 6 response types • p8 – p2 = (R11 – R01) – (R10 – R00) • Effect of Z (effect modifier) when X=1 – Effect of Z when X=0 • This means you will not be able to specify the biologic interaction (p8) exactly from the 4 known probabilities, but you can find the boundaries.

  40. Summary notes on synergy and parallelism • Can only be partially determined from the data at hand • Example of synergy (assuming the factors are causal ): if the gene and environment factors acted together, infants would only get the congenital disorder if exposed to both gene and environment • Example of parallelism (assuming the factors are causal ): infants would only get the congenital disorder if exposed to either gene or environment but would not get the congenital disorder if exposed to neither. • If synergy - parallelism or R(AB) - R(AB) - R(A) - R(B) + R is a positive number the result is consistent with the presence of more synergy than parallelism in the population studied • The public health approach would be to prevent exposure to either genes or environment • Greater than an additive relationship is consistent with superadditivity and multiplicativity but inconsistent with the single hit model of disease causation • If synergy – parallelism or R(AB) - R(A) - R(B) + R is a negative number it is an indication that there is more parallelism than synergy in the population • Less than an additive relationship is consistent with subaddivitity and inconsistent with the no hit and multistage models of disease • The public health approach would be to prevent exposure to both genes and environment. • If there is no additive interaction there may be no synergism or the proportion of individuals for whom the exposures work synergistically may be the same for whom the exposures work in a parallel manner

  41. Example from Darroch 1997

  42. Darroch vs. R&G • p8 = (R11 – R01) – (R10 – R00) R R 8 R R 2 6 R R R R R R

  43. Darroch vs. R&G p8 = (20.7 – 5.1) – (7.2 – 1) = 9.4 > 0 - superadditivity p8 = (R11 – R01) – (R10 – R00) R R R R R R 8 R R 6 2 R R R R R R 8

  44. An additive model with a “twist” • Additive model with a “twist” allows the best representation of synergy • An additive model assumes that risks add in their effects • Positive deviations from additivity (superadditivity) indicates the presence of synergy • The “twist” is that risks do something slightly less than add (parallelism – some individuals can develop disease from either one of the two exposures under study) • What we see as the combined effect of two exposures reflects the balance of synergy and parallelism • In summary, although superadditivity indicates synergy, a failure to find superadditivity does not imply the absence of synergy

  45. Estimating synergy in epi studies • If there is positive interaction on the multiplicative scale, there will be positive interaction on the additive scale (supermultiplicativity implies superadditivity) • We can assess interaction on the additive scale from the multiplicative model by calculating: • Interaction contrast and interaction contrast ratio • Relative excess risk due to interaction (RERI) – for multivariate models or models with continuous exposures • Attributable proportion • Synergy index

  46. In a 2x2 table • IC=0.33-0.17-0.10+0.10=0.16 >0  synergy Dunedin Child-Development Study, Caspi et al. 2002, 2003

  47. Estimation of IC and ICR • Cohort studies • Intercept provides the baseline odds of disease • OR for risk factors could be used to obtain the odds of disease under the other conditions • Odds could be converted to risks (odds=p/ (1-p)) • Case-controls studies • Intercept may be biased • Odds for those exposed to both factors: 0.33/0.67; odds for those exposed to life events only: 0.17/0.83; odds for those with short allele only: 0.10/0.90; odds for those exposed to neither: 0.10/0.90 • ICR=ORboth/neither-ORlife events/neither-ORshort allele/neither + baseline ICR=((0.33/0.67)/(0.10/0.90)) –((0.17/0.83)/(0.10/0.90)) – –((0.10/0.90)/(0.10/0.90)) +1=2.6 ICR/ORboth/neither=2.6/4.4=0.59 – the proportion of disease among those with both risk factors that is attributable to interaction RG Ch 16

  48. Bringing it all together: From synergy to its mathematical representation Brown and Harris 1978

  49. Causes of depression: Theory about life events and their interaction with intimacy problems

More Related