Why the Response Scale Matters
When two drugs are studied together, the combination response can differ from what would be expected by combining their separate effects. For a continuous outcome, that comparison can be expressed as a factorial interaction contrast. The contrast is not a scale-free property of the drugs: it asks whether effects combine additively on a specified response scale. A conclusion on the original scale may therefore differ from a conclusion on the log scale.
Consider a two-by-two factorial design. Drug A is either absent or present, as is drug B. Let \(\mu_{ab}\) be the mean outcome for the cell with A status \(a\) and B status \(b\), where \(a,b\in\{0,1\}\). The outcome could be a continuous biomarker, symptom score, or another measured endpoint. The four means must refer to the same outcome definition, measurement time, and target population.
The sign has a direct interpretation only after defining what larger outcomes mean. If larger values indicate greater benefit, a positive contrast means the combination exceeds the additive prediction; if larger values indicate greater harm, its clinical meaning reverses. A nonzero contrast is evidence of nonadditivity on the chosen scale, not by itself proof of a biological mechanism or a universal claim of synergy.
From the Four Cell Means to a Product Term
A linear model with binary indicators \(a\) and \(b\) represents the four cell means as
The product term \(ab\) is zero in three cells and one only when both drugs are present. Substituting each pair of indicator values gives \(\mu_{00}=\beta_0\), \(\mu_{10}=\beta_0+\beta_A\), \(\mu_{01}=\beta_0+\beta_B\), and \(\mu_{11}=\beta_0+\beta_A+\beta_B+\beta_{AB}\). Subtracting the cell means in the interaction contrast yields
Thus the product coefficient is exactly the factorial interaction contrast when the model uses 0/1 indicators and is saturated for the four cell means. The main-effect coefficients are conditional: \(\beta_A\) is A's mean effect when B is absent, and \(\beta_B\) is B's mean effect when A is absent. They are not necessarily average effects across both levels of the other drug.
Worked Example: Estimating an Additive Interaction
Worked Example: Estimating an Additive Interaction
In a hypothetical randomized factorial study, the outcome is a continuous improvement score in points, with higher values indicating greater improvement. Each of the four treatment cells has 25 participants. The cell means are 10 points for neither drug, 14 for A alone, 13 for B alone, and 20 for the combination. Suppose each cell has sample standard deviation 4 points.
The interaction contrast is
The additive prediction for the combination is \(14+13-10=17\) points: start from the 10-point mean with neither drug, then add A's 4-point effect and B's 3-point effect. The observed combination mean is 20, which is 3 points above that prediction. Equivalently, A's effect is \(20-13=7\) points with B present and \(14-10=4\) points with B absent; their difference is 3 points.
For a conventional pooled-variance analysis, assume independent observations, approximately normal within-cell errors, and a common within-cell variance. The four observed standard deviations are all 4 points, consistent with the common-variance assumption for this example. The pooled variance estimate is \(s_p^2=16\) square points. With equal cell sizes, the contrast standard error is
To test \(H_0:\Delta_Y=0\) against \(H_a:\Delta_Y\ne 0\), the statistic is \(t=3/1.60=1.875\). With \(df=4(25-1)=96\), the two-sided p-value is approximately 0.064. At \(\alpha=0.05\), this does not provide sufficient evidence to reject zero interaction on the additive points scale. A 95% confidence interval is \(3\pm t_{0.975,96}(1.60)\), or \(3\pm1.985(1.60)=(-0.176,\ 6.18)\) points. Under the stated model, this interval is compatible with effects ranging from a small negative contrast to a positive contrast of about 6.18 points; it is not a 95% probability statement about the fixed parameter.
In an unbalanced study, the same contrast of estimated cell means remains the estimand, but its standard error must reflect the cell-specific sample sizes and variance estimates. A model-based standard error can be used if its variance assumptions are suitable; robust or design-based methods may be preferable when they are not. Randomization supports a causal interpretation of cell-mean differences for the assigned treatments when treatment implementation, outcome measurement, and follow-up preserve the intended comparison. Missing outcomes, nonadherence, or interference between participants can complicate that interpretation.
Interaction on an Additive or Multiplicative Scale
For a positive outcome, a common alternative is to model its logarithm. Let \(\nu_{ab}=E(\log Y\mid a,b)\), and define the interaction contrast on that transformed scale as \(\Delta_{\log Y}=\nu_{11}-\nu_{10}-\nu_{01}+\nu_{00}\). Exponentiating gives a ratio of ratios:
Here \(\exp(\nu_{ab})\) is the cell's geometric mean, when that quantity is defined. A ratio of ratios equal to 1 corresponds to no interaction on the log scale: the relative effect of A is the same with and without B. This is different from a zero additive contrast in the original outcome units. A log model requires positive outcomes unless a different transformation or model is justified; adding an arbitrary constant to accommodate zero values changes the estimand and needs substantive justification.
Worked Example: The Interaction Conclusion Changes With Scale
Suppose hypothetical positive biomarker outcomes have cell means of 1, 10, 10, and 50 units for neither drug, A alone, B alone, and the combination, respectively. For illustration, suppose the same four numbers also represent the corresponding geometric means.
On the original scale, the factorial contrast is
The contrast is positive: the combination mean is 31 units above the additive prediction \(10+10-1=19\) units. On the log scale, however,
The ratio of ratios is \(0.5\), not 1. On this scale, A's relative effect is smaller when B is present: the A-to-no-A geometric-mean ratio is \(50/10=5\) with B and \(10/1=10\) without B. Thus the additive contrast is positive while the log-scale contrast is negative. These are not contradictory calculations: they answer different questions about the four means.
Affine transformations behave differently from nonlinear ones. If \(Y^*=c+dY\), then
The constant cancels, and a positive change of units preserves the sign and zero status of the contrast. A nonlinear transformation such as a logarithm does not generally preserve additivity, so both the magnitude and whether the contrast is zero can change. The scale should therefore be chosen to match the scientific question and the outcome's measurement properties, rather than selected only because it produces a preferred interaction result.
Worked Example: Recovering the Contrast From a Regression
Worked Example: Recovering the Contrast From a Regression
Use the improvement-score means from the first example: \(\mu_{00}=10\), \(\mu_{10}=14\), \(\mu_{01}=13\), and \(\mu_{11}=20\) points. Under the 0/1 indicator model, the coefficients are obtained from the cell means:
The combination mean then determines the product coefficient:
This equals the direct factorial contrast, \(20-14-13+10=3\) points. If the indicator coding or reference levels change, individual regression coefficients can change their interpretation, but the four-cell contrast can still be reconstructed from the fitted cell means. In models that impose additional restrictions, such as omitting the product term, the fitted cell means need not equal the observed cell means, and the coefficient no longer represents the unrestricted four-cell contrast.
Assumptions, Diagnostics, and Alternatives
The contrast is defined by four means and can be calculated without assuming normality. Statistical uncertainty and tests require an appropriate sampling or modeling framework. For the pooled \(t\) analysis illustrated above, observations should be independent, errors should be reasonably normal within cells for small-sample inference, and the cell variances should be sufficiently similar. With larger samples, inference on means is often less sensitive to moderate nonnormality, but severe outliers, dependence, or strong variance differences can still undermine a simple pooled analysis.
Check the design and data before interpreting the result. Confirm that all four treatment combinations are represented, the same endpoint and assessment time are used in every cell, and the outcome coding has the intended direction. Inspect cell-specific distributions, residuals, and variances; look for influential observations, ceiling or floor effects, and uneven missingness. A missing or very small cell makes the contrast unstable or impossible to estimate precisely. A model can return a number despite poor support in the data, so successful computation is not evidence that assumptions hold.
The contrast is useful because it is transparent: it compares observed combination response with an additive prediction on a stated scale. Simple-effect comparisons can show which conditional drug effect differs, but several separate tests do not replace a direct test of the prespecified interaction contrast. A log-scale contrast instead targets multiplicative combination of geometric means. Bliss, Loewe, and Highest-Single-Agent additivity use different reference expectations for combination response; their conclusions should not be treated as interchangeable with a factorial contrast on raw or log-transformed continuous outcomes. None of these comparisons alone establishes a molecular mechanism.
A four-cell contrast also has limits. It summarizes interaction at the particular doses, population, endpoint, and time studied; it does not describe how interaction varies across dose combinations. That broader dose-response question requires a model appropriate to the design and is a separate modeling task. Finally, an estimated contrast near zero does not prove exact additivity: its interval may still include effects large enough to matter clinically. Precision and a clinically meaningful threshold should be considered alongside a p-value.
Common Mistake
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- For cell means \(\mu_{00}=6\), \(\mu_{10}=9\), \(\mu_{01}=8\), and \(\mu_{11}=13\) units, calculate the additive interaction contrast and the additive prediction for the combination.
- In the model \(E(Y\mid a,b)=\beta_0+\beta_Aa+\beta_Bb+\beta_{AB}ab\), what does \(\beta_A\) mean, and when does it equal the factorial interaction contrast?
- Four positive cell geometric means are 4, 8, 6, and 15, in the order \(00,10,01,11\). Calculate the log-scale interaction contrast and ratio of ratios.
- Why can a zero interaction contrast on the raw scale differ from a zero contrast on the log scale?
Solutions. (1) \(\Delta_Y=13-9-8+6=2\) units. The additive prediction is \(9+8-6=11\) units, so the observed combination mean is 2 units above it. (2) \(\beta_A\) is the effect of A when B is absent, since it compares \(a=1\) with \(a=0\) at \(b=0\). In a saturated 0/1 two-factor model, \(\beta_{AB}\), not generally \(\beta_A\), equals the factorial interaction contrast. (3) \(\Delta_{\log Y}=\log(15)-\log(8)-\log(6)+\log(4)=\log(60/48)=\log(1.25)\approx0.2231\); the ratio of ratios is \(1.25\). (4) The raw contrast measures departures from additivity in outcome units, whereas the log contrast measures departures from multiplicativity in geometric means. A nonlinear transformation does not generally preserve the zero-contrast condition.