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.
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:
The patient contributes approximately one quarter of the information that would be available after completing the entire observation window.
CRM Before TITE-CRM
To understand TITE-CRM, first consider the standard CRM.
Suppose there are \(J\) dose levels:
For each dose level, the design specifies an initial estimate of the probability of DLT.
These values form the CRM skeleton:
The skeleton should be monotonically increasing:
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:
For example, suppose the target toxicity probability is 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.
The CRM Dose-Toxicity Model
A common one-parameter CRM model uses the power model:
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:
the model reproduces the skeleton:
If \(\theta>0\), the probabilities are shifted downward because
\(\exp(\theta)>1\) and \(0
If \(\theta<0\), the probabilities are shifted upward.
CRM uses a prior distribution for the unknown parameter.
For this tutorial, suppose:
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.
The fundamental CRM model remains the same.
The key modification is the treatment of incomplete follow-up.
For patient \(i\), define:
For the basic TITE-CRM weighting scheme:
The no-DLT weight is typically bounded at 1:
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.
The weight is therefore a measure of the amount of available follow-up under
this simple weighting approach.
Suppose patient \(i\) is treated at dose \(j\), and let the modeled DLT
probability at that dose be:
For a conventional binary outcome, the likelihood contribution would be:
TITE-CRM raises this contribution to the patient's follow-up weight:
Equivalently:
The total weighted likelihood is:
The posterior is proportional to the prior multiplied by this likelihood:
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:
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.
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:
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.
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:
This can substantially reduce trial duration when toxicity outcomes are late
relative to the desired enrollment rate.
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:
The investigators specify the following skeleton:
The working model is:
and the prior is:
Suppose eight patients have entered the study.
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.
For Patient 3:
For Patient 4:
For Patient 7:
Patient 6 experienced a DLT, so:
The complete vector of weights is therefore:
The five skeleton probabilities are:
For any value of \(\theta\), the modeled probabilities are:
Patient 1 was treated at Dose 1 and had no DLT with complete follow-up.
His contribution is:
Patient 3 was treated at Dose 2 and had no DLT after 21 days.
His contribution is:
Patient 4 was treated at Dose 3 and had no DLT after 14 days.
His contribution is:
Patient 6 experienced a DLT at Dose 4.
His contribution is:
Patient 7 was treated at Dose 4 and had no DLT after 7 days.
His contribution is:
Multiplying all eight contributions gives the complete weighted likelihood.
The posterior distribution is proportional to:
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:
For the illustrative dataset above, numerical integration gives approximately:
This is the posterior mean used in the subsequent dose-toxicity calculation.
Using:
the fitted DLT probabilities are approximately:
The target toxicity probability is:
Dose 4 has an estimated DLT probability of approximately 26.2%, which is very
close to the target.
The standard CRM decision rule is:
For the example:
Therefore:
Under this illustrative model and dataset, Dose 4 is the dose whose posterior
estimated DLT probability is closest to the 25% target.
This provides a useful way to understand the weighting mechanism.
At seven days:
After completing the full observation window:
The patient's contribution would therefore become:
rather than:
The statistical information contributed by the patient increases as follow-up
accumulates.
Suppose Patient 7 develops a DLT on day 12.
The patient is no longer treated as a pending DLT-free observation.
Instead:
The contribution becomes:
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.
The conventional CRM uses a complete binary outcome:
TITE-CRM supplements this with:
and therefore uses:
within the weighted likelihood.
The conceptual transition is:
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:
Therefore, TITE-CRM is particularly attractive when the toxicity observation
window is long relative to the desired accrual interval.
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.
The simple weight:
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.
Instead of a linear weight, one can define:
where \(g(\cdot)\) reflects an assumed pattern of information accumulation.
A simple linear approach is:
More sophisticated approaches can incorporate assumptions about the hazard or
timing distribution of DLTs.
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:
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.
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.
Whether dose skipping is allowed is a design choice and must be explicitly
specified.
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:
and investigators define an overdose threshold:
A possible safety criterion could require:
where \(c\) is a prespecified probability cutoff.
This is conceptually different from simply choosing the dose with toxicity
probability closest to the target.
One of the most important practical questions in TITE-CRM is how pending
patients are represented.
Suppose the trial currently has:
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.
TITE-CRM is called a continual reassessment method because the model is
re-estimated as new information becomes available.
For example:
A TITE-CRM protocol should include explicit stopping rules.
Examples include:
The exact stopping criteria are design-specific.
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:
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.
There are three important components:
The Bayesian updating mechanism is therefore:
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.
Suppose a DLT occurs at Dose 4.
The DLT contributes a complete event:
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.
Suppose a patient at Dose 4 remains DLT-free for 14 days.
The weight is:
The contribution is:
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.
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:
This is one reason CRM can be more statistically efficient than independent
dose-by-dose rules.
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 and TITE-BOIN both address the same operational problem:
DLT outcomes may not be immediately available.
Their underlying statistical approaches differ.
Both approaches can be useful, but they should not be treated as equivalent
algorithms.
The following R code establishes the example design.
The resulting weights are:
For example, the Dose 4 probability is:
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.
The posterior is proportional to:
One simple approach is to evaluate the posterior on a sufficiently fine grid
of possible \(\theta\) values.
For the illustrative dataset:
The resulting probabilities are approximately:
Therefore, the model recommends Dose 4 in this illustrative example.
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.
This separates two concepts:
The important outputs are:
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:
and:
so Dose 4 is closest to:
A useful diagnostic is to display each patient's weight.
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.
The protocol should define what constitutes valid observation time.
Potential considerations include:
The exact handling of treatment discontinuation or loss to follow-up is a
protocol-level decision.
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:
The operational schedule should be selected together with the accrual plan and
safety monitoring process.
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.
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 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:
Unlike a simple fixed-sample calculation, a TITE-CRM design has many interacting
components.
A simulation study can estimate:
Suppose the true DLT probabilities are:
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.
For a simulated patient treated at Dose \(j\), generate:
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.
A simulation report might include:
Several choices must be made before the study begins.
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.
The skeleton should represent plausible prior toxicity probabilities across
the dose range.
For example:
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.
Suppose the skeleton were:
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.
The CRM model is a working model.
The true dose-toxicity relationship may not exactly follow:
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:
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:
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:
before assigning another patient to Dose \(j\).
Alternatively, a minimum amount of cumulative follow-up may be required before
certain escalation decisions.
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:
A cohort size of one provides the most frequent potential reassessment, while
larger cohorts reduce the number of dose-assignment decisions.
Suppose the current posterior estimates are:
and the target is:
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.
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:
For the illustrative posterior probabilities:
the highest dose not exceeding 25% would be Dose 3.
This differs from the closest-to-target rule, which selects Dose 4.
A Bayesian model also permits probability-based safety summaries.
For dose \(j\), investigators may calculate:
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.
Using the posterior grid:
The resulting quantity can be compared with a prespecified safety cutoff.
The exact cutoff should be selected during design development and evaluated
through simulation.
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:
The statistical model should use the prespecified composite DLT definition
rather than changing the definition during the trial.
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.
A TITE-CRM protocol should specify at minimum:
The Statistical Analysis Plan should make the adaptive algorithm reproducible.
A useful TITE-CRM SAP section should state exactly how:
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.
The next dose is then selected using the prespecified target-toxicity rule,
subject to the trial's safety constraints.
Cheung, Y.K. & Chappell, R. (2000).
Sequential designs for Phase I clinical trials with late-onset
toxicities.
Biometrics, 56(4), 1177–1182.The Bayesian Prior
What Changes in TITE-CRM?
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
How the Weighted Likelihood Works
Why a DLT Gets Weight 1
Why a DLT-Free Patient Gets a Partial Weight
Staggered Patient Entry
A Complete Worked Example
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
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
Patient-Level Data at an Interim Look
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
Step 1: Calculate the Follow-Up Weights
Step 2: Write the Dose-toxicity Model
Step 3: Construct the Weighted Likelihood
Step 4: Estimate the Posterior Parameter
Step 5: Convert the Posterior Parameter into DLT Probabilities
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%
Step 6: Select the Dose Closest to the Target
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
What If Patient 7 Had Completed 28 Days?
What If Patient 7 Develops a DLT?
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
Why TITE-CRM Can Shorten a Trial
But TITE-CRM Does Not Eliminate Uncertainty
Time-to-Event Information Is Not Necessarily Linear
Generalized Follow-Up Weights
Safety Restrictions Are Essential
No Dose Skipping
Overdose Control
Handling Pending Patients
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
Repeated Model Updating
Stopping Rules
Final MTD Selection
TITE-CRM Is a Bayesian Adaptive Design
Interpreting the Posterior Dose-Toxicity Curve
What Happens After a DLT?
What Happens After a DLT-Free Patient?
Why the Model Borrows Information Across Doses
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
TITE-CRM vs. TITE-BOIN
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
R Implementation: Define the 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
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)
}
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)
)
)
}
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)
}
Numerically Integrate the Posterior
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
posterior_mean
# approximately 0.106
Calculate Posterior DLT Probabilities
posterior_probs <- tox_prob(
posterior_mean,
1:length(skeleton)
)
posterior_probs
# 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
Implement a No-Skipping Restriction
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
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
result$posterior_mean
result$toxicity_probabilities
result$recommended_dose
Interpreting the R Output
Calculating the Information Contribution of Each Patient
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
What Counts as Follow-Up?
Interim Dose Assignment Timing
Accrual Rate Matters
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
The Price of Staggered Enrollment
Simulation Is Essential
Example Simulation Scenario
Generating a DLT Outcome
Operating Characteristics
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
Choosing the Starting Dose
Choosing the Skeleton
Why the Dose Range Matters
Model Misspecification
The Long-Memory Property
Pending Patients and Safety Decisions
TITE-CRM and Cohort Size
Interim Decision Example
A More Conservative Decision Rule
Posterior Probability of Excessive Toxicity
Posterior Probability Calculation
overdose_limit <- 0.35
dose4_probs <-
tox_prob(
theta_grid,
dose = 4
)
posterior_overdose_prob <-
sum(
post[
dose4_probs >
overdose_limit
]
)
posterior_overdose_prob
Handling Multiple DLT Types
Time-to-Event CRM Is Not a Survival Model
Protocol Specification
Statistical Analysis Plan Considerations
Common Mistakes
A Practical TITE-CRM Workflow
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
References
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.