Tutorials › Model-Informed Drug Development (MIDD) › Representing Nonadherence and Treatment Discontinuation

Clinical Trial Simulation · Tutorial 282 of 596

Representing Nonadherence and Treatment Discontinuation

Build a simulation that links dose-taking behavior and discontinuation to drug exposure and to the treatment strategy being evaluated.

Advanced 12 min read

What You'll Learn

  • Represent missed doses, temporary interruptions, and permanent discontinuation as distinct processes.
  • Specify adherence mechanisms that may vary with participant characteristics and history.
  • Propagate actual dose histories through a pharmacokinetic model.
  • Keep follow-up outcomes after discontinuation separate from missing outcomes.
  • Distinguish a treatment-policy outcome from a hypothetical outcome under sustained exposure.
  • Diagnose implausible adherence patterns and exposure consequences in simulated trials.

Why Dose-Taking Behavior Belongs in a Trial Simulation

A prescribed regimen is not necessarily the regimen a participant receives. A participant may miss an occasional dose, pause treatment for several days, reduce doses informally, or permanently discontinue the assigned drug. These behaviors change the input to the pharmacokinetic and pharmacodynamic models, and therefore may change the outcomes a trial is designed to measure. A simulation that gives every participant every scheduled dose can misrepresent both exposure and the performance of a trial design when departures from the regimen are plausible.

The simulation should represent what happens to dosing separately from whether an outcome is observed. Discontinuation is not itself dropout: participants may stop study drug but continue scheduled assessments, and their later outcomes remain relevant to the treatment strategy they were assigned. The earlier material on dropout and informative missingness addresses missing outcome records; here the focus is on dose histories and their consequences.

Definition: For participant \(i\) and scheduled dose occasion \(k\), let \(A_{ik}=1\) if the dose is taken and \(A_{ik}=0\) if it is not. The administered dose is \(D_{ik}^{\mathrm{admin}}=A_{ik}D_{ik}^{\mathrm{prescribed}}\) while the participant remains on treatment. A missed dose is an isolated or occasional \(A_{ik}=0\); an interruption is a temporary period during which scheduled doses are not taken; and discontinuation is permanent cessation of the assigned study drug. These are different events and need not have the same causes or consequences.

An interruption can be represented by a run of zero dose indicators, with a separate process determining when the run begins and ends. Discontinuation can be represented by a stopping time \(T_i^{\mathrm{disc}}\): after that time, assigned-drug doses are zero. If the protocol permits rescue medication, background therapy, or a restart, those actions require their own rules. A safety hold or restart should follow the relevant protocol rule rather than being silently treated as ordinary nonadherence.

Specify the Behavior Process Before Simulating It

The simplest model makes dose-taking decisions independently with a constant probability \(q\): \(A_{ik}\sim\mathrm{Bernoulli}(q)\). This may be a useful starting scenario, but it assumes that each dose decision is independent of prior decisions and that the probability of taking a dose is the same for every participant and occasion. Those assumptions are often too strong. Adherence may depend on treatment burden, participant characteristics, adverse events, perceived benefit, prior missed doses, or access to medication.

A more realistic specification models adherence conditionally on information available at that time. For example, a probability model may use participant covariates and recent history to determine the chance of taking the next dose. An alternative is a transition model: the probability of a missed dose may be higher after a missed dose, allowing missed doses to cluster. Temporary interruptions can be simulated with start and restart rules, while permanent discontinuation can have a separate event model. The key is temporal order: information generated in the simulation should influence only subsequent behavior, not earlier dose decisions.

The event model should also distinguish causes when they imply different future histories. A safety-related hold may end after a stated recovery condition; a participant who stops because of perceived lack of benefit may be less likely to restart. If the simulation includes these mechanisms, it must state how they operate and how they affect both dosing and follow-up. A model that makes adherence independent of response and safety is a defensible scenario only when that independence is intended, not a default truth.

Worked Example: Probability of a Given Number of Missed Doses

Suppose, hypothetically, that a participant has 10 scheduled doses and takes each dose independently with probability \(q=0.90\). Calculate the probability of exactly two missed doses. This example assumes a constant adherence probability, no interruptions, and independence across occasions.

The probability of taking exactly 8 of the 10 doses is binomial:

$$ P(\text{exactly 2 missed})=\binom{10}{2}(0.90)^8(0.10)^2 $$

Substituting and calculating gives

$$ 45(0.43046721)(0.01)=0.1937102445\approx 0.1937. $$

Thus, under these assumptions, the probability of exactly two missed doses is about 19.37%. This result concerns the count of missed doses; it says nothing about their timing. Two adjacent missed doses and two widely separated missed doses have the same probability in this independent model, even though their exposure consequences may differ.

To represent dependence, one can instead define probabilities conditional on the preceding state. For instance, let the next dose be missed with probability 0.10 after a taken dose and 0.40 after a missed dose. The resulting process allows missed doses to cluster. A simulation should be explicit about whether it uses independent dose decisions, a history-dependent model, an interruption process, or a combination.

Translate the Dose History Into Exposure

The dose-taking process matters only through its consequences for treatment and outcomes. For each simulated participant, pass the administered dose history—not merely the prescribed regimen—into the PK model. Then use the resulting exposure history in the PD or exposure-response model. A missed dose does not necessarily cause an immediate proportional reduction in concentration: the magnitude and timing of the change depend on elimination, absorption, dosing interval, and the model’s state at the time of the omission.

Worked Example: A Two-Day Interruption in a Simple PK Model

Consider a hypothetical one-compartment model sampled once per day after each dose. Let \(C_k\) be the concentration immediately after dose occasion \(k\), and let \(r=0.8187\) be the fraction of the previous concentration remaining one day later. Each administered dose adds \(1.00\ \mathrm{mg/L}\) to concentration. The recurrence is \(C_{k+1}=rC_k+1.00A_{k+1}\), with \(C_0=0\). Assume the participant takes the first two doses, misses the next two during an interruption, and takes the fifth dose.

The administered-dose indicators are \(1,1,0,0,1\). Applying the recurrence in order:

$$ \begin{aligned} C_1&=0.8187(0)+1.00=1.000\ \mathrm{mg/L},\\ C_2&=0.8187(1.000)+1.00=1.819\ \mathrm{mg/L},\\ C_3&=0.8187(1.819)+0=1.489\ \mathrm{mg/L},\\ C_4&=0.8187(1.489)+0=1.219\ \mathrm{mg/L},\\ C_5&=0.8187(1.219)+1.00=1.998\ \mathrm{mg/L}. \end{aligned} $$

For comparison, if all five doses had been taken, the fifth concentration would be \(1+r+r^2+r^3+r^4\) in \(\mathrm{mg/L}\):

$$ 1+0.8187+0.6703+0.5488+0.4493=3.487\ \mathrm{mg/L}. $$

The interruption therefore produces \(1.998\ \mathrm{mg/L}\) rather than \(3.487\ \mathrm{mg/L}\) at this occasion under the stated model. This is a deliberately simple model, not a general prediction: real PK models may include absorption compartments, nonlinear clearance, multiple dosing compartments, and between-participant variability. The calculation illustrates why the placement of missed doses matters.

A short justification follows directly from the model’s state update. At each occasion, the previous concentration contributes the fraction \(rC_k\) that remains, and the new administered dose contributes \(1.00A_{k+1}\). When \(A_{k+1}=0\), the new-dose term is zero but the existing concentration does not disappear. Repeating that update through an interruption gives the concentration trajectory; restarting adds doses to the residual concentration left at restart.

Discontinuation, Follow-Up, and the Outcome Strategy

After permanent discontinuation, assigned-drug dosing stops, but the participant’s clinical course need not stop. Drug concentrations may decline gradually; symptoms or disease status may persist or change; rescue treatment may begin; and scheduled measurements may continue. A trial simulation should therefore specify post-discontinuation PK/PD and care as appropriate to the intended outcome. Setting all outcomes after discontinuation to missing confuses an exposure event with a data-collection event.

Two outcome strategies are particularly important to distinguish. Under a treatment-policy strategy, the outcome is evaluated regardless of missed doses or discontinuation, including the consequences of subsequent care as specified by the trial and model. Under a hypothetical sustained-exposure strategy, the outcome is defined under a hypothetical intervention in which the assigned treatment is taken as specified throughout the relevant period. The second is not the observed outcome among only those who adhered: selecting adherent participants can select a different mix of prognosis, response, and tolerability.

FeatureTreatment-policy outcomeHypothetical sustained-exposure outcome
Doses after a missed dose or stopUse the simulated actual dose historyImpose the stated sustained regimen
Follow-up after discontinuationRetain outcomes and relevant subsequent carePredict the outcome under the hypothetical sustained regimen
InterpretationOutcome under assignment with departures representedOutcome under the specified no-departure exposure scenario
Main cautionDoes not isolate the pharmacologic effect of perfect adherenceDepends on assumptions about outcomes under unobserved sustained exposure

Worked Example: Comparing Two Outcome Strategies

In a hypothetical simulation, four participants are assigned active treatment and four are assigned control. The model generates week-4 benefit scores, with larger values indicating greater benefit. For the active-treatment participants, outcomes under their simulated actual dose histories are 8, 6, 5, and 1 units. Their model-predicted outcomes under sustained exposure are 8, 8, 7, and 7 units. Control outcomes are 2, 2, 3, and 1 units. Assume all outcomes are observed, and the treatment-policy comparison uses the actual-history active outcomes and the control outcomes.

The active-arm treatment-policy mean is

$$ \bar{x}_{\mathrm{active,policy}}=\frac{8+6+5+1}{4}=5.00\ \text{benefit units}. $$

The control mean is

$$ \bar{x}_{\mathrm{control}}=\frac{2+2+3+1}{4}=2.00\ \text{benefit units}. $$

The simulated treatment-policy contrast is therefore \(5.00-2.00=3.00\) benefit units. For the hypothetical sustained-exposure strategy, the active-arm mean is

$$ \bar{x}_{\mathrm{active,sustained}}=\frac{8+8+7+7}{4}=7.50\ \text{benefit units}, $$

and its contrast against the same control mean is \(7.50-2.00=5.50\) benefit units. The 2.50-unit difference between active-arm means reflects the model’s hypothetical change from actual dosing to sustained dosing in this example. It is not, by itself, an estimate of a causal effect among observed adherers. The sustained-exposure result relies on the model’s assumptions about how each outcome would change if dosing were maintained.

Simulate Discontinuation as an Event Process

A discrete-time discontinuation model can assign a probability of stopping in each interval among participants still on treatment. Let \(h_k\) be the conditional probability of discontinuation during interval \(k\), given that the participant has not previously discontinued. If the hazards are fixed for this example, the probability of remaining on treatment through \(K\) intervals is the product of the conditional probabilities of not stopping:

$$ P(T^{\mathrm{disc}}>K)=\prod_{k=1}^{K}(1-h_k). $$

This product follows by applying the definition of conditional probability successively: survival through interval \(k\), conditional on survival through \(k-1\), has probability \(1-h_k\). Multiplying these conditional probabilities gives the probability of surviving all \(K\) intervals. If hazards depend on evolving safety, response, or adherence history, use the corresponding conditional \(h_{ik}\) for each participant and interval rather than a fixed common value.

Worked Example: A Time-Varying Discontinuation Risk

Suppose, hypothetically, that among participants still receiving treatment, the conditional probabilities of discontinuing in four successive weeks are 0.05, 0.10, 0.15, and 0.20. Assume these values apply to every participant and that stopping is permanent. The probability of remaining on treatment through week 4 is

$$ (1-0.05)(1-0.10)(1-0.15)(1-0.20) =(0.95)(0.90)(0.85)(0.80)=0.5814. $$

Thus the probability of discontinuing by the end of week 4 is \(1-0.5814=0.4186\), or 41.86%, under this model. The interval probabilities are conditional on still being on treatment at the start of each interval; they are not four unconditional probabilities that can simply be added. If discontinuation is driven by simulated adverse events or lack of response, the hazards should be updated using those histories, and the resulting event rate need not equal this illustrative calculation.

Diagnostics, Alternatives, and Limitations

Diagnostics should examine the behavior process as well as its downstream implications. Review the proportions of doses taken, missed, and interrupted; the distribution of run lengths for missed doses; time to discontinuation; and these summaries by relevant participant characteristics and simulated histories. Check whether simulated probabilities are in range and whether events occur in the intended order. Then inspect resulting exposure trajectories, including accumulation, decline during interruption, and washout after stopping. Implausible exposure patterns may indicate an implementation error or an inconsistent link between dose records and the PK model.

A constant independent probability is transparent and easy to stress-test, but it cannot represent clustering or behavior linked to response and tolerability. A history-dependent model can represent clustering, but may require more assumptions and data to calibrate. A rule-based protocol model is often appropriate for mandated safety holds and restarts; it should not be used to imply that all participant-driven nonadherence follows a deterministic rule. More detailed event-time models can represent variable timing, but their additional complexity is useful only if the relevant parameters and consequences are supportable.

The principal limitation is that adherence is difficult to predict exactly, particularly when it depends on unmeasured or changing motivations. A single “best guess” process can therefore create false precision. Use plausible alternative adherence and discontinuation scenarios where uncertainty could affect the design decision. Also distinguish the model’s treatment-policy predictions from predictions under sustained exposure: the latter require counterfactual assumptions and may be sensitive to how response, safety, and post-discontinuation care are represented.

  • Pitfall: Treating discontinuation as dropout and deleting every later outcome. Correction: Stop assigned-drug dosing at discontinuation, but continue outcome simulation and observation when the trial collects follow-up data.
  • Pitfall: Multiplying the prescribed dose by an average adherence percentage and using that reduced dose at every occasion. Correction: Simulate dose occasions or interruptions when timing matters; the same average number of doses can produce different exposure trajectories.
  • Pitfall: Calling outcomes among simulated adherers “the sustained-exposure effect.” Correction: Define the hypothetical sustained regimen and simulate outcomes under that intervention; adherers may differ systematically from nonadherers.
Key takeaway: Simulate the dose history, discontinuation, and outcome follow-up as related but distinct processes. Feed actual administered doses into the exposure model for treatment-policy outcomes, and represent sustained exposure as a separate, explicitly hypothetical strategy.

Check Your Understanding

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

  1. In a simulation with a temporary three-day interruption, what distinguishes it from permanent discontinuation?
  2. Under independent dose-taking with probability 0.80, what is the probability of taking all three scheduled doses?
  3. Why should a missed dose not generally be represented by setting concentration to zero?
  4. A participant permanently stops study drug but attends every later assessment. Should those outcomes be missing under a treatment-policy strategy?
  5. Why can the mean outcome among simulated adherers differ from the hypothetical sustained-exposure outcome?

Solutions. (1) The interruption is temporary and requires a modeled restart or end condition; permanent discontinuation sets future assigned-drug doses to zero. (2) By independence, \(0.80^3=0.512\), so the probability is 51.2%. (3) Drug remaining from earlier doses may persist; concentration evolves according to the PK model even when the current dose is omitted. (4) No. If assessments are collected, retain those outcomes and represent the actual post-discontinuation course specified by the model. (5) Adherers may differ in prognosis, response, or tolerability, whereas the sustained-exposure outcome is defined under an intervention imposing the regimen, not by selecting participants who happened to adhere.