Why the Reason a Biomarker Is Unavailable Matters
A biomarker measurement can be absent for different reasons. A scheduled blood sample may have been missed, a specimen may have been lost, or a participant may have withdrawn before the visit. In these cases the value may still be conceptually defined at the scheduled time, even though it was not recorded. By contrast, a biomarker value at week 12 may not be defined for a participant who died at week 8. The issue is not merely that the measurement is absent: the clinical event prevents the scheduled measurement from existing as an outcome under the participant’s observed course.
These cases call for different estimands and analyses. Treating every unavailable value as if it were an unrecorded version of the same underlying measurement can silently change the scientific question. The analysis should first state what outcome is being compared, for whom, and how the event is handled. The earlier material on Post-Baseline Biomarkers, Time Ordering, and Selection Bias is relevant because analyses restricted to people with observed post-baseline values can select different kinds of participants in the treatment groups.
An estimand is the precise treatment-effect quantity the study aims to estimate. At minimum, its specification identifies the target population, the treatment conditions being compared, the outcome variable, how relevant events are handled, and the population-level summary. For a biomarker measured at a fixed visit, possible summaries include an arm-specific mean and the difference between those means. If a clinical event can make the biomarker undefined, the event-handling choice is part of the outcome definition, not a technical detail to be decided after seeing the data.
Estimands for Defined but Unobserved Measurements
Suppose the target is the mean week-12 biomarker under each randomized treatment among the original trial population, and the biomarker remains defined for every participant at week 12. Some values are unrecorded, but the outcome still exists conceptually. An analysis may estimate the full-population mean by using observed information, provided its assumptions about the measurement process are credible.
A common assumption is missing at random conditional on recorded information: within treatment arm, the probability that a value is observed may depend on measured history, but not on the unobserved biomarker itself after conditioning on that history. If \(H\) denotes the recorded history used in the analysis, this assumption can be written as \(R(t)\) independent of \(Y(t)\) conditional on \(A\) and \(H\), among participants for whom \(Y(t)\) is defined. It is not established merely because the data set includes many covariates. It requires that the relevant predictors of both measurement and biomarker value are adequately recorded and modeled. Positivity is also needed: each history pattern targeted by the analysis must have a nonzero chance of a recorded measurement.
Under these assumptions, multiple imputation can draw plausible values from a model for the missing biomarker conditional on treatment and recorded history; the treatment contrast is then estimated across completed data sets, with uncertainty reflecting both within- and between-imputation variation. Inverse-probability weighting instead gives more weight to observed participants whose measured histories predict a lower chance of observation. Either approach depends on the model and assumptions. Neither turns an implausible missingness assumption into a fact.
Worked Example: Standardizing for a Predictor of Measurement
Consider a hypothetical randomized study in which a lower week-12 biomarker concentration is favorable. In each arm, 20 participants are divided equally between low- and high-baseline-risk strata. Assume the outcome is defined for everyone, and that within each arm and stratum the observed participants represent the stratum’s mean under a missing-at-random assumption conditional on stratum.
In the treatment arm, 8 of 10 low-risk participants are measured, with an observed mean of 2 ng/mL; 4 of 10 high-risk participants are measured, with an observed mean of 8 ng/mL. The estimated full-arm mean standardizes these stratum means to the full randomized group:
In the control arm, 8 of 10 low-risk participants are measured, with a mean of 3 ng/mL; 8 of 10 high-risk participants are measured, with a mean of 7 ng/mL. Thus:
The estimated treatment-minus-control contrast is \(5-5=0\) ng/mL. In contrast, the unadjusted observed-case means are \((8(2)+4(8))/12=4\) ng/mL for treatment and \((8(3)+8(7))/16=5\) ng/mL for control, a contrast of \(-1\) ng/mL. That observed-case contrast compares different mixtures of baseline-risk strata. The standardized result is justified only if the within-stratum missing-at-random assumption is credible; it is not proof that the unobserved values equal the stratum means.
When a Clinical Event Makes the Biomarker Undefined
If death or another event makes the scheduled biomarker undefined, imputing a value as though it were merely unrecorded answers a different and usually incoherent question. For example, a post-death biomarker concentration is not an unknown measurement that could be recovered by asking for a blood sample. A model can generate a number, but that number is an extrapolation to a hypothetical outcome, not a value observed or defined under the participant’s actual clinical course.
Several estimand strategies are possible, each with a different interpretation:
- Composite strategy: define an outcome that includes both the clinical event and the biomarker, such as being alive with a biomarker response at the scheduled visit. This retains everyone in the randomized population, but the result concerns the combined outcome, not the biomarker alone.
- While-alive summary: compare the biomarker distribution among participants alive and eligible for measurement in each arm. This describes each arm’s observed survivors. Because treatment may affect who survives, the groups need not represent the same underlying people, so the contrast is not automatically a causal effect in a common population.
- Principal-stratum estimand: target participants who would be event-free under either treatment. This “always-event-free” group is defined by the pair of potential event outcomes under both treatments. Randomization alone does not reveal each participant’s membership in that group, so estimating its treatment effect requires additional assumptions.
- Hypothetical strategy: ask what the biomarker would have been under a specified hypothetical condition in which the event did not occur. This requires a meaningful hypothetical intervention or scenario and a model for outcomes that are not observed after the event. It must not be described as an estimate of the actual post-event biomarker.
A composite can be useful when the clinical event itself is unfavorable and the combined endpoint matches the decision of interest. Its definition must be prespecified: for instance, “alive and at least a 30% biomarker reduction” makes death a failure, whereas “alive with a measured biomarker” would mix survival and measurement availability. A while-alive comparison can be descriptively useful, but it should report survival and measurement patterns alongside the biomarker summary. A principal-stratum result may have a focused causal interpretation, but the stratum is latent and often not identified without strong assumptions.
Worked Example: Defining a Composite When Death Precedes Measurement
In a hypothetical randomized trial, 10 participants receive treatment and 10 receive control. At week 12, a biomarker response is defined as a reduction of at least 30%. Two treatment participants and one control participant died before week 12. Among the survivors, 6 of 8 treatment participants and 7 of 9 control participants met the response threshold. The estimand is the difference in the proportion of all randomized participants who are alive and have a response at week 12; death counts as not achieving this composite outcome.
The treatment-arm composite success proportion is \(6/10=0.600\), and the control-arm proportion is \(7/10=0.700\). Therefore:
The estimated difference is \(-10.0\) percentage points: in this hypothetical sample, the proportion alive with a biomarker response is 10.0 percentage points lower under treatment. Among survivors alone, response proportions are \(6/8=0.750\) and \(7/9=0.778\), a difference of about \(-2.78\) percentage points. That survivor-only contrast answers a different question and conditions on a post-randomization event. Neither calculation imputes biomarker values after death. A randomized comparison supports a treatment-group contrast for the prespecified composite, subject to ordinary concerns about sampling uncertainty and trial conduct.
Sensitivity Analyses and Diagnostics
For defined but missing values, a primary analysis under missing at random should generally be accompanied by sensitivity analysis for departures from that assumption. A transparent approach is a delta adjustment: imputed missing values are shifted by a specified amount from their missing-at-random predictions. The shift should be expressed in biomarker units and chosen to represent clinically plausible departures, not selected to produce a preferred conclusion.
Worked Example: A Delta-Adjusted Missing-Not-at-Random Analysis
Suppose a hypothetical trial has 20 participants per arm, and the week-12 biomarker is defined for everyone. Treatment has 15 observed values with mean 6.0 ng/mL, while control has 16 observed values with mean 5.0 ng/mL. For illustration, assume the missing-at-random model predicts each arm’s missing-value mean to equal its observed mean. Let \(\delta_T\) and \(\delta_C\) be shifts, in ng/mL, applied to the predicted missing-value means. The adjusted arm means are:
Under the missing-at-random reference case, \(\delta_T=\delta_C=0\), so the estimated treatment-minus-control difference is \(6.0-5.0=1.0\) ng/mL. Suppose a sensitivity scenario assigns \(\delta_T=-1.0\) ng/mL and \(\delta_C=0.5\) ng/mL. Then the adjusted means are \(6.0+0.25(-1.0)=5.75\) ng/mL and \(5.0+0.20(0.5)=5.10\) ng/mL, giving a difference of \(0.65\) ng/mL. The result has changed by \(0.35\) ng/mL under these specified assumptions.
The formula for the adjusted contrast is \(1.0+0.25\delta_T-0.20\delta_C\) ng/mL. If \(\delta_C=0\), the point estimate reaches zero at \(\delta_T=-4.0\) ng/mL. This is a tipping-point calculation, not a claim that the missing treatment values actually differ by 4.0 ng/mL. Investigators should judge whether that departure is plausible and assess uncertainty as well as the point estimate.
Diagnostics for ordinary missingness include tabulating observation rates by treatment, visit, and recorded prognostic history; checking whether dropout or missed visits relate to prior biomarker values; and examining overlap in estimated observation probabilities. Very small observation probabilities make weighting unstable and signal that the full-population estimate relies heavily on a few participants. Imputation models should respect the biomarker’s scale and include treatment, relevant baseline and post-baseline history, and predictors of missingness. These checks can reveal model weaknesses, but observed data cannot verify missing-at-random for values that were never observed.
For truncation by a clinical event, report event incidence and timing by treatment, and make clear which estimand is being presented. Do not label an analysis “missing-data imputation” if it assigns biomarker values after death without defining a hypothetical estimand. In a while-alive analysis, differences in survival can change the composition of the groups being compared. In a composite analysis, inspect each component separately as well as the overall result, since a contrast can be driven primarily by event frequency rather than biomarker response.
Choosing an Approach and Avoiding Common Errors
| Situation | Possible target | Key requirement or limitation |
|---|---|---|
| Biomarker exists but is unrecorded | Full-population mean or treatment contrast | Requires a credible missingness model, such as missing at random conditional on recorded history, plus adequate overlap |
| Event makes biomarker undefined | Composite including the event | Interpretation is for the prespecified combined outcome, not biomarker alone |
| Biomarker among observed survivors | Arm-specific while-alive summary | Survivor groups may differ because treatment can affect survival |
| Biomarker effect among always-event-free participants | Principal-stratum contrast | Membership is unobserved; identification needs additional assumptions |
| Biomarker under a specified no-event scenario | Hypothetical estimand | Requires a coherent hypothetical scenario and extrapolation beyond observed outcomes |
A frequent mistake is to say, “The week-12 value is missing, so we will impute it under missing at random,” without first asking whether the value is defined. The correction is to classify each unavailable measurement by its reason. If the participant was alive and eligible for measurement, a missing-data analysis may be appropriate under explicit assumptions. If death made the value undefined, choose and state an event strategy; do not disguise an unobserved hypothetical value as an ordinary missing observation.
A second mistake is to interpret the mean among observed survivors as the treatment effect on the biomarker for the original randomized population. The correction is to name the survivor-conditioned target and acknowledge that the treatment arms may contain different survivor populations. A third mistake is to treat a composite result as if it were a pure biomarker effect. The correction is to state the composite rule precisely and report its components, so the clinical interpretation follows the quantity actually estimated.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- A participant is alive and eligible for a scheduled week-8 blood draw, but the sample is lost. Is the biomarker missing or truncated by an event?
- A participant dies before the scheduled biomarker visit. Why is ordinary missing-at-random imputation not, by itself, a suitable treatment of the post-death value?
- In a trial, 12 of 20 treated participants and 10 of 20 control participants are alive with a prespecified biomarker response. Calculate the composite risk difference, treatment minus control.
- What additional interpretation problem arises when comparing biomarker means only among survivors if treatment affects survival?
- In a delta-adjusted analysis, 30% of an arm’s values are missing. If the missing-value mean is shifted by \(-2\) biomarker units from the missing-at-random prediction, by how much does the full-arm mean shift?
Solutions. (1) It is ordinary missingness: the biomarker remains defined, but \(R=0\). (2) Death makes the scheduled value undefined under the observed course; an imputed number would represent a hypothetical outcome and would require a separately defined hypothetical estimand and assumptions. (3) The difference is \(12/20-10/20=0.600-0.500=0.100\), or 10.0 percentage points in favor of treatment for this composite. (4) The survivor groups may represent different populations because treatment can change who survives; their contrast is not automatically a causal effect for a common population. (5) The full-arm mean shifts by the missing fraction times the shift: \(0.30(-2)=-0.600\) biomarker units.