Why Intermittent Missingness Needs Special Care
A missing visit is intermittent when an outcome is unobserved at one time but observed again at a later time. For example, a participant may have outcomes recorded at visits 1 and 3, but not at visit 2. This pattern differs from dropout: the later measurement remains available and contains information about the unobserved visit. An imputation for visit 2 should therefore respect the observed outcome at visit 3 rather than treating the sequence as though follow-up ended at visit 1.
Delta adjustment represents a sensitivity assumption about the distribution of outcomes that are missing, beyond the reference model used to predict them. In the intermittent setting, the reference distribution should generally condition on all relevant observed information, including observed later outcomes. The delta then modifies that conditional distribution without replacing or changing those observed values.
The direction of the shift must be stated in context. If larger values indicate worse health, a positive delta means the missing outcomes are assumed to be worse, on average, than the reference model predicts. A delta is a sensitivity parameter, not a quantity identified by the observed data alone. Its value should be justified by substantive knowledge or explored across plausible values.
Condition on the Observed Future
For a participant with \(Y_1\) observed, \(Y_2\) missing, and \(Y_3\) observed, a bridge-style reference model uses \(f_0(Y_2\mid Y_1,Y_3,X,R)\). The future value \(Y_3\) is part of the participant’s observed data. It can inform the likely value at the missing visit through the longitudinal association between visits. A delta adjustment changes the assumed distribution of \(Y_2\) conditional on those observations; it does not impute a new \(Y_3\).
The future observation does not necessarily impose a deterministic rule such as \(Y_2\) having to lie between \(Y_1\) and \(Y_3\). Such a rule would require a substantive constraint, such as a known monotone trajectory. In ordinary longitudinal data, outcomes can rise and fall. The later value instead constrains the imputation probabilistically through the conditional model: values of \(Y_2\) that fit poorly with the observed endpoints receive less probability under the reference distribution.
Worked Example: Shifting a Missing Visit Given Both Endpoints
Suppose a hypothetical study records a continuous log-biomarker, in log units, at three visits. A participant has \(Y_1=1.2\), \(Y_2\) missing, and \(Y_3=2.0\). Assume a fitted reference model gives \(Y_2\mid Y_1,Y_3,X,R\sim N(m_0,0.25^2)\), where \(m_0=0.50+0.40Y_1+0.50Y_3\). For this example, suppose the reference model is adequate for the conditional distribution and the participant’s covariates and missingness pattern are represented appropriately.
Reference conditional mean. \(m_0=0.50+0.40(1.2)+0.50(2.0)=0.50+0.48+1.00=1.98\) log units. The observed later value contributes to this prediction through the \(0.50Y_3\) term.
Suppose the sensitivity assumption is \(\delta=0.30\) log units, meaning missing visit-2 values are shifted upward relative to the reference distribution, conditional on the observed data. Under an additive shift,
Thus, the adjusted conditional mean is \(1.98+0.30=2.28\) log units, and the conditional standard deviation remains \(0.25\) log units under this location-shift specification. The observed \(Y_3=2.0\) remains exactly 2.0; it is not shifted. The calculation assumes that a constant additive shift is meaningful on the log scale and that the reference conditional model adequately captures how the endpoints relate to the missing visit.
Why the Conditional Construction Is Coherent
A useful justification comes from factorizing the joint distribution into an observed-data distribution and a conditional distribution for the missing components. For simplicity, condition throughout on \(X\) and \(R\). Let \(f_0(Y_O)\) be the reference distribution of the observed outcomes and \(f_0(Y_M\mid Y_O)\) the reference conditional distribution for the missing outcomes. Define the sensitivity model by keeping the observed-data distribution fixed and replacing the conditional distribution:
This is a proper joint distribution whenever \(f_\delta(Y_M\mid Y_O)\) is a normalized conditional distribution for each \(Y_O\). Indeed, integrating over the missing values gives \(\int f_\delta(Y_O,Y_M)\,dY_M=f_0(Y_O)\). Integrating over the observed values then gives 1. Consequently, the construction changes the assumed distribution of missing values while retaining the reference distribution for the observed values.
For an unbounded continuous outcome, an additive shift can be written as \(f_\delta(y_M\mid y_O)=f_0(y_M-\delta\mid y_O)\). A change of variable \(u=y_M-\delta\), with \(dy_M=du\), shows that the shifted density integrates to 1. For several missing visits, the corresponding vector shift is \(f_{\boldsymbol{\delta}}(\mathbf y_M\mid\mathbf y_O)=f_0(\mathbf y_M-\boldsymbol{\delta}\mid\mathbf y_O)\), provided the shifted distribution is well-defined on the outcome’s support. The joint conditional distribution matters: shifting components separately while ignoring their dependence can produce a different and potentially incoherent longitudinal model.
Worked Example: Why the Later Visit Changes the Imputation
Return to a hypothetical participant with \(Y_1=1.2\), missing \(Y_2\), and observed \(Y_3=2.0\). Suppose the reference model conditioned only on past history would give mean \(1.50+0.60Y_1\), whereas a bridge model using both observed endpoints gives mean \(0.50+0.40Y_1+0.50Y_3\). These are stipulated fitted models for illustration; in an application, each would need adequate data and diagnostic checking.
History-only prediction: \(1.50+0.60(1.2)=1.50+0.72=2.22\) log units.
Bridge prediction: \(0.50+0.40(1.2)+0.50(2.0)=1.98\) log units.
With the same \(\delta=0.30\) log-unit sensitivity shift, the adjusted conditional means are \(2.52\) and \(2.28\) log units, respectively. The 0.24-unit difference arises from conditioning on the observed later outcome in the bridge model; it is not a difference in the delta. Neither calculation is automatically preferable in every setting, but when \(Y_3\) is observed and relevant, omitting it discards information that can help constrain \(Y_2\).
Several Intermittently Missing Visits
When multiple visits are missing but later and earlier visits are observed, specify the missing outcomes jointly conditional on the observed outcomes. This makes the dependence among missing visits explicit. For example, if visits 2 and 4 are missing while visits 1, 3, and 5 are observed, use a reference distribution for \((Y_2,Y_4)\) conditional on \((Y_1,Y_3,Y_5)\), then define a sensitivity shift for each missing visit. The observed values at visits 1, 3, and 5 remain fixed.
A vector of deltas allows the sensitivity assumptions to differ by visit. They may also depend on covariates or the missingness pattern if that is scientifically justified. Such flexibility does not remove the need for a coherent joint conditional model: the adjusted distributions should retain reasonable dependence across missing visits and should be compatible with any constraints on the outcome.
Worked Example: Two Missing Visits With Dependence Preserved
Suppose a hypothetical longitudinal log-biomarker is observed at visits 1, 3, and 5, while visits 2 and 4 are missing. Given the three observed values and covariates, suppose a reference joint conditional model is \((Y_2,Y_4)^\mathsf{T}\sim N(\mathbf m,\Sigma)\), with \(\mathbf m=(1.70,2.10)^\mathsf{T}\) and \(\Sigma=\begin{pmatrix}0.040&0.012\\0.012&0.090\end{pmatrix}\). Assume this conditional model is adequate and that an additive shift on the log scale is scientifically interpretable.
Choose visit-specific sensitivity shifts \(\boldsymbol{\delta}=(0.20,-0.10)^\mathsf{T}\) log units. The adjusted mean vector is \(\mathbf m+\boldsymbol{\delta}=(1.70+0.20,2.10-0.10)^\mathsf{T}=(1.90,2.00)^\mathsf{T}\). Under a location shift, the conditional covariance remains \(\Sigma\); in particular, the covariance between visits 2 and 4 remains \(0.012\), rather than being discarded.
For illustration, if a draw from the reference conditional distribution were \((1.50,2.40)^\mathsf{T}\), the corresponding shifted draw would be \((1.50+0.20,2.40-0.10)^\mathsf{T}=(1.70,2.30)^\mathsf{T}\). The observed visit values used to form \(\mathbf m\) are unchanged. This example assumes the shifted joint distribution remains appropriate on the log scale; if a shift conflicts with known support or clinical constraints, the model must be reformulated rather than used mechanically.
Assumptions, Diagnostics, and Alternatives
Delta adjustment does not make an inadequate reference model adequate. The reference conditional distribution must represent the relevant observed history, future outcomes, covariates, and missingness pattern, with sufficient overlap to estimate the needed conditional relationships. A common reference is based on a missing-at-random assumption: conditional on observed data included in the model, the missingness does not further depend on the unobserved outcome. The delta then expresses a chosen departure from that reference assumption. The observed data generally cannot determine whether the delta is correct.
Check the reference model before interpreting adjusted imputations. Useful diagnostics include comparing predicted and observed outcome distributions at visits with data, checking conditional residual patterns, and examining whether simulated trajectories reproduce the dependence among observed visits. After applying delta, inspect adjusted distributions and complete trajectories for implausible values, abrupt changes, or contradictions with known clinical constraints. If the outcome has bounded support, a simple additive shift of an unbounded model may assign probability outside the valid range; use a distribution and adjustment that respect the support. A numerical delta that pushes many imputations into impossible values is a model warning, not a reason to ignore those values.
A history-only imputation model predicts a missing visit from information available before that visit. It can be appropriate when the future outcome is unavailable at the time of prediction or when the analysis specifically requires that conditioning set. For retrospective imputation of an intermittently missing outcome, however, an observed later outcome is usually relevant information and should not be omitted without reason. Last observation carried forward is not a substitute for a conditional model: it fixes the missing value to a prior observed value and does not represent uncertainty or the information in a later observed visit.
A selection model is another standard framework for sensitivity analysis, but it parameterizes the relationship between outcome values and their observation probabilities rather than directly specifying the missing-outcome distribution conditional on observed outcomes. The choice of framework does not make the sensitivity assumption empirically identified. Delta adjustment is useful when it is easier to state a clinically meaningful departure in outcome units, but its results depend on both the reference model and the chosen delta values.
Common Pitfall
A second error is to shift an entire participant trajectory, including observed visits. Delta adjustment here concerns the distribution of missing outcomes. Changing observed values would alter the data rather than express a sensitivity assumption about unobserved outcomes.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- A participant has \(Y_1=4\), \(Y_2\) missing, and \(Y_3=8\). A reference conditional mean is \(0.5+0.25Y_1+0.50Y_3\). With an additive delta of \(-0.4\) outcome units, what is the adjusted conditional mean?
- Why is it generally inappropriate to shift an observed later outcome when applying a delta to a missing visit?
- For two missing visits, why might a joint conditional distribution be preferable to two unrelated univariate imputations?
- Name one diagnostic concern that arises when an additive shift is used for a bounded outcome.
Solutions. (1) The reference mean is \(0.5+0.25(4)+0.50(8)=0.5+1+4=5.5\). The adjusted mean is \(5.5-0.4=5.1\) outcome units. (2) The delta represents an assumption about missing outcomes; the observed later value is part of the data and remains fixed. It can inform the conditional distribution of the missing visit without itself being changed. (3) A joint distribution represents dependence between the missing visits. Separate univariate imputations may discard that dependence and generate implausible combinations of values. (4) A shifted distribution may assign probability to values outside the permitted range. The analyst should inspect support violations and use a model and adjustment that respect the outcome scale.