Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Interaction Contrasts on Continuous Response Scales

Drug-Drug Interaction PK/PD · Tutorial 303 of 598

Interaction Contrasts on Continuous Response Scales

Learn to calculate and interpret continuous-response interaction contrasts, connect them to product terms, and assess how changing the outcome scale changes the estimand.

Advanced 11 min read

What You'll Learn

  • Define the four-cell factorial interaction contrast for a continuous response.
  • Derive its equivalence to the product coefficient in a two-factor regression.
  • Calculate an interaction contrast, standard error, test statistic, and confidence interval.
  • Interpret interaction on additive and log response scales.
  • Explain why nonlinear outcome transformations can change the interaction conclusion.
  • Distinguish a response-scale interaction contrast from other additivity reference models.

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.

Definition: On the scale of the outcome \(Y\), the factorial interaction contrast is $$ \Delta_Y=\mu_{11}-\mu_{10}-\mu_{01}+\mu_{00}. $$ Equivalently, \(\Delta_Y=(\mu_{11}-\mu_{01})-(\mu_{10}-\mu_{00})\): the effect of A when B is present minus the effect of A when B is absent. A zero contrast means no departure from additivity of the four cell means on this particular scale.

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

$$ E(Y\mid a,b)=\beta_0+\beta_A a+\beta_B b+\beta_{AB}ab. $$

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

$$ \Delta_Y=(\beta_0+\beta_A+\beta_B+\beta_{AB}) -(\beta_0+\beta_A)-(\beta_0+\beta_B)+\beta_0 =\beta_{AB}. $$

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.

Justification: Each cell mean enters \(\Delta_Y\) with coefficients \(+1,-1,-1,+1\). Under the saturated model, the intercept and the two single-drug terms cancel; only the coefficient multiplying \(ab\) remains. This cancellation also shows that adding a constant to every outcome does not change the contrast.

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

$$ \widehat{\Delta}_Y=\bar{Y}_{11}-\bar{Y}_{10}-\bar{Y}_{01}+\bar{Y}_{00} =20-14-13+10=3\text{ points}. $$

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

$$ SE(\widehat{\Delta}_Y) =\sqrt{s_p^2\left(\frac{1}{25}+\frac{1}{25}+\frac{1}{25}+\frac{1}{25}\right)} =\sqrt{16(0.16)}=1.60\text{ points}. $$

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:

$$ \exp(\Delta_{\log Y}) =\frac{\exp(\nu_{11})/\exp(\nu_{01})}{\exp(\nu_{10})/\exp(\nu_{00})}. $$

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

$$ \Delta_Y=50-10-10+1=31\text{ units}. $$

The contrast is positive: the combination mean is 31 units above the additive prediction \(10+10-1=19\) units. On the log scale, however,

$$ \Delta_{\log Y}=\log(50)-\log(10)-\log(10)+\log(1) =\log\left(\frac{50(1)}{10(10)}\right) =\log(0.5)\approx-0.6931. $$

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

$$ \Delta_{Y^*} =(c+d\mu_{11})-(c+d\mu_{10})-(c+d\mu_{01})+(c+d\mu_{00}) =d\Delta_Y. $$

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:

$$ \beta_0=\mu_{00}=10,\qquad \beta_A=\mu_{10}-\mu_{00}=4,\qquad \beta_B=\mu_{01}-\mu_{00}=3. $$

The combination mean then determines the product coefficient:

$$ \beta_{AB}=\mu_{11}-\beta_0-\beta_A-\beta_B =20-10-4-3=3\text{ points}. $$

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

Pitfall: “The combination mean is greater than either single-drug mean, so the drugs interact.” This compares the combination with each drug separately, not with the additive prediction. The correction is to calculate \(\mu_{11}-\mu_{10}-\mu_{01}+\mu_{00}\) on the chosen scale, quantify uncertainty, and state which outcome scale gives the contrast its interpretation. A combination can outperform both single drugs and still have a zero additive interaction contrast.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. 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.
  2. 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?
  3. 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.
  4. 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.

Clinical Biostats