Tutorials › Biostatistics › Time-to-Event CRM (TITE-CRM)

Phase I & Adaptive Study Design

Time-to-Event CRM (TITE-CRM): A Worked Example

A computational guide to the Time-to-Event Continual Reassessment Method (TITE-CRM), including partial follow-up weighting, Bayesian dose-toxicity modeling, posterior updating, dose assignment, safety considerations, and a complete worked example.

Advanced 18 min read

What You'll Learn

  • Why conventional CRM can become inefficient when DLTs occur late
  • How TITE-CRM uses partial follow-up information from patients still at risk
  • How the Bayesian dose-toxicity model is updated as patients accumulate
  • How partial follow-up weights enter the likelihood
  • How a TITE-CRM dose recommendation is obtained from the posterior model
  • How to implement the core calculations in R

Introduction

Phase I dose-finding trials are designed to learn about the relationship between dose and toxicity while identifying a dose whose probability of dose-limiting toxicity (DLT) is close to a prespecified target.

The Continual Reassessment Method (CRM) is a model-based Bayesian approach that updates a dose-toxicity model as patient outcomes accumulate. Rather than relying on a fixed rule such as the traditional 3+3 algorithm, CRM uses information from the entire accumulated trial to determine the dose for the next patient or cohort.

The original CRM, however, generally requires toxicity outcomes to be observed before they can fully contribute to the dose assignment process. This can become problematic when DLTs occur relatively late in the predefined observation window.

For example, suppose every patient must be followed for 28 days before the DLT outcome is considered complete. If patients are enrolled one at a time and the trial waits 28 days after every patient, the dose-escalation process can become extremely slow.

Key idea: TITE-CRM extends CRM by allowing patients to enter the trial in a staggered fashion while incorporating the information available from patients who have not yet completed their full DLT observation period.

Why Time-to-Event Information Matters

Consider a patient who has received treatment at a particular dose and has been followed for only 7 days of a planned 28-day DLT window.

The patient has not experienced a DLT.

A conventional CRM analysis might treat this patient as having an unresolved outcome and therefore delay the next dose decision.

TITE-CRM instead recognizes that the patient has already contributed some information.

If the patient remains DLT-free through 7 days, that provides evidence against the hypothesis that the dose is highly toxic, although it is obviously less information than 28 complete days of follow-up.

A simple TITE-CRM implementation therefore assigns this patient a follow-up weight of:

$$ w_i=\frac{7}{28}=0.25 $$

The patient contributes approximately one quarter of the information that would be available after completing the entire observation window.

Important: The weight is not a statement that the patient has a 25% chance of having a DLT. It is an information weight representing the amount of follow-up currently available under the TITE-CRM approximation.

CRM Before TITE-CRM

To understand TITE-CRM, first consider the standard CRM.

Suppose there are \(J\) dose levels:

$$ d_1

For each dose level, the design specifies an initial estimate of the probability of DLT.

These values form the CRM skeleton:

$$ (p_1,p_2,\ldots,p_J) $$

The skeleton should be monotonically increasing:

$$ p_1

because the working assumption is that higher doses are associated with greater toxicity probability.

The skeleton is not necessarily treated as the true toxicity curve. Instead, it provides the initial shape of the dose-toxicity relationship that is subsequently updated using observed data.

The Target Toxicity Level

The investigator specifies a target probability of DLT, denoted by:

$$ \phi $$

For example, suppose the target toxicity probability is 25%:

$$ \phi=0.25 $$

The goal is to identify the dose whose estimated DLT probability is closest to 25%.

This dose is often referred to as the estimated maximum tolerated dose, or MTD, although the precise terminology and decision rule should be defined in the protocol.

Design principle: The target toxicity level is a clinical decision, not a number that should be selected simply because it produces a convenient statistical design.

The CRM Dose-Toxicity Model

A common one-parameter CRM model uses the power model:

$$ p_j(\theta) = p_j^{\exp(\theta)} $$

where:

  • \(p_j\) is the skeleton probability at dose \(j\)
  • \(\theta\) is an unknown model parameter
  • \(p_j(\theta)\) is the modeled probability of DLT at dose \(j\)

The exponential transformation ensures that the model remains monotone in the dose index when the skeleton is monotone.

When:

$$ \theta=0 $$

the model reproduces the skeleton:

$$ p_j(0)=p_j $$

If \(\theta>0\), the probabilities are shifted downward because \(\exp(\theta)>1\) and \(0

If \(\theta<0\), the probabilities are shifted upward.

Do not interpret \(\theta\) as a direct toxicity probability. It is a model parameter controlling the location of the dose-toxicity curve. The clinically relevant quantities are the resulting \(p_j(\theta)\) values.

The Bayesian Prior

CRM uses a prior distribution for the unknown parameter.

For this tutorial, suppose:

$$ \theta\sim N(0,1.34^2) $$

The posterior distribution is then obtained by combining this prior with the likelihood contributed by the observed patient data.

In practice, prior specification and skeleton calibration should be performed carefully because the prior can influence early dose assignments, particularly when little information has accumulated.

What Changes in TITE-CRM?

The fundamental CRM model remains the same.

The key modification is the treatment of incomplete follow-up.

For patient \(i\), define:

  • \(Y_i=1\) if a DLT has occurred
  • \(Y_i=0\) if no DLT has yet occurred
  • \(u_i\) = amount of follow-up currently observed
  • \(\tau\) = planned DLT observation window

For the basic TITE-CRM weighting scheme:

$$ w_i= \begin{cases} 1, & Y_i=1\\[6pt] \frac{u_i}{\tau}, & Y_i=0 \end{cases} $$

The no-DLT weight is typically bounded at 1:

$$ w_i=\min\left(\frac{u_i}{\tau},1\right) $$

A patient who experiences a DLT receives weight 1 because the event has already been observed.

A patient with no DLT receives a weight reflecting the fraction of the planned observation period that has elapsed.

Examples of TITE Weights

Follow-Up Observation Window Weight if DLT-Free
0 days 28 days 0.000
7 days 28 days 0.250
14 days 28 days 0.500
21 days 28 days 0.750
28 days 28 days 1.000
35 days 28 days 1.000

The weight is therefore a measure of the amount of available follow-up under this simple weighting approach.

How the Weighted Likelihood Works

Suppose patient \(i\) is treated at dose \(j\), and let the modeled DLT probability at that dose be:

$$ p_j(\theta) $$

For a conventional binary outcome, the likelihood contribution would be:

$$ p_j(\theta)^{Y_i} \left[1-p_j(\theta)\right]^{1-Y_i} $$

TITE-CRM raises this contribution to the patient's follow-up weight:

$$ L_i(\theta) = \left[ p_j(\theta)^{Y_i} \left\{1-p_j(\theta)\right\}^{1-Y_i} \right]^{w_i} $$

Equivalently:

$$ L_i(\theta) = p_j(\theta)^{w_iY_i} \left[1-p_j(\theta)\right]^{w_i(1-Y_i)} $$

The total weighted likelihood is:

$$ L(\theta) = \prod_{i=1}^{n} p_{j_i}(\theta)^{w_iY_i} \left[ 1-p_{j_i}(\theta) \right]^{w_i(1-Y_i)} $$

The posterior is proportional to the prior multiplied by this likelihood:

$$ \pi(\theta\mid D) \propto L(\theta)\pi(\theta) $$
This is the central mathematical idea behind TITE-CRM: patients with incomplete DLT follow-up are allowed to contribute information without being treated as if their full observation period were complete.

Why a DLT Gets Weight 1

Suppose a patient develops a DLT on day 10 of a 28-day observation window.

The DLT is already known.

There is no need to discount the event simply because the patient did not complete the remaining 18 days.

The relevant outcome for the dose-toxicity model is already observed:

$$ Y_i=1 \qquad\Rightarrow\qquad w_i=1 $$

The remaining follow-up may matter for other safety endpoints, but it does not make the already-observed DLT more or less of a DLT.

Why a DLT-Free Patient Gets a Partial Weight

Now suppose a patient has completed only 7 days without a DLT.

There are still 21 days during which the patient could experience a DLT.

Under the simple linear weighting approach:

$$ w_i=\frac{7}{28}=0.25 $$

The patient therefore contributes a fractional likelihood contribution.

This allows the trial to use the accumulating information without pretending that the patient's outcome is already fully known.

Staggered Patient Entry

The major operational advantage of TITE-CRM is that patients do not necessarily have to wait for previous patients to complete the entire DLT observation window.

For example:

1
Patient 1 begins treatment at the current dose.
2
Patient 2 can enter before Patient 1 completes the 28-day observation period.
3
Patient 3 can enter while Patients 1 and 2 are still being followed.
4
The TITE-CRM model uses each patient's available follow-up when reassessing the dose.

This can substantially reduce trial duration when toxicity outcomes are late relative to the desired enrollment rate.

A Complete Worked Example

Consider an oncology Phase I trial evaluating five dose levels.

The primary endpoint is the occurrence of a dose-limiting toxicity during a 28-day observation period.

The target toxicity probability is:

$$ \phi=0.25 $$

The investigators specify the following skeleton:

Dose Level Skeleton DLT Probability
Dose 1 0.05
Dose 2 0.10
Dose 3 0.20
Dose 4 0.30
Dose 5 0.45

The working model is:

$$ p_j(\theta)=p_j^{\exp(\theta)} $$

and the prior is:

$$ \theta\sim N(0,1.34^2) $$

Trial Configuration

Parameter Value
Number of dose levels 5
Target DLT probability 25%
DLT observation window 28 days
CRM model Power model
Prior \(\theta\sim N(0,1.34^2)\)
Planned maximum sample size 12 patients
Illustrative cohort size 1 patient
This is an illustrative statistical example, not a clinical trial recommendation. Actual TITE-CRM studies require simulation-based design evaluation, clinical review, safety constraints, dose-escalation restrictions, and prespecified stopping rules.

Patient-Level Data at an Interim Look

Suppose eight patients have entered the study.

Patient Dose DLT? Follow-Up TITE Weight
1 1 No 28 days 1.00
2 2 No 28 days 1.00
3 2 No 21 days 0.75
4 3 No 14 days 0.50
5 3 No 28 days 1.00
6 4 Yes 10 days 1.00
7 4 No 7 days 0.25
8 3 No 28 days 1.00

The important feature is that several patients are still being followed.

Patients 3, 4, and 7 have not completed their full 28-day DLT windows.

Nevertheless, they contribute information through their partial follow-up weights.

Step 1: Calculate the Follow-Up Weights

For Patient 3:

$$ w_3=\frac{21}{28}=0.75 $$

For Patient 4:

$$ w_4=\frac{14}{28}=0.50 $$

For Patient 7:

$$ w_7=\frac{7}{28}=0.25 $$

Patient 6 experienced a DLT, so:

$$ w_6=1 $$

The complete vector of weights is therefore:

$$ (1,\;1,\;0.75,\;0.50,\;1,\;1,\;0.25,\;1) $$

Step 2: Write the Dose-toxicity Model

The five skeleton probabilities are:

$$ (0.05,0.10,0.20,0.30,0.45) $$

For any value of \(\theta\), the modeled probabilities are:

$$ \begin{aligned} p_1(\theta)&=0.05^{\exp(\theta)}\\ p_2(\theta)&=0.10^{\exp(\theta)}\\ p_3(\theta)&=0.20^{\exp(\theta)}\\ p_4(\theta)&=0.30^{\exp(\theta)}\\ p_5(\theta)&=0.45^{\exp(\theta)} \end{aligned} $$

Step 3: Construct the Weighted Likelihood

Patient 1 was treated at Dose 1 and had no DLT with complete follow-up. His contribution is:

$$ 1-p_1(\theta) $$

Patient 3 was treated at Dose 2 and had no DLT after 21 days. His contribution is:

$$ [1-p_2(\theta)]^{0.75} $$

Patient 4 was treated at Dose 3 and had no DLT after 14 days. His contribution is:

$$ [1-p_3(\theta)]^{0.50} $$

Patient 6 experienced a DLT at Dose 4. His contribution is:

$$ p_4(\theta) $$

Patient 7 was treated at Dose 4 and had no DLT after 7 days. His contribution is:

$$ [1-p_4(\theta)]^{0.25} $$

Multiplying all eight contributions gives the complete weighted likelihood.

$$ \begin{aligned} L(\theta) ={}& [1-p_1(\theta)] [1-p_2(\theta)] [1-p_2(\theta)]^{0.75}\\ &\times [1-p_3(\theta)]^{0.50} [1-p_3(\theta)] p_4(\theta)\\ &\times [1-p_4(\theta)]^{0.25} [1-p_3(\theta)] \end{aligned} $$

The posterior distribution is proportional to:

$$ \pi(\theta\mid D) \propto L(\theta) \exp\left( -\frac{\theta^2}{2(1.34)^2} \right) $$

Step 4: Estimate the Posterior Parameter

The posterior distribution does not generally have a convenient closed-form expression under this model.

Numerical integration can therefore be used to calculate the posterior mean:

$$ \hat{\theta} = E(\theta\mid D) $$

For the illustrative dataset above, numerical integration gives approximately:

$$ \boxed{\hat{\theta}\approx0.106} $$

This is the posterior mean used in the subsequent dose-toxicity calculation.

Step 5: Convert the Posterior Parameter into DLT Probabilities

Using:

$$ \hat{\theta}=0.106 $$

the fitted DLT probabilities are approximately:

Dose Skeleton Posterior Estimated DLT Probability
1 0.05 3.6%
2 0.10 7.7%
3 0.20 16.7%
4 0.30 26.2%
5 0.45 41.2%

The target toxicity probability is:

$$ \phi=0.25 $$

Dose 4 has an estimated DLT probability of approximately 26.2%, which is very close to the target.

Step 6: Select the Dose Closest to the Target

The standard CRM decision rule is:

$$ d_{\text{next}} = \arg\min_j \left| p_j(\hat{\theta})-\phi \right| $$

For the example:

Dose Estimated DLT Probability Distance From 25%
1 3.6% 21.4 percentage points
2 7.7% 17.3 percentage points
3 16.7% 8.3 percentage points
4 26.2% 1.2 percentage points
5 41.2% 16.2 percentage points

Therefore:

$$ \boxed{d_{\text{next}}=4} $$

Under this illustrative model and dataset, Dose 4 is the dose whose posterior estimated DLT probability is closest to the 25% target.

Notice what TITE-CRM accomplished: Patient 7 had only seven days of follow-up, yet that patient still contributed information to the model. The trial did not have to wait 21 additional days before performing the statistical update.

What If Patient 7 Had Completed 28 Days?

This provides a useful way to understand the weighting mechanism.

At seven days:

$$ w_7=0.25 $$

After completing the full observation window:

$$ w_7=1 $$

The patient's contribution would therefore become:

$$ 1-p_4(\theta) $$

rather than:

$$ [1-p_4(\theta)]^{0.25} $$

The statistical information contributed by the patient increases as follow-up accumulates.

What If Patient 7 Develops a DLT?

Suppose Patient 7 develops a DLT on day 12.

The patient is no longer treated as a pending DLT-free observation.

Instead:

$$ Y_7=1 \qquad\text{and}\qquad w_7=1 $$

The contribution becomes:

$$ p_4(\theta) $$

This is an important property of time-to-event designs: the exact timing of the event determines when the event becomes known, while a DLT-free patient accumulates information progressively over time.

Conventional CRM vs. TITE-CRM

Feature Conventional CRM TITE-CRM
Model-based dose finding Yes Yes
Bayesian updating Commonly used Commonly used
Uses a dose-toxicity model Yes Yes
Requires complete DLT follow-up before updating Typically No
Uses partial follow-up No Yes
Allows staggered enrollment Limited by outcome availability Yes
Designed for late-onset toxicity Less efficient Yes

The Core Difference in One Equation

The conventional CRM uses a complete binary outcome:

$$ Y_i\in\{0,1\} $$

TITE-CRM supplements this with:

$$ w_i\in[0,1] $$

and therefore uses:

$$ Y_iw_i $$

within the weighted likelihood.

The conceptual transition is:

A
CRM: wait until the binary toxicity outcome is known.
B
TITE-CRM: use the toxicity outcome if known and partial follow-up if it is not yet known.

Why TITE-CRM Can Shorten a Trial

Suppose the DLT window is 28 days.

A traditional sequential design that waits for every patient to complete the window could require approximately 28 days between dose decisions.

TITE-CRM can allow patients to enter while previous patients are still being followed.

The actual trial duration then depends on:

  • Patient accrual rate
  • DLT observation-window length
  • Timing of observed DLTs
  • Frequency of interim model updates
  • Operational safety rules
  • Maximum sample size

Therefore, TITE-CRM is particularly attractive when the toxicity observation window is long relative to the desired accrual interval.

But TITE-CRM Does Not Eliminate Uncertainty

A patient with seven days of DLT-free follow-up is not equivalent to a patient with 28 days of DLT-free follow-up.

The weighting scheme attempts to reflect this difference.

Consequently, the model is making decisions with incomplete information.

This is a trade-off: TITE-CRM gains enrollment efficiency by using partial follow-up, but the statistical decision process must therefore be carefully calibrated through simulation under plausible toxicity and event-time scenarios.

Time-to-Event Information Is Not Necessarily Linear

The simple weight:

$$ w_i=\frac{u_i}{\tau} $$

assumes that information accumulates linearly with follow-up.

That assumption is convenient but need not be clinically realistic.

If DLTs are much more likely early in the observation window, seven days of follow-up may provide more information than one quarter of the total information.

Conversely, if DLT risk is concentrated near the end of the window, linear weighting may overstate the information obtained early.

This motivates alternative time-to-event weighting approaches.

Generalized Follow-Up Weights

Instead of a linear weight, one can define:

$$ w_i=g(u_i) $$

where \(g(\cdot)\) reflects an assumed pattern of information accumulation.

A simple linear approach is:

$$ g(u)=\frac{u}{\tau} $$

More sophisticated approaches can incorporate assumptions about the hazard or timing distribution of DLTs.

Design implication: The weighting function is part of the statistical design. It should therefore be specified and justified before the trial rather than selected after seeing the observed toxicity timing.

Safety Restrictions Are Essential

The original CRM framework attracted criticism because an unconstrained model could recommend aggressive dose escalation when little toxicity information was available.

Modern implementations commonly incorporate safety modifications.

Examples include:

  • No dose skipping during escalation
  • Maximum permitted escalation per decision
  • Minimum number of patients before escalation
  • Requirements for sufficient follow-up at a dose
  • Posterior overdose-control rules
  • Trial stopping rules
  • Maximum number of patients treated at a dose

These modifications are not interchangeable with the underlying TITE-CRM model. They are additional design components intended to control the practical risks of adaptive dose assignment.

No Dose Skipping

A common safety restriction is that the trial may not skip an untested dose during escalation.

For example, if the current dose is Dose 2, the next patient may be permitted to receive Dose 3 but not Dose 4.

1
Current dose = Dose 2.
2
CRM posterior recommends Dose 4.
3
The safety restriction prevents direct escalation to Dose 4.
4
The next patient is instead assigned Dose 3.

Whether dose skipping is allowed is a design choice and must be explicitly specified.

Overdose Control

Another approach is to require that a dose have sufficiently low posterior probability of being excessively toxic before it can be selected.

Suppose the target is:

$$ \phi=0.25 $$

and investigators define an overdose threshold:

$$ \phi_{\text{over}}=0.35 $$

A possible safety criterion could require:

$$ P\left( p_j(\theta)>\phi_{\text{over}} \mid D \right) < c $$

where \(c\) is a prespecified probability cutoff.

This is conceptually different from simply choosing the dose with toxicity probability closest to the target.

Important: The numerical overdose-control threshold and posterior probability cutoff must be prespecified and calibrated through simulation. They should not be invented after observing interim data.

Handling Pending Patients

One of the most important practical questions in TITE-CRM is how pending patients are represented.

Suppose the trial currently has:

Patient DLT Status Follow-Up Weight
A DLT 10 days 1.00
B No DLT 7 days 0.25
C No DLT 14 days 0.50
D No DLT 28 days 1.00

The model incorporates all four patients simultaneously.

Patient A provides a complete DLT event.

Patients B and C provide partial information.

Patient D provides a complete DLT-free observation.

This is what makes TITE-CRM different from simply analyzing completed patients.

Repeated Model Updating

TITE-CRM is called a continual reassessment method because the model is re-estimated as new information becomes available.

For example:

1
Enroll patient.
2
Observe DLT status and current follow-up for all enrolled patients.
3
Calculate TITE weights.
4
Update the posterior distribution of the dose-toxicity parameter.
5
Calculate posterior DLT probabilities for every dose.
6
Apply prespecified safety restrictions.
7
Assign the next patient or cohort.

Stopping Rules

A TITE-CRM protocol should include explicit stopping rules.

Examples include:

  • Maximum total sample size reached
  • Maximum number of patients treated at a dose reached
  • Excessive posterior probability that all doses exceed the target
  • Insufficient evidence that any dose is acceptably safe
  • Predefined safety signal
  • Operational stopping criteria

The exact stopping criteria are design-specific.

Do not confuse dose selection with trial stopping. Selecting the dose closest to the target is one decision. Determining whether the trial should continue at all is another.

Final MTD Selection

At the end of the trial, the estimated MTD is generally selected using the same prespecified model and decision rule used during dose escalation.

For example:

$$ \widehat{MTD} = \arg\min_j \left| p_j(\hat{\theta})-\phi \right| $$

where \(\hat{\theta}\) is the final posterior summary specified by the protocol, such as the posterior mean.

The final selection rule should be prespecified rather than chosen after examining the final toxicity curve.

TITE-CRM Is a Bayesian Adaptive Design

There are three important components:

1
Prior: initial uncertainty about the dose-toxicity relationship.
2
Likelihood: observed DLTs and partial follow-up from enrolled patients.
3
Posterior: updated dose-toxicity relationship used for the next decision.

The Bayesian updating mechanism is therefore:

$$ \text{Prior} \times \text{Weighted Likelihood} \propto \text{Posterior} $$

Interpreting the Posterior Dose-Toxicity Curve

The skeleton provides the initial curve.

The patient data shift the curve according to the evidence accumulated during the trial.

A DLT generally pushes the fitted toxicity curve upward.

DLT-free follow-up generally pushes the curve downward, with the magnitude of the contribution depending on the available follow-up.

This produces a continually evolving estimate of the toxicity probability at each dose.

What Happens After a DLT?

Suppose a DLT occurs at Dose 4.

The DLT contributes a complete event:

$$ L_i(\theta)=p_4(\theta) $$

The posterior toxicity curve will generally shift upward.

The next dose recommendation may therefore remain at Dose 4, move downward, or—depending on the complete accumulated data and safety restrictions—remain near the same dose.

The important point is that the decision uses the entire accumulated dataset, not simply the most recent patient.

What Happens After a DLT-Free Patient?

Suppose a patient at Dose 4 remains DLT-free for 14 days.

The weight is:

$$ w_i=0.50 $$

The contribution is:

$$ [1-p_4(\theta)]^{0.50} $$

This provides evidence that the toxicity probability at Dose 4 may not be as high as initially assumed.

However, because the patient is still within the DLT window, the evidence is only partial.

Why the Model Borrows Information Across Doses

A major difference between CRM and simple rule-based designs is that the model connects the dose levels.

Suppose there are no observed DLTs at Dose 3.

That information does not affect only Dose 3.

Because the dose-toxicity model is shared across doses, the posterior estimate of \(\theta\) changes, which changes the estimated toxicity probabilities at all dose levels:

$$ p_1(\theta),p_2(\theta),\ldots,p_J(\theta) $$

This is one reason CRM can be more statistically efficient than independent dose-by-dose rules.

TITE-CRM vs. 3+3

Feature 3+3 TITE-CRM
Model-based No Yes
Uses all accumulated dose-toxicity data No Yes
Bayesian model No Commonly yes
Uses partial follow-up Generally no Yes
Staggered enrollment Usually constrained Yes
Explicit target toxicity probability Usually implicit Yes
Flexible safety modifications Limited Yes

The important distinction is not that one design is universally superior. Rather, the designs solve the dose-finding problem using fundamentally different statistical mechanisms.

TITE-CRM vs. TITE-BOIN

TITE-CRM and TITE-BOIN both address the same operational problem: DLT outcomes may not be immediately available.

Their underlying statistical approaches differ.

Feature TITE-CRM TITE-BOIN
Underlying framework Model-based CRM Interval-based BOIN
Continuous dose-toxicity model Yes No
Skeleton Yes No in the same CRM sense
Bayesian posterior parameter Yes Design-specific Bayesian calculations
Partial follow-up Yes Yes
Dose assignment Based on fitted toxicity curve Based on interval decision rules

Both approaches can be useful, but they should not be treated as equivalent algorithms.

R Implementation: Define the Design

The following R code establishes the example design.

skeleton <- c(
  0.05,
  0.10,
  0.20,
  0.30,
  0.45
)

target <- 0.25

tau <- 28

sigma <- 1.34

max_n <- 12

Create the Patient-Level Dataset

data <- data.frame(
  patient = 1:8,

  dose = c(
    1, 2, 2, 3,
    3, 4, 4, 3
  ),

  dlt = c(
    0, 0, 0, 0,
    0, 1, 0, 0
  ),

  followup = c(
    28, 28, 21, 14,
    28, 10, 7, 28
  )
)

data

Calculate TITE Weights

data$weight <- ifelse(
  data$dlt == 1,
  1,
  pmin(data$followup / tau, 1)
)

data

The resulting weights are:

data$weight

# 1.00
# 1.00
# 0.75
# 0.50
# 1.00
# 1.00
# 0.25
# 1.00

Define the CRM Model

tox_prob <- function(theta, dose) {

  skeleton[dose]^exp(theta)

}

For example, the Dose 4 probability is:

tox_prob(
  theta = 0.106,
  dose = 4
)

# approximately 0.262

Define the Weighted Log-Likelihood

log_likelihood <- function(theta, data) {

  p <- tox_prob(
    theta,
    data$dose
  )

  sum(
    data$weight *
      (
        data$dlt * log(p) +
        (1 - data$dlt) * log(1 - p)
      )
  )
}

The weighting enters directly into the likelihood.

A fully observed DLT-free patient contributes weight 1. A patient with only 25% of the planned follow-up contributes weight 0.25.

Define the Normal Prior

log_prior <- function(theta) {

  dnorm(
    theta,
    mean = 0,
    sd = sigma,
    log = TRUE
  )

}

Construct the Log Posterior

log_posterior <- function(theta, data) {

  log_likelihood(theta, data) +
    log_prior(theta)

}

The posterior is proportional to:

$$ \exp\{ \log L(\theta)+\log\pi(\theta) \} $$

Numerically Integrate the Posterior

One simple approach is to evaluate the posterior on a sufficiently fine grid of possible \(\theta\) values.

theta_grid <- seq(
  -3,
  1.5,
  length.out = 18001
)

log_post <- sapply(
  theta_grid,
  log_posterior,
  data = data
)

post <- exp(
  log_post -
    max(log_post)
)

post <- post /
  sum(post)

posterior_mean <- sum(
  theta_grid * post
)

posterior_mean

For the illustrative dataset:

posterior_mean

# approximately 0.106

Calculate Posterior DLT Probabilities

posterior_probs <- tox_prob(
  posterior_mean,
  1:length(skeleton)
)

posterior_probs

The resulting probabilities are approximately:

# Dose 1   0.0358
# Dose 2   0.0773
# Dose 3   0.1671
# Dose 4   0.2622
# Dose 5   0.4116

Select the Dose Closest to the Target

next_dose <- which.min(
  abs(
    posterior_probs -
      target
  )
)

next_dose

# 4

Therefore, the model recommends Dose 4 in this illustrative example.

Implement a No-Skipping Restriction

Suppose the current dose is Dose 2. Even if the unrestricted CRM recommends Dose 4, a no-skipping rule can restrict the next dose to Dose 3.

current_dose <- 2

allowed_doses <- c(
  current_dose - 1,
  current_dose,
  current_dose + 1
)

allowed_doses <- allowed_doses[
  allowed_doses >= 1 &
  allowed_doses <= length(skeleton)
]

restricted_dose <-
  allowed_doses[
    which.min(
      abs(
        posterior_probs[
          allowed_doses
        ] - target
      )
    )
  ]

restricted_dose

This separates two concepts:

  • Statistical recommendation: what the model would choose.
  • Operational assignment: what the prespecified safety rules permit.

A Simple TITE-CRM Algorithm in R

tite_crm <- function(
  data,
  skeleton,
  target = 0.25,
  tau = 28,
  sigma = 1.34
) {

  data$weight <- ifelse(
    data$dlt == 1,
    1,
    pmin(data$followup / tau, 1)
  )

  tox_prob <- function(
    theta,
    dose
  ) {
    skeleton[dose]^exp(theta)
  }

  log_likelihood <- function(
    theta
  ) {

    p <- tox_prob(
      theta,
      data$dose
    )

    sum(
      data$weight *
        (
          data$dlt * log(p) +
          (1 - data$dlt) *
          log(1 - p)
        )
    )
  }

  theta_grid <- seq(
    -3,
    1.5,
    length.out = 18001
  )

  log_post <- sapply(
    theta_grid,
    function(theta) {
      log_likelihood(theta) +
        dnorm(
          theta,
          0,
          sigma,
          log = TRUE
        )
    }
  )

  post <- exp(
    log_post -
      max(log_post)
  )

  post <- post /
    sum(post)

  theta_hat <-
    sum(
      theta_grid * post
    )

  fitted_probs <-
    tox_prob(
      theta_hat,
      1:length(skeleton)
    )

  recommended_dose <-
    which.min(
      abs(
        fitted_probs -
          target
      )
    )

  list(
    posterior_mean = theta_hat,
    toxicity_probabilities =
      fitted_probs,
    recommended_dose =
      recommended_dose
  )
}

Run the Example

result <- tite_crm(
  data = data,
  skeleton = skeleton,
  target = 0.25
)

result

The important outputs are:

result$posterior_mean

result$toxicity_probabilities

result$recommended_dose

Interpreting the R Output

The posterior mean summarizes the updated dose-toxicity curve.

The toxicity probability vector gives the estimated DLT probability for each dose.

The recommended dose identifies the dose closest to the target toxicity probability.

For the worked example:

$$ \hat{\theta}\approx0.106 $$

and:

$$ \hat{\mathbf p} \approx (0.036,0.077,0.167,0.262,0.412) $$

so Dose 4 is closest to:

$$ \phi=0.25 $$

Calculating the Information Contribution of Each Patient

A useful diagnostic is to display each patient's weight.

Patient Dose Outcome Follow-Up Weight Interpretation
1 1 No DLT 28 days 1.00 Complete information
2 2 No DLT 28 days 1.00 Complete information
3 2 No DLT 21 days 0.75 Substantial partial information
4 3 No DLT 14 days 0.50 Half-window information
5 3 No DLT 28 days 1.00 Complete information
6 4 DLT 10 days 1.00 Observed event
7 4 No DLT 7 days 0.25 Early partial information
8 3 No DLT 28 days 1.00 Complete information

Why Follow-Up Time Is Different From Calendar Time

TITE-CRM is driven by patient-level follow-up, not simply by the number of days since the trial started.

For example, two patients could enter the study on the same calendar day but have different observation histories because of treatment delays or different event timing.

The statistical weight should therefore be calculated from the prespecified follow-up definition used by the protocol.

What Counts as Follow-Up?

The protocol should define what constitutes valid observation time.

Potential considerations include:

  • Date of first treatment
  • Date of DLT onset
  • Date of treatment discontinuation
  • Date of withdrawal
  • Date of loss to follow-up
  • Date of data cutoff

The exact handling of treatment discontinuation or loss to follow-up is a protocol-level decision.

Do not automatically equate "lost to follow-up" with "DLT-free." The statistical treatment of incomplete follow-up must be specified in advance and should reflect the endpoint definition and assumptions underlying the design.

Interim Dose Assignment Timing

TITE-CRM does not imply that the model must be updated continuously after every minute of follow-up.

The protocol can specify decision times based on:

  • Each new patient
  • Each cohort
  • Scheduled interim reviews
  • Availability of new DLT information

The operational schedule should be selected together with the accrual plan and safety monitoring process.

Accrual Rate Matters

The benefit of TITE-CRM depends heavily on how quickly patients enter the trial.

Suppose the DLT window is 28 days.

If a new patient enters every 30 days, there is little need for a time-to-event design because the previous patient will generally have completed follow-up.

If a new patient enters every three days, however, many patients may be simultaneously pending.

TITE-CRM becomes much more relevant in the second situation.

Accrual and Follow-Up Timeline

Patient Enrollment Day Decision-Relevant Follow-Up at Day 21
1 0 21 days
2 3 18 days
3 6 15 days
4 9 12 days
5 12 9 days
6 15 6 days
7 18 3 days

At calendar day 21, every patient contributes some information.

A conventional approach that required complete follow-up might still have only one fully evaluable patient.

The Price of Staggered Enrollment

The advantage of faster accrual is accompanied by greater dependence on the statistical model.

At any interim decision, some patients may have very little follow-up.

Consequently, the design should be calibrated under realistic combinations of:

  • True dose-toxicity curves
  • DLT probabilities
  • DLT event-time distributions
  • Patient accrual rates
  • Dropout rates
  • Starting dose
  • Cohort size
  • Safety restrictions

Simulation Is Essential

Unlike a simple fixed-sample calculation, a TITE-CRM design has many interacting components.

A simulation study can estimate:

  • Probability of selecting the true target dose
  • Probability of overdosing patients
  • Average number of patients treated at each dose
  • Average trial duration
  • Probability of early stopping
  • Number of dose escalations
  • Number of dose de-escalations
  • Probability of selecting a dose above the target
  • Probability of selecting a dose below the target
Key point: The final TITE-CRM design should be evaluated by simulation rather than judged only from the mathematical form of the model.

Example Simulation Scenario

Suppose the true DLT probabilities are:

$$ (0.03,0.08,0.16,0.25,0.40) $$

The target is 25%, so Dose 4 is the true target dose.

A simulation can repeatedly generate trials under this toxicity curve and record how frequently the design selects Dose 4.

Generating a DLT Outcome

For a simulated patient treated at Dose \(j\), generate:

$$ Y_i\sim\operatorname{Bernoulli}(p_j) $$

If \(Y_i=1\), generate a DLT time from the assumed event-time distribution.

If \(Y_i=0\), the patient remains DLT-free through the current follow-up time.

The simulation then advances calendar time according to the planned accrual process.

Operating Characteristics

A simulation report might include:

Metric Purpose
Probability of selecting Dose 4 Measures target-dose selection accuracy
Probability of selecting an overly toxic dose Measures safety
Mean number treated at target dose Measures dose allocation efficiency
Mean trial duration Measures operational efficiency
Mean number of DLTs Measures patient exposure to toxicity
Probability of stopping early Measures stopping-rule behavior

Common TITE-CRM Design Choices

Several choices must be made before the study begins.

  • Target toxicity probability
  • Number and spacing of dose levels
  • Skeleton probabilities
  • Working dose-toxicity model
  • Prior distribution
  • Follow-up weighting function
  • Starting dose
  • Cohort size
  • Accrual schedule
  • Maximum sample size
  • No-skipping restrictions
  • Overdose-control rules
  • Minimum follow-up requirements
  • Stopping criteria
  • Final MTD selection rule

Choosing the Starting Dose

The starting dose should not simply be the dose that happens to be closest to the target according to the skeleton.

The initial dose should reflect the total clinical context, including preclinical evidence, prior human experience, pharmacology, and safety considerations.

The statistical model should not be treated as a substitute for clinical judgment about a safe starting dose.

Choosing the Skeleton

The skeleton should represent plausible prior toxicity probabilities across the dose range.

For example:

$$ (0.05,0.10,0.20,0.30,0.45) $$

implies increasing toxicity with increasing dose.

The target of 25% lies between the skeleton probabilities for Dose 3 and Dose 4.

This is a sensible structure for an illustrative example because the dose range contains levels on both sides of the target.

Why the Dose Range Matters

Suppose the skeleton were:

$$ (0.01,0.02,0.04,0.07,0.10) $$

If the true target is 25%, every dose is far below the target.

The trial would have difficulty identifying a dose near the target because the available dose range is poorly positioned.

Conversely, if every skeleton probability were far above the target, the dose range could be excessively toxic.

Design principle: The dose range must contain enough plausible toxicity probabilities around the target to make meaningful dose-finding possible.

Model Misspecification

The CRM model is a working model.

The true dose-toxicity relationship may not exactly follow:

$$ p_j^{\exp(\theta)} $$

This is one reason simulation studies should consider multiple possible true toxicity curves rather than only the curve assumed by the working model.

For example, simulations might consider scenarios in which:

  • Toxicity increases slowly at low doses and rapidly at high doses.
  • Toxicity increases rapidly at low doses and then plateaus.
  • The skeleton is approximately correct.
  • The skeleton is substantially misspecified.

The Long-Memory Property

CRM methods can exhibit a phenomenon sometimes described as long memory: early observations can continue to influence the posterior model throughout the trial.

This is a consequence of the model-based approach.

A few early DLTs can substantially affect the posterior curve, particularly when the sample size is small.

This is not necessarily a flaw, but it highlights the importance of:

  • Prior calibration
  • Starting-dose selection
  • Safety restrictions
  • Simulation under alternative scenarios

Pending Patients and Safety Decisions

One of the most important operational questions is whether the trial should allow escalation when many patients are still incompletely observed.

A statistical model can incorporate partial information, but the protocol may still impose additional restrictions.

For example, a design might require:

$$ P( p_j>\phi_{\mathrm{over}} \mid D )

before assigning another patient to Dose \(j\).

Alternatively, a minimum amount of cumulative follow-up may be required before certain escalation decisions.

TITE-CRM and Cohort Size

Although CRM can operate patient-by-patient, cohorts can also be used.

For cohort size \(m\), all patients in the cohort may be assigned the same dose while the model is updated after the cohort's available information has been incorporated.

The choice of cohort size affects:

  • How frequently the model is updated
  • How quickly dose changes occur
  • Operational complexity
  • Patient exposure to a potentially incorrect dose

A cohort size of one provides the most frequent potential reassessment, while larger cohorts reduce the number of dose-assignment decisions.

Interim Decision Example

Suppose the current posterior estimates are:

$$ (0.036,0.077,0.167,0.262,0.412) $$

and the target is:

$$ \phi=0.25 $$

Dose 4 is closest to the target.

If the current dose is Dose 3 and no-skipping is required, the model may permit escalation from Dose 3 to Dose 4.

If the current dose is Dose 2, the same model recommendation could be constrained to Dose 3.

A More Conservative Decision Rule

Some CRM implementations use a conservative rule that selects the highest dose whose estimated toxicity probability does not exceed the target.

That rule can be expressed as:

$$ d_{\text{next}} = \max \left\{ j: p_j(\hat{\theta})\le\phi \right\} $$

For the illustrative posterior probabilities:

$$ (0.036,0.077,0.167,0.262,0.412) $$

the highest dose not exceeding 25% would be Dose 3.

This differs from the closest-to-target rule, which selects Dose 4.

Important: The dose-selection rule must be prespecified. "Closest to target" and "highest dose below target" are different statistical rules and can lead to different assignments.

Posterior Probability of Excessive Toxicity

A Bayesian model also permits probability-based safety summaries.

For dose \(j\), investigators may calculate:

$$ P\left( p_j(\theta)>\phi_{\mathrm{over}} \mid D \right) $$

This asks: Given the current data, what is the posterior probability that this dose's true DLT probability exceeds the predefined overdose threshold?

Such probabilities can be incorporated into prespecified safety rules.

Posterior Probability Calculation

Using the posterior grid:

overdose_limit <- 0.35

dose4_probs <-
  tox_prob(
    theta_grid,
    dose = 4
  )

posterior_overdose_prob <-
  sum(
    post[
      dose4_probs >
        overdose_limit
    ]
  )

posterior_overdose_prob

The resulting quantity can be compared with a prespecified safety cutoff.

The exact cutoff should be selected during design development and evaluated through simulation.

Handling Multiple DLT Types

The basic TITE-CRM formulation assumes a clearly defined binary DLT endpoint.

Real studies may have multiple toxicity events or competing safety outcomes.

The protocol therefore needs a precise definition of what constitutes a DLT.

Examples include:

  • Grade-specific toxicities
  • Clinically significant laboratory abnormalities
  • Hospitalization-related events
  • Treatment interruptions
  • Specific organ toxicities

The statistical model should use the prespecified composite DLT definition rather than changing the definition during the trial.

Time-to-Event CRM Is Not a Survival Model

The term "time-to-event" can be misleading.

TITE-CRM uses the timing of toxicity information to determine how much information is available during the dose-finding process.

It does not necessarily mean that the design is fitting a full survival model for the DLT time distribution.

The basic weighted-likelihood implementation is better understood as a time-weighted extension of CRM.

Do not confuse TITE-CRM with a conventional time-to-event analysis. The primary objective remains dose finding. The follow-up time is used to handle incompletely observed DLT outcomes during adaptive dose assignment.

Protocol Specification

A TITE-CRM protocol should specify at minimum:

  • Primary DLT definition
  • DLT observation window
  • Target toxicity probability
  • Dose levels
  • Starting dose
  • CRM skeleton
  • Working dose-toxicity model
  • Prior distribution
  • Follow-up weighting function
  • Patient or cohort size
  • Accrual assumptions
  • Dose-escalation rules
  • Dose-de-escalation rules
  • Dose-skipping restrictions
  • Overdose-control criteria
  • Minimum follow-up requirements
  • Stopping rules
  • Final MTD selection rule
  • Handling of withdrawals and missing data

Statistical Analysis Plan Considerations

The Statistical Analysis Plan should make the adaptive algorithm reproducible.

A useful TITE-CRM SAP section should state exactly how:

  • DLT events are coded
  • Follow-up time is calculated
  • Weights are assigned
  • Pending patients enter the likelihood
  • The posterior distribution is calculated
  • The posterior summary is selected
  • Doses are ranked relative to the target
  • Safety restrictions are applied
  • Stopping criteria are evaluated
  • The final MTD is selected

Common Mistakes

  1. Treating a pending patient as fully DLT-free. A patient who has completed only part of the observation window should not automatically receive a full weight of one.
  2. Ignoring partial follow-up entirely. That defeats the primary purpose of TITE-CRM.
  3. Treating a DLT as partial information. Once a DLT is observed, the event itself is known and receives the event contribution specified by the design.
  4. Changing the weighting function after seeing the data. The weighting scheme is part of the design.
  5. Ignoring accrual rate. The practical value of TITE-CRM depends heavily on how quickly patients enter the trial relative to the DLT window.
  6. Ignoring event-time distributions. Simulation should consider plausible timing of DLTs, not just their overall probability.
  7. Using an unsafe unrestricted CRM. Practical safety restrictions should be evaluated during design development.
  8. Assuming the skeleton is the final toxicity curve. The skeleton is a starting point for the working model.
  9. Using the model without simulation. The operating characteristics of the complete adaptive design should be evaluated under realistic scenarios.
  10. Confusing the target toxicity level with an observed toxicity rate. The target is a design parameter; it is not necessarily the observed DLT rate at the end of the trial.

A Practical TITE-CRM Workflow

1
Define the DLT endpoint.
2
Define the DLT observation window.
3
Specify the target toxicity probability.
4
Select the dose levels and starting dose.
5
Construct and calibrate the CRM skeleton.
6
Specify the working dose-toxicity model.
7
Specify the prior distribution.
8
Specify the time-to-event weighting function.
9
Specify cohort size and accrual assumptions.
10
Specify escalation, de-escalation, and no-skipping rules.
11
Specify overdose-control and stopping rules.
12
Evaluate the complete design through simulation.
13
Prespecify the final MTD selection rule.

Worked Example Summary

Component Value
Design TITE-CRM
Number of doses 5
DLT observation window 28 days
Target toxicity probability 25%
Skeleton 0.05, 0.10, 0.20, 0.30, 0.45
Working model \(p_j(\theta)=p_j^{\exp(\theta)}\)
Prior \(\theta\sim N(0,1.34^2)\)
Patients at interim analysis 8
Observed DLTs 1
Posterior mean \(\theta\) Approximately 0.106
Posterior DLT probability, Dose 1 Approximately 3.6%
Posterior DLT probability, Dose 2 Approximately 7.7%
Posterior DLT probability, Dose 3 Approximately 16.7%
Posterior DLT probability, Dose 4 Approximately 26.2%
Posterior DLT probability, Dose 5 Approximately 41.2%
Target dose by closest-to-target rule Dose 4

The Most Important Concept

The most important conceptual difference between conventional CRM and TITE-CRM is not the dose-toxicity model.

The model can be essentially the same.

The critical change is that TITE-CRM allows partially observed patients to contribute information before their complete DLT observation windows have ended.

A patient with a known DLT contributes the event.

A patient who remains DLT-free contributes information proportional to the prespecified amount of follow-up under the basic weighting approach.

The model then combines these contributions with the prior to estimate the current dose-toxicity relationship.

$$ \boxed{ \text{Prior} + \text{Observed DLTs} + \text{Partial Follow-Up} \longrightarrow \text{Posterior Dose-Toxicity Curve} } $$

The next dose is then selected using the prespecified target-toxicity rule, subject to the trial's safety constraints.

Bottom line: Time-to-Event CRM is a model-based Bayesian dose-finding design that extends the continual reassessment method to settings in which DLTs may occur late in the observation window. Rather than waiting for every patient's complete follow-up, TITE-CRM incorporates partial follow-up through prespecified weights. This permits staggered enrollment and more frequent dose reassessment. The design combines a dose-toxicity model, prior distribution, weighted likelihood, target toxicity level, and prespecified safety rules. Because the behavior of the complete adaptive design depends on toxicity probabilities, event-time distributions, accrual, prior assumptions, and operational restrictions, simulation is an essential part of design development.

References

Cheung, Y.K. & Chappell, R. (2000). Sequential designs for Phase I clinical trials with late-onset toxicities. Biometrics, 56(4), 1177–1182.

O'Quigley, J., Pepe, M. & Fisher, L. (1990). Continual reassessment method: a practical design for Phase I clinical trials in cancer. Biometrics, 46, 33–48.

Lee, S.M. & Cheung, Y.K. (2011). Model calibration in the continual reassessment method. Clinical Trials, 8, 584–591.

Cheung, Y.K. (2011). Stochastic approximation and modern model-based designs for dose-finding clinical trials. Statistical Science, 25(4), 476–493.

Iasonos, A. & O'Quigley, J. (2014). Adaptive dose-finding studies: a review of model-based designs. Clinical Trials.

Ji, Y. & Wang, S.J. (2013). Modified toxicity probability interval design: a safer and more reliable method for dose-finding trials. Clinical Trials, 10(3), 321–331.

Babb, J., Rogatko, A. & Zacks, S. (1998). Cancer Phase I clinical trials: efficient dose escalation with overdose control. Statistics in Medicine, 17, 1103–1120.