Tutorials › Causal Inference › Prespecifying Subgroups in Trial Protocols

Effect Modification and Subgroup Effects · Tutorial 867 of 972

Prespecifying Subgroups in Trial Protocols

Build a prespecified subgroup analysis plan that makes the target comparisons and their interpretation clear before trial outcomes are available.

Advanced 11 min read

What You'll Learn

  • Define subgroups using baseline variables, measurement rules, and prespecified cut points.
  • State a causal subgroup estimand and its target population.
  • Choose an effect scale that answers the clinical question.
  • Align the estimation model with the estimand and specify it in advance.
  • Distinguish confirmatory subgroup questions from exploratory analyses.
  • Recognize how outcome access and post hoc choices can undermine interpretation.

Why a Subgroup Plan Must Come Before Outcomes

A trial may show an average treatment effect while clinicians still need to know whether the effect differs across patient groups. That question is easy to phrase after results are visible: try a different age threshold, inspect another biomarker, or report whichever subgroup looks most striking. The flexibility is the problem. When subgroup definitions or analysis choices are selected in response to the observed outcomes, the resulting estimates and tests no longer have the same interpretation as a plan fixed in advance.

A protocol-ready plan is more than a list of variables. It specifies who belongs to each subgroup, what causal contrast will be estimated, on what scale, how it will be estimated, and which questions are confirmatory rather than exploratory. It also records decisions that could otherwise be made after seeing the results, such as cut points, category handling, and the primary outcome scale.

Definition: A prespecified subgroup analysis is an analysis whose subgroup variable and definition, target contrast, effect scale, and principal analysis method are documented before trial outcomes are examined. The subgroup characteristic should ordinarily be measured before randomization. A confirmatory analysis addresses a stated, prioritized scientific question under a prespecified analysis plan; an exploratory analysis is intended to generate or investigate hypotheses and is interpreted with appropriate caution.

Specify the Subgroup and Its Target Contrast

Begin by naming the subgroup variable precisely. State how it is measured, when it is measured, its units, and how it is coded. For a categorical variable, list the categories and their reference or comparison group. For a continuous variable, say whether it remains continuous or is divided into categories. If a cut point is used, give its numerical value, units, rationale, and rule for observations exactly at the boundary. A label such as “older patients” is not a reproducible definition.

Let \(A\) denote randomized treatment assignment, with \(A=1\) for the intervention and \(A=0\) for control. Let \(S\) be a baseline subgroup characteristic, and let \(Y^a\) denote the outcome that would be observed under assignment \(a\). For a subgroup defined by \(S=s\), a common causal estimand is the conditional average treatment effect:

$$ \tau(s)=E[Y^1-Y^0\mid S=s]. $$

This expression compares the mean potential outcome under intervention with the mean potential outcome under control among people in the specified subgroup. For a binary outcome, the plan must say whether the contrast is, for example, a risk difference or a risk ratio. The estimand should also state the outcome and follow-up time, the treatment conditions being compared, and the target population. If the analysis concerns treatment assignment regardless of adherence, state that intention-to-treat target explicitly.

A subgroup effect is not itself evidence that effects differ across subgroups. To describe additive effect modification between two groups, the relevant contrast is the difference between their treatment effects, such as \(\tau(s_1)-\tau(s_0)\). The plan should distinguish this contrast from simply reporting a treatment effect separately in each group. It should also name the scale: effect modification on the risk-difference scale need not agree with effect modification on the risk-ratio or odds-ratio scale.

Identification and Estimation in a Randomized Trial

For a baseline subgroup in a randomized trial, randomization supports a causal comparison within that subgroup. In precise terms, if assignment is independent of the potential outcomes conditional on \(S\), both treatment arms have positive probability within each target subgroup, and the observed outcome under the assigned treatment equals the corresponding potential outcome, then:

$$ E[Y^a\mid S=s]=E[Y\mid A=a,S=s]. $$

The assumptions are conditional exchangeability from randomization, positivity, and consistency. Positivity requires that each treatment being compared can occur for members of the subgroup. Consistency links the outcome actually observed to the potential outcome under the treatment actually assigned. Randomization may be stratified on the subgroup variable, but stratification does not replace the need to specify the estimand and analysis.

The equality above justifies estimating subgroup-specific causal effects from the randomized-arm outcomes, subject to the stated assumptions and adequate outcome measurement. It does not guarantee that a particular model is correctly specified, that estimates will be precise, or that missing outcomes are harmless. A protocol should describe how the estimand will be estimated and how uncertainty will be quantified. For a categorical subgroup and a binary outcome, for example, one approach is to estimate each arm’s risk within each subgroup and subtract the risks to estimate a risk difference. A regression model with treatment-by-subgroup terms is another possible approach, but the model and the scale of the reported contrasts must match the target estimand.

Worked Example: Prespecifying an Age Subgroup on the Risk-Difference Scale

Consider a hypothetical randomized trial of a treatment to prevent a binary complication by 30 days. Before outcomes are available, the protocol defines age as age in years at randomization and sets the subgroups to under 65 years and 65 years or older. The primary subgroup estimand is the intervention-minus-control difference in 30-day complication risk in each age group. The prespecified interaction contrast is the difference between those two risk differences.

Suppose the hypothetical results are 12 complications among 100 treated participants and 20 among 100 control participants under age 65; among participants aged 65 or older, there are 24 among 100 treated and 30 among 100 control participants. The under-65 risk difference is:

$$ \widehat{RD}_{<65}=\frac{12}{100}-\frac{20}{100} =0.12-0.20=-0.08. $$

Thus the estimated intervention risk is 8 percentage points lower than the control risk in the under-65 subgroup. For the older subgroup:

$$ \widehat{RD}_{\geq65}=\frac{24}{100}-\frac{30}{100} =0.24-0.30=-0.06. $$

The estimated additive interaction contrast, defined as the under-65 risk difference minus the 65-or-older risk difference, is \(-0.08-(-0.06)=-0.02\), or a 2-percentage-point more negative risk difference in the younger group. These are descriptive point estimates from hypothetical data; whether the contrast is estimated precisely requires its uncertainty interval. The protocol would additionally specify the estimator, interval method, and whether this contrast is confirmatory. The age cut point and direction of the contrast are fixed independently of these outcomes.

Choose a Scale and an Analysis Model in Advance

The scale is part of the scientific question, not a cosmetic reporting choice. An absolute difference in event risks answers how many percentage points the risks differ. A risk ratio answers how the risks compare proportionally. An odds ratio compares odds and is not generally interchangeable with a risk ratio. The protocol should identify one primary scale for each confirmatory subgroup estimand and explain why it is useful for the clinical decision. Other scales may be reported as secondary or descriptive, but they should not quietly replace the prespecified primary scale when results are known.

The model should be chosen to estimate the stated contrast. For a continuous outcome, a linear model may directly target conditional mean differences if its mean structure is appropriate. For a binary outcome, a model for the outcome probability or a model followed by marginal standardization can be used to obtain risks and risk differences. A model coefficient is not automatically the desired causal contrast: for example, a logistic-regression interaction coefficient is on the log-odds scale. If the target is a risk difference, the plan should explain how predicted risks will be converted to that scale.

For continuous effect modifiers, prespecify the functional form rather than searching for a convenient split. If the plan uses a linear treatment-by-marker interaction, state the marker’s units, centering value, and the range over which the interaction will be interpreted. If a non-linear relationship is scientifically plausible, specify the planned form or a limited sensitivity analysis in advance. Do not choose the cut point or curve after comparing which version produces the strongest apparent treatment difference.

Worked Example: Keeping a Biomarker Continuous

A hypothetical trial measures a baseline biomarker in units per litre and uses a continuous symptom score, where lower scores indicate fewer symptoms. The protocol retains the biomarker continuously, centers it at 50 units per litre, and reports changes in treatment effect per 10-unit increase. A prespecified linear mean model is:

$$ E[Y\mid A,M]=\beta_0+\beta_A A+\beta_M\left(\frac{M-50}{10}\right) +\beta_{AM}A\left(\frac{M-50}{10}\right). $$

Here \(M\) is the baseline biomarker value, and \(A=1\) denotes intervention. Under this model, the intervention-minus-control mean difference at marker value \(M\) is \(\beta_A+\beta_{AM}(M-50)/10\). Suppose, purely hypothetically, the fitted coefficients are \(\hat{\beta}_A=-2.0\) score points and \(\hat{\beta}_{AM}=-0.6\) score points per 10-unit increase. At \(M=40\), the centered and scaled marker is \((40-50)/10=-1\), giving:

$$ -2.0+(-0.6)(-1)=-1.4 \text{ score points}. $$

At \(M=60\), it is \(1\), giving \(-2.0+(-0.6)(1)=-2.6\) score points. Under the stated model, the estimated mean score reduction is larger at 60 than at 40 units per litre. These calculations illustrate the model’s implied contrasts; they do not establish that the linear interaction is adequate. The protocol should state that the interaction is on the additive mean-score scale and specify how model fit and uncertainty will be assessed.

Separate Confirmatory Questions From Exploratory Work

A useful plan identifies a small set of prioritized subgroup questions and labels each one. For every confirmatory question, record the subgroup definition, outcome and time point, estimand and scale, model or estimator, and the planned uncertainty summary. State the scientific rationale and whether the question concerns an effect within a subgroup or a contrast of effects between subgroups. A large subgroup estimate in one group and a small estimate in another do not, by themselves, show that the effects differ; the contrast between effects is the relevant quantity.

Exploratory analyses may examine additional baseline characteristics or alternative functional forms, but their status should be explicit. Report how those subgroups were defined and identify analyses that were not part of the confirmatory plan. Avoid presenting a data-driven threshold as though it had been chosen in advance. Prespecification improves transparency; it does not turn an underpowered or imprecise analysis into strong evidence.

The protocol or statistical analysis plan should be dated and finalized before access to unblinded outcome comparisons. Where feasible, analysts can work with treatment labels masked while finalizing data rules and code. The plan should also cover missing subgroup measurements, values outside the expected range, and missing outcomes. State whether the analysis will use complete cases, a specified imputation approach, or another method, and identify any assumptions required. The later interpretation depends on these choices just as much as on the subgroup label.

Worked Example: Prespecifying the Primary Scale and a Secondary Scale

A hypothetical randomized trial studies a treatment for a binary 90-day event and prespecifies baseline kidney function below 60 versus 60 or above, measured in the study’s stated units. The protocol names the risk difference as the primary scale because the investigators want to describe absolute event reduction. The risk ratio is listed as a secondary descriptive scale, not as an alternative to be selected after results are seen.

Suppose the hypothetical below-60 subgroup has 18 events among 120 treated participants and 24 among 120 controls. Its risk difference is \(18/120-24/120=0.15-0.20=-0.05\), a 5-percentage-point lower intervention risk. The risk ratio is \((18/120)/(24/120)=0.15/0.20=0.75\).

Suppose the 60-or-above subgroup has 9 events among 180 treated participants and 18 among 180 controls. Its risk difference is \(9/180-18/180=0.05-0.10=-0.05\), again a 5-percentage-point lower intervention risk. Its risk ratio is \((9/180)/(18/180)=0.05/0.10=0.50\). The point estimates therefore agree on the absolute scale but differ on the relative scale. This reinforces why the protocol must name the scale for its target question. It does not, by itself, show that the effects differ statistically or clinically; that requires a prespecified comparison and an assessment of its uncertainty.

Common Pitfalls and a Protocol Checklist

The most common mistake is to treat prespecification as a label rather than a set of reproducible decisions. Writing “subgroups will be examined” does not define what will be estimated or how. Likewise, a cutoff chosen after viewing outcome patterns is not made confirmatory merely by being described in the final report as clinically plausible. Keep a dated record of the plan and disclose deviations.

  • Vague subgroup definitions. Replace “high risk” with a measured baseline variable, units, threshold, and boundary rule.
  • Post-outcome cut-point searching. Do not try multiple thresholds and report only the most favorable one as the planned analysis.
  • Mixing scales. State whether the primary effect is an absolute difference, relative effect, or another contrast, and interpret estimates on that scale.
  • Confusing subgroup estimates with heterogeneity. Prespecify the contrast between subgroup effects if the question is whether the effect differs across groups.
  • Overstating exploratory findings. Mark unplanned analyses as exploratory and avoid treating their apparent patterns as confirmation.

A complete protocol entry can be checked by asking: Is the baseline subgroup variable and its coding reproducible? Is the target population and outcome time point explicit? Are the treatment contrast, estimand, and effect scale stated? Does the model or estimator target that estimand? Are confirmatory and exploratory questions distinguished? Are missing data and deviations addressed? Prespecification cannot eliminate uncertainty, but it makes the analysis auditable and limits opportunities for outcome-driven choices. The detailed plan for handling the combined evidence from multiple subgroups and outcomes should also be stated rather than left implicit.

Key takeaway: Prespecify the subgroup definition, causal estimand, effect scale, estimator, and confirmatory status before outcomes are examined. A reproducible plan makes clear both what the trial can answer and how cautiously its subgroup findings should be interpreted.

Check Your Understanding

Use these questions to test whether a subgroup analysis plan is precise enough to implement.

  1. What details should a protocol give when it defines a subgroup using a continuous baseline characteristic?
  2. For a binary outcome, what is the difference in interpretation between a risk difference and a risk ratio?
  3. Why is a treatment effect estimated in one subgroup not, by itself, evidence that treatment effects differ between subgroups?
  4. A protocol says it will examine “patients with high biomarker values.” What information is missing for a reproducible subgroup definition?
  5. What makes an analysis exploratory, and how should that status affect how its findings are described?
Clinical Biostats