<figure><img src="image_1.png" /></figure>

HAL
open science

# GUIP1: a R package for dose escalation strategies in phase I cancer clinical trials

Derek Dinart, Julien Fraisse, Diego Tosi, Audrey Mauguen, Célia Touraine, Sophie Gourgou, Marie Cécile Le Deley, Carine A. Bellera, Caroline Mollevi

► To cite this version:

Derek Dinart, Julien Fraisse, Diego Tosi, Audrey Mauguen, Célia Touraine, et al.. GUIP1: a R package for dose escalation strategies in phase I cancer clinical trials. BMC Medical Informatics and Decision Making, 2020, 20 (1), 10.1186/s12911-020-01149-3. hal-02889984

HAL Id: hal-02889984

https://hal.science/hal-02889984

Submitted on 6 Jul 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SOFTWARE

Open Access

<figure><img src="image_3.png" /></figure>

Check for
updates

# GUIP1: a R package for dose escalation strategies in phase I cancer clinical trials

D. Dinart¹†, J. Fraisse<fnref n="2" />†, D. Tosi<fnref n="2" />, A. Mauguen<fnref n="3" />, C. Touraine<fnref n="2" />, S. Gourgou<fnref n="2" />, M. C. Le Deley⁴,⁵, C. Bellera¹ and C. Mollevi<fnref n="2" />,⁶*

<figure><img src="image_4.png" /></figure>

## Abstract

**Background:** The main objective of phase I cancer clinical trials is to identify the maximum tolerated dose, usually defined as the highest dose associated with an acceptable level of severe toxicity during the first cycle of treatment. Several dose-escalation designs based on mathematical modeling of the dose-toxicity relationship have been developed. The main ones are: the continual reassessment method (CRM), the escalation with overdose control (EWOC) method and, for late-onset and cumulative toxicities, the time-to-event continual reassessment method (TITE-CRM) and the time-to-event escalation with overdose control (TITE-EWOC) methods. The objective of this work was to perform a user-friendly R package that combines the latter model-guided adaptive designs.

**Results:** GUIP1 is an R Graphical User Interface for dose escalation strategies in Phase 1 cancer clinical trials. It implements the CRM (based on Bayesian or maximum likelihood estimation), EWOC and TITE-CRM methods using the dfCRM and bCRM R packages, while the TITE-EWOC method has been specifically developed. The program is built using the TCL/TK programming language, which can be compiled via R software libraries (tcltk, tkplot, tcltk2). GUIP1 offers the possibility of simulating and/or conducting and managing phase I clinical trials in real-time using file management options with automatic backup of study and/or simulation results.

**Conclusions:** GUIP1 is implemented using the software R, which is widely used by statisticians in oncology. This package simplifies the use of the main model-based dose escalation methods and is designed to be fairly simple for beginners in R. Furthermore, it offers multiple possibilities such as a full traceability of the study. By including multiple innovative adaptive methods in a free and user-friendly program, we hope that GUIP1 will promote and facilitate their use in designing future phase I cancer clinical trials.

**Keywords:** Phase 1 clinical trial, Graphical user interface, Dose escalation design

## Background

The development of cancer drugs in a clinical setting requires a well-codified procedure. Phase I studies of a new treatment are usually the first to involve human subjects, and their aim is to select doses according to

acceptable toxicity. The drug efficacy is then preliminarily tested in a phase II trial. Finally, a phase III trial compares the safety and efficacy of the new treatment against the current standard treatment prior to its approval for general use.

First-in-human (FIH) phase I studies of cancer drugs are performed on patients for whom no other therapeutic option is available. The main objective of phase I cancer clinical trials is to identify the maximum tolerated dose (MTD) of an experimental drug. Commonly, the MTD is associated with a predefined probability of unacceptable toxicity (20–33%), called the dose-limiting

<footnote>* Correspondence: Caroline.Mollevi@icm.unicancer.fr</footnote>

<footnote>†D. Dinart and J. Fraisse contributed equally to this work.</footnote>

<footnote n="2">²Institut du Cancer Montpellier (ICM), Université de Montpellier, Montpellier, France</footnote>

<footnote n="3">³Institut de Recherche en Cancérologie de Montpellier INSERM U1194, Université de Montpellier, Montpellier, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
toxicity (DLT), which is evaluated in general during the first cycle of treatment. The recommended phase II dose (RP2D) is then suggested for subsequent phase II trials. Most of the available statistical methods used to design phase I clinical trials in oncology have been developed for cytotoxic conventional agents.

Two types of dose-escalation designs have been proposed: algorithm-based designs, such as the standard 3+3 design, and model-based designs. The 3+3 design is the simplest to interpret and implement [1], and it is the most commonly used phase I dose-escalation method. This method defines the MTD as the highest dose at which less than one-third of the treated patients experience intolerable toxicity during the first cycle of treatment. However, the 3+3 method underuses available trial data and is based only on empirical evidence without including statistical considerations; thus, this method allocates more patients to sub-therapeutic doses and presents lower operating characteristics than model-based methods [2, 3].

Model-based methods, being more computationally intensive and complex, are still underused in practice, although their use has gradually increased in recent years [4]. A barrier to their use is the need for software that requires specific skills. They are based on mathematical modeling of the dose-toxicity relationship. Specifically, the dose-toxicity curve is estimated using a parametric mathematical model that determines the probability of a DLT for the next patient to be included in the study, based on the doses and responses obtained for all previous patients. The continual reassessment method with either Bayesian estimation (CRMB) [5] or maximum likelihood estimation (CRML) [6] were the first two adaptive designs proposed. The CRMB and CRML models only incorporate complete observations, leading to a staggered accrual. In certain frameworks, such as the assessment of late-onset toxicities in radiotherapy-based trial, the observation period, and thus the trial duration, increases. Indeed, if one wants to wait for the full assessment of all previous observations period, then it is necessary to temporarily stop accrual, even though eligible patients who have been informed about the trial are immediately available. An extension of the CRM method, called time-to-event (TITE)-CRM [7], takes into account incomplete observations in a sequential reestimation of the initial CRM model. Specifically, the likelihood is weighted by a function of the actual observation time of the patients and the total observation period defined for the study. For example, if a linear function is chosen, assuming 2 months of follow-up out of 4, the weight would be equal to 0.5.

The escalation with overdose control (EWOC) method [8] is the first dose-finding procedure that directly incorporates an ethical constraint to limit the probability of

treating patients with excessively high doses. This constraint is applied to the CRM model via a parameter $\alpha$, which represents the expected proportion of patients treated with doses higher than the MTD. The EWOC method approaches the MTD as quickly as possible while keeping the expected proportion of overdosed patients below the $\alpha$ value.

Note that others robust versions of the CRM method were also developed such Bayesian model averaging CRM [9] which proposes to model several parallel CRM with different prespecified toxicity probabilities or Bayesian data augmentation CRM [10] which addresses the late-onset toxicity problem. At last, CRM was extended to take into account ordinal toxicity outcomes [11, 12] or two competing outcomes [13].

More recently, a new class of models, the model-assisted designs, combining the simplicity of the algorithm-based designs and the performance of model based designs was introduced [14–16] such the modified toxicity probability interval (mTPI) design [17] or Bayesian optimal interval (BOIN) design [18] for example. Zhou et al. [15] showed that model-based designs, such as CRM, and model-assisted designs, such as BOIN, had similar performances and that we could favor one or the other depending on prior knowledge of the true dose-toxicity curve. However, despite good operating characteristics model-assisted designs are still underutilized.

Several software programs that implement dose-escalation designs, such as np1 [19], R packages: dfcrm, bcrm, and CRM ( [20], SAS macro TITE-CRM (https://sph.umich.edu/ccb/tite-resources.html), and biostatistics software proposed by MD Anderson: https://biostatistics.mdanderson.org/SoftwareDownload/ and the website named Trial Designs (https://www.trialdesign.org/), are available to conduct a phase I trial. Here, we propose a free and user-friendly R graphical interface, named GUIP1, that combines different model-guided adaptive designs (CRMB, CRML, EWOC, TITE-CRM, TITE-EWOC) for simulating and conducting phase I cancer clinical trials. This software aims to facilitate the design and analysis of phase I trials by offering an easy-to-use interface for entering the trial's characteristics, a simulation setting to assess each model's performance for the given trial, and a graphical output of the results. We used published R packages (dfcrm and bcrm) and implemented the TITE-EWOC method.

## Implementation

Model-based dose escalation designs in phase I clinical trials are generally based on the following two stages (for additional details, see [5–8, 21]).

First, clinical inputs are required, such as the number of dose levels and the first dose level to be investigated. Similarly, the DLT must be carefully defined, including
the observation period. Based on previous publications on treatments used in similar clinical settings, the toxicity target level and the prior toxicity probabilities must be established for each dose level. Finally, a stopping rule needs to be chosen.

Second, the dose-toxicity relationship is modeled with the underlying assumption that the probability of DLT monotonically increases with dose levels.

## Definitions and notations

Let $n$ be the number of patients included in the study. $Y_1, Y_2, \dots, Y_n$ are the binary random variables defined for any patient $j \in [1:n]$ such that $Y_j = 1$ if patient j has DLT and $Y_j = 0$ otherwise. Let $X$ be the dose variable. We consider the space $V$ of the $k$ dose levels administered: $V = \{x_1, \le x_2 \le \dots x_k\}$, where $x_i$ is the dose associated with the $i^{th}$ level. Let $\Omega_j = \{(x_1, y_1), (x_2, y_2), \dots, (x_{j-1}, y_{j-1})\}$ be the history of the $j-1$ first patients' toxic responses and $\theta$ the target toxicity level.

## Continual reassessment method (CRM)

This method models the probability of toxic response for the $j^{th}$ patient conditionally to the dose $x_i$ by $\psi(x_i, a) = P(Y_j = 1 | X = x_i)$, where $i \in [1..k]$, $j \in [1..n]$, and $a$ is the model parameter. It is assumed that there is an $a_0 > 0$ such that $\psi(x^*, a_0) = \theta$ where $x^*$ is the dose with the target toxicity level. In practice, once $a_0$ obtained, the administered dose corresponds to the dose $x_0 \in V$, such that $\psi(x_0, a_0)$ is the closest to the targeted toxicity level $\theta$. The CRM was initially developed in a Bayesian setting, and subsequently a frequentist approach was proposed [6]. The empiric power model ($\psi(x_0, a_0) = x_0^{\exp(a_0)}$) presented in the original paper is one option of modeling, but it is not the only option. Logistic regression modeling can also be used (one-parameter logistic or two-parameter logistic) with similar operating characteristics, especially when one parameter is fixed and only one parameter is estimated [22].

## Bayesian approach (CRMB)

Before the first inclusion ($j=1$), a prior distribution is chosen for $a_0$ which is denoted by $f(a, \Omega_1)$ (for example, $f(a, \Omega_1) = \exp(-a)$). At the $j^{th}$ patient enrollment ($j > 1$), the distribution of the parameter $a_0$ is reassessed using Bayes' theorem. We obtain the updated posterior $f(a, \Omega_j)$ from

$$f(a, \Omega_j) = \frac{f(a, \Omega_{j-1}) \times \prod_{l=1}^{j-1} [\psi(x_l, a)]^{y_l} [1-\psi(x_l, a)]^{1-y_l}}{\int_0^\infty \prod_{l=1}^{j-1} [\psi(x_l, u)]^{y_l} [1-\psi(x_l, u)]^{1-y_l} du}$$

where

$$\prod_{l=1}^{j-1} [\psi(x_l, a)]^{y_l} [1-\psi(x_l, a)]^{1-y_l} \quad (1)$$

(1) is the current likelihood of the model.

An updated estimate of parameter $a_0$ is given by the conditional expectation of this distribution:

$$\hat{a}_j = E[a|\Omega_j] = \int_0^\infty af(a, \Omega_j)da$$

Then, the probability of DLT for the $j^{th}$ patient at dose $x_i$ is estimated by $\psi(x_i, \hat{a}) = \int_0^\infty \psi(x_i, a) \times f(a, \Omega_j) da$.

## Frequentist approach (CRML)

The aim of the frequentist approach is to address the need for estimating the Bayesian prior information. The parameter $a$ of the distribution is estimated by maximizing the likelihood provided in Eq. (1) above. In contrast to the Bayesian approach, the maximum likelihood method requires an additional condition to get started. This is needed because the likelihood equation has no solution until at least one DLT and one non-DLT have been observed. Therefore, this method generally starts with an algorithm-based approach (for example, 3+3) until the first DLT is observed. The choice of a Bayesian approach versus the frequentist approach has a very little impact on the operating characteristics of the method [22].

## Time-to-event CRM (TITE-CRM)

In both TITE approaches, we choose to use a linear weight. Let $T$ be the planned observation period, and we suppose that patient $j$ was followed until $u_j \le T$. We define the weighting factor $w_j$ for this patient by

$$w_j = \begin{cases} \frac{u_j}{T} & \text{if patient j is still on follow-up without DLT} \\ 1 & \text{if patient j has a DLT} \end{cases}$$

The weighted likelihood of the model would thus simply be written as $\prod_{l=1}^{j-1} [w_l \times \psi(x_l, a)]^{y_l} [1-w_l \times \psi(x_l, a)]^{1-y_l}$.

The weights are now part of the likelihood of the model, which affect the $\hat{a}$-parameter posterior distribution. Note that the TITE-CRM method is only proposed in a Bayesian setting in GUIP1.

## Escalation with overdose control (EWOC)

EWOC is a Bayesian dose-finding design that produces a consistent sequence of doses while controlling for the probability of overdose. This ethical constraint is the main difference between the EWOC and CRM methods. In other words, EWOC ensures that the expected proportion of patients receiving a dose greater than the MTD (overdose) does not exceed a pre-
specified value α, the feasibility bound. The EWOC method selects, at each new inclusion j, the dose level such that the conditional probability of exceeding the MTD, defined by $π_j(x) = P(MTD ≤ x | Ω_j)$, is inferior or equal to the parameter α.

A two-parameter logistic model was initially considered to model the dose-toxicity relationship [8]. By denoting $a = (a_0, a_1)$ as the unknown parameter vector, the probability of toxic response for the $j^{th}$ patient at dose $x_i$ is estimated to be $\psi(x_i, a) = \frac{\exp(a_0+a_1x_i)}{1+\exp(a_0+a_1x_i)}$, where $a_1 \in \mathbb{R}^+$ to ensure that the probability of DLT monotonically increases with dose levels. For additional information on parameter estimation, please see the article by Babb et al. (1998), which includes a re-parametrization of this approach in terms of MTD and the probability of DLT at the initial dose level.

## Time-to-event EWOC (TITE-EWOC)

TITE-EWOC is a hybrid design that introduces the time-to-event approach in the EWOC method. The aim of this design is to enable continuous recruitment and to decrease the dose-finding trial duration without impairing the characteristics of the EWOC design, especially its ability to control against overdose.

## GUIP1 implementation

The different dose escalation methods discussed above are implemented in different libraries in the R software. The dfcrm library allows the CRM method to be applied with the crm function. The choice of CRMB or CRML is made using the “method” parameter. Two dose-toxicity models are proposed: a one-parameter logistic model with fixed intercept and an empirical model. The titecrm function of this library implements the TITE-CRM [7]. The EWOC method [8] can be applied with the bcrm function of the bcrm library, which allows the dose-toxicity relationship to be modeled with several functional forms. We implemented the TITE-EWOC method and integrated it into GUIP1, as it was only available upon request to the authors.

GUIP1 was built using the TCL/TK programming language, which can be compiled via the R software libraries tcltk, tkrplot, and tcltk2. It implements the methods previously described using five main functions: CRML, CRMB, TITE-CRM, EWOC, and TITE-EWOC. Using the GUIP1 package, the interface can be launched using the command GUIP1() (Fig. 1a, b). A long form of the documentation (vignette) was also developed to facilitate the use of GUIP1.

## Results

### Motivating example

A motivating example is used to describe the program. This section describes the required clinical inputs for implementing CRMB (Fig. 1c, d) as well as CRML, TITE-CRM, EWOC, or TITE-EWOC (supplementary Fig. 1a, b, c, d, respectively). Six dose levels {$x_1, ..., x_6$} were chosen, where $x_3$ was the assumed dose level associated with the target toxicity level of 20% (i.e., $\theta = 0.2$). The prior estimates of the probability of patients who would experience a DLT, at each dose level from 1 to 6, were 0.049, 0.111, 0.2, 0.308, 0.423, and 0.534, respectively. It was assumed that patients were followed over 12 weeks for time-to-event designs and over 4 weeks for the other designs. A one-parameter model (empiric or logistic) was chosen in both cases. For the EWOC method, the probability of exceeding the target dose $\alpha$ was set to 0.25. The stopping rule of the study corresponded to the maximum number of patients who could be enrolled in the study and was fixed at 25 (Fig. 1c, d).

## The CRMB and CRML functions

The views obtained for both the CRMB and CRML functions are similar. The “Interactive CRM” tab allows a new study to be started or an ongoing study to be opened (Fig. 1b). Using the “Prior calibration” tab (Fig. 1c), the following different clinical input parameters are defined: number of dose levels, target toxicity level, desired halfwidth of the indifference intervals, dose-toxicity model, parameters related to the chosen model, and the starting dose level (assumed MTD level). The acceptable toxicity probabilities associated with each dose level were estimated by clicking on the “Run” button using the getprior function (package dfcrm) according to the model calibration approach [23]. The acceptable toxicity probabilities can be changed manually by the users. Once the toxicity probabilities are validated, we complete the additional trial parameters: the stopping rule and its related number of patients which will define the maximum cohort size, as well as the study name (“Input parameters” tab, Fig. 1d). Two stopping rules are currently implemented. The first rule stops the inclusions if the maximum sample size is reached, and the second rule causes a stop if a predefined number of patients has been treated to a certain dose level. On the same tab, the “Save study” button generates two files in the working directory: “study_name.crmb” or “study_name.crml” contains different characteristics of the study, and “study_name-crml.Rdata” or “study_name-crmb.Rdata” contains the updated data after each patient assessment. Note that the extensions (.crmb or .crml) differ according to the function. In the “Include” tab (Fig. 2a), the “New patient” button allows patients to be enrolled in the trial. At each new inclusion, the user
<figure><img src="image_5.png" /></figure>

must provide the given dose level and the observed DLT response (0 = no and 1 = yes) or indicate “pending” if appropriate. The default dose level corresponds to the estimated dose closest to the MTD, taking into account the data from previous patients. The “Pending patient” button allows the patient data to be updated (DLT yes/no) once the response is known. After the first patient is included, the study parameters can no longer be modified. A summary of the results of the inclusions can be viewed at any time via the “Results” tab. Two tables summarizing the results per patient and dose level and a graph showing the dose escalation are displayed (Fig. 2b). These different results can be exported via the “Export” button, which creates an .xlsx file in the current directory that contains these two tables as well as the graphic in png format.

Our motivating example treated up to 25 patients. The starting dose was defined as dose level 3 and was given to the first enrolled patient. The CRMB ran sequentially

for each patient. Each new patient was treated at the proposed dose level (closest to the MTD) and, at the end of the observation period, the occurrence (or not) of DLT was recorded. This step was repeated until the study stopped. Results are shown in Fig. 2b and summarized by patient and dose level.

## TITE-CRM

The operation of the function TITE-CRM is similar to that of the CRMB and CRML functions except for the “Input parameters” and “Include” tabs. Indeed, since this method accounts for the time until DLT observation, it is necessary to define the maximum observation period for all patients. This parameter is indicated in the “Observation Window” field of the “Input parameters” tab. Consequently, for each new inclusion, a variable recording the follow-up was added. An update of the follow-up is performed for all patients with incomplete observations. The default value is equal to the maximum
<figure><img src="image_6.png" /></figure>

a

Snapshot for the step to define the dose to be allocated to the 12th patient, after observing only 1 DLT on dose level 5 and 2 DLTs on dose level 6. Note that the DLTs are marked as red points in Fig. 2b.

<figure><img src="image_7.png" /></figure>

b

Snapshot displaying the summary results at the end of the trial (here, after the inclusion of the planned 25 patients). In this example, skipping a dose is allowed. The users can manually decide whether to apply the “no skipping dose” constraint.

<table><caption>**Fig. 2** Conducting a study with GUIP1 / CRMB (motivating example)</caption><thead><tr><th>Level</th><th>Prior</th><th>n.patient</th><th>total.tox</th><th>ptox</th><th>ptoxL</th><th>ptoxU</th></tr></thead><tbody><tr><td>1</td><td>0.049</td><td>0</td><td>0</td><td>0.003</td><td>0.000</td><td>0.023</td></tr><tr><td>2</td><td>0.111</td><td>0</td><td>0</td><td>0.016</td><td>0.002</td><td>0.065</td></tr><tr><td>3</td><td>0.200</td><td>1</td><td>0</td><td>0.047</td><td>0.010</td><td>0.135</td></tr><tr><td>4</td><td>0.308</td><td>3</td><td>0</td><td>0.107</td><td>0.033</td><td>0.231</td></tr><tr><td>5</td><td>0.423</td><td>16</td><td>3</td><td>0.196</td><td>0.084</td><td>0.343</td></tr><tr><td>6</td><td>0.534</td><td>5</td><td>2</td><td>0.305</td><td>0.164</td><td>0.458</td></tr></tbody></table>

observation period. Note that the value of patient follow-up must be smaller than the duration of the observation period. Supplementary Fig. 1a shows the view obtained for the TITE-CRM method.

## EWOC

Similar to the CRMB, CRML, and TITE-CRM methods, the parameters for the EWOC function are defined in the “Input parameters” tab. However, running this method using the bcrm function (package bcrm) [20] requires additional variables. This function allows the prior distribution of the model parameter to be selected as

either a gamma, uniform, or lognormal distribution along with its respective associated parameters. Two supplementary parameters are added according to the EWOC method. The “pointest” parameter, between 0 and 0.5, provides overdose control and represents the maximum acceptable proportion of overdosed patients. A parameter “Dose skipping constraint” was also added to avoid skipping doses in escalation immediately after a toxic outcome. The “Include” and “Results” tabs are similar to those for the previous methods. Supplementary Fig. 1b shows the view obtained for the EWOC method.
TITE-EWOC

The TITE-EWOC and EWOC views are quite similar. However, for TITE-CRM, in addition to basic input parameters, the “Observation window” is required to record the follow-up. The user can also enter an initial value for the parameter model with the field “Prior alpha value”. Given that the TITE-EWOC method accounts for the follow-up of each patient and some updated follow-ups are expected at the inclusion of new patients, we proposed some new functionalities for the “Include” tab. A “Modif patient data” button was added to update, if needed, the data of the previous patient included. The possibility of entering multiple patient data at the same time is now available by clicking on “Yes” in the pop-up display after selecting “New patient”. The data to be entered are “dose level,” “DLT,” and “Follow-up” of each new patient as comma-separated values. Note that it is always possible to enter patient data sequentially by clicking “No” in the pop-up display. Supplementary

Fig. 1c shows the view obtained for the TITE-EWOC method.

Simulating a study

We simulated 1000 trials of $n = 25$ patients according to the parameters used in the motivating example (Fig. 3). In addition, the probabilities of toxicity that would be observed in a real-life setting for each dose level were assumed to be the following: 0.003, 0.016, 0.047, 0.107, 0.196, and 0.305. In practice, it is usually recommended that several scenarios be considered. Dose level 5 (0.196) is the closest to the MTD (0.2).

The simulations revealed that the expected dose (dose level 3) is recommended (Fig. 3c, column “MTD dist.”) as the MTD (target DLT rate = 20%) in 3.3% of the experiments. Similarly, the dose levels 4, 5, and 6 are recommended in 26.4, 45.5, and 24.6% of the experiments, respectively. Dose levels 1 and 2 were shown to be poor or not recommended. In accordance with the true probabilities of toxicity and despite the wrong guess for the

<figure><img src="image_8.png" /></figure>
prior MTD, dose level 5 was the most recommended (45.5%). The number of patients treated at each dose level (Fig. 3c, column “Av. patient”) was distributed as follows: 0.373 patients (1.49%) at dose level 1, and 0.604 (2.42%), 3.114 (12.45%), 6.102 (24.41%), 7.524 (30.1%), and 7.283 (29.13%) at dose levels 2, 3, 4, 5, and 6, respectively. Four out of ten patients were treated at doses lower than the MTD, three out of ten at the MTD, and three out of ten at dose level 6 (overdose). The average results might be improved with better choices of input parameters or with another dose-escalation design. For instance, the overdose rate can be decreased using an EWOC design.

## Discussion

Several model-based dose-finding methods exist in oncology drug development. Despite a gradual increase in their use, these methods remain underutilized in practice due to their complexity.

Unlike the software already available for implementing dose-escalation designs for a phase I clinical trial, we developed an interface, GUIP1, available on GitHub (https://github.com/ddinart/GUIP1), that facilitates the use of these adaptive mathematical methods based on the modeling of the dose-toxicity relationship. GUIP1 is implemented using the free software R, which is widely used by statisticians in oncology. GUIP1 simplifies the use of these dose escalation methods and is designed to be fairly simple for beginners in R. Furthermore, its user-friendly interface is an easy-to-use solution that offers multiple possibilities.

GUIP1 provides the possibility of managing real clinical trials using file management options with automatic backup of study and/or simulation results. It uses libraries already published, tested, and validated by the scientific community and allows access to the data at all times for further analysis if necessary. In addition to being an interface for running statistical dose-escalation models, GUIP1 allows the real-time management of phase I clinical trials with an automatic backup of the results and the ability to export .xlx files. The export module returns the results of a study or simulations. This service provides full traceability of the study and is very convenient. The tools used to develop GUIP1 are very flexible, and this makes it possible to update the software at any time by adding new methods not yet implemented under R or by adding stopping criteria recently published in the literature. The stopping rules implemented in our interface do not extend beyond the existing rules in the package bcrm. Other stopping rules, including those based on binary trees, have been suggested [24]; however, they have not been made available in the interface. Therefore, adding more existing rules would be an interesting enhancement.

## Conclusions

One challenge with the conception of GUIP1 is a lack of information on the use of the tcltck package. Hence, developing the attractiveness and appeal of the interface as well as adding new functionality might be a future goal. We hope that this interface will help spread the adoption of new methods and increase their use in future therapeutic trials.

## Availability and requirements

**Project name:** GUIP1 Package

**Project home page:** https://github.com/ddinart/GUIP1

**Operating system(s):** Platform independent

**Programming language:** R

**Other requirements:** R v3.5.2

**License:** GPL (>= 2)

**Any restrictions to use by non-academics:** None

## Supplementary information

**Supplementary information** accompanies this paper at https://doi.org/10.1186/s12911-020-01149-3.

**Additional file 1:** Figure 1. Implementation of GUIP1 (a) CRML, (b) TITE-CRM, (c) EWOC, and (d) TITE-EWOC (motivating example).

## Abbreviations

**BOIN:** Bayesian optimal interval; CRM: Continual reassessment method;

**CRMB:** Continual reassessment method with Bayesian estimation;

**CRML:** Continual reassessment method with maximum likelihood estimation;

**DLT:** Dose-limiting toxicity; EWOC: Escalation with overdose control;

**FIH:** First-in-human; INCa: French National Cancer Institute; MTD: Maximum

tolerated dose; mTPI: Modified toxicity probability interval;

**RP2D:** Recommended phase II dose; TITE-CRM: Time-to-event-CRM; TITE-

EWOC: Time-to-event EWOC

## Acknowledgements

The authors wish to thank Dr. Ravi Nookala of Institut Bergonié for medical writing services.

## Authors' contributions

DD and JF developed the interface. CB and CM supervised the project. DD, JF, CT, CB and CM wrote the manuscript. DD, JF, DT, AM, CT, SG, MCLD, CB and CM discussed the results and contributed to the final manuscript. All authors have read and approved the manuscript.

## Funding

The presented research was supported by a grant of the French National Cancer Institute (INCA-Grant n°SHS-ESP 2015–164). Note: The funding body had no other role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.)

## Availability of data and materials

The interface GUIP1 is available on GitHub (https://github.com/ddinart/GUIP1).

## Ethics approval and consent to participate

Not applicable.

## Consent for publication

Not applicable.

## Competing interests

None of the authors have any competing interests.
## Author details

¹Inserm CIC1401, Module Epidémiologie Clinique, Institut Bergonié, Bordeaux, France. ²Institut du Cancer Montpellier (ICM), Université de Montpellier, Montpellier, France. ³Memorial Sloan Kettering Cancer Center, New York, USA. ⁴Centre Oscar Lambret, Lille, France. ⁵Université Paris-Saclay, Université Paris-Sud, UVSQ, CESP, INSERM, Villejuif, France. ⁶Institut de Recherche en Cancérologie de Montpellier INSERM U1194, Université de Montpellier, Montpellier, France.

Received: 18 October 2019 Accepted: 5 June 2020

Published online: 24 June 2020

## References

1. Storer BE. Design and analysis of phase I clinical trials. Biometrics. 1989;45(3): 925-37.

2. Le Tourneau C, Gan HK, Razak ARA, Paoletti X. Efficiency of New Dose Escalation Designs in Dose-Finding Phase I Trials of Molecularly Targeted Agents. PLoS One. 2012;7(12):e51039. Stemmer SM, editor.

3. Reiner E, Paoletti X, O'Quigley J. Operating characteristics of the standard phase I clinical trial design. Comput Stat Data Anal. 1999;30(3):303-15.

4. Rogatko A, Schoeneck D, Jonas W, Tighiouart M, Khuri FR, Porter A. Translation of innovative designs into phase I trials. J Clin Oncol. 2007; 25(31):4982-6.

5. O'Quigley J, Pepe M, Fisher L. Continual reassessment method: a practical design for phase 1 clinical trials in cancer. Biometrics. 1990;46(1):33-48.

6. O'Quigley J, Shen LZ. Continual reassessment method: a likelihood approach. Biometrics. 1996;52(2):673-84.

7. Cheung YK, Chappell R. Sequential designs for phase I clinical trials with late-onset toxicities. Biometrics. 2000;56(4):1177-82.

8. Babb J, Rogatko A, Zacks S. Cancer phase I clinical trials: efficient dose escalation with overdose control. Stat Med. 1998;17(10):1103-20.

9. Yin G, Yuan Y. Bayesian model averaging continual reassessment method in phase I clinical trials. J Am Stat Assoc. 2009;104(487):954-68.

10. Liu S, Yin G, Yuan Y. Bayesian data augmentation dose finding with continual reassessment method and delayed toxicity. Ann Appl Stat. 2013; 7(4):1837-2457.

11. Van Meter EM, Garrett-Mayer E, Bandyopadhyay D. Proportional odds model for dose-finding clinical trial designs with ordinal toxicity grading. Stat Med. 2011;30(17):2070-80.

12. Van Meter EM, Garrett-Mayer E, Bandyopadhyay D. Dose-finding clinical trial design for ordinal toxicity grades using the continuation ratio model: an extension of the continual reassessment method. Clin Trials. 2012;9(3):303-13.

13. Braun TM. The bivariate continual reassessment method. Extending the CRM to phase I trials of two competing outcomes. Control Clin Trials. 2002; 23(3):240-56.

14. Yan F, Mandrekar SJ, Yuan Y. Keyboard: a novel Bayesian toxicity probability interval Design for Phase I Clinical Trials. Clin Cancer Res. 2017;23(15):3994-4003.

15. Zhou H, Murray TA, Pan H, Yuan Y. Comparative review of novel model-assisted designs for phase I clinical trials. Stat Med. 2018;37(14):2208-22.

16. Yuan Y, Lee JJ, Hilsenbeck SG. Model-assisted designs for early-phase clinical trials: simplicity meets superiority. JCO Precis Oncol. 2019;3:1-12.

17. Ji Y, Liu P, Li Y, Bekele BN. A modified toxicity probability interval method for dose-finding trials. Clin Trials. 2010;7(6):653-63.

18. Yuan Y, Hess KR, Hilsenbeck SG, Gilbert MR. Bayesian optimal interval design: a simple and well-performing design for phase I oncology trials. Clin Cancer Res. 2016;22(17):4291-301.

19. Kramar A, Houédé N, Paoletti X. np1: a computer program for dose escalation strategies in phase I clinical trials. Comput Methods Prog Biomed. 2007;88(1):8-17.

20. Sweeting M, Mander A, Sabin T. bcrm: Bayesian Continual Reassessment Method Designs for Phase I Dose-Finding Trials. J Stat Softw. 2013;54(13) [cited 2017 Jun 7]. Available from: http://www.jstatsoft.org/v54/i13/.

21. Mauguen A, Le Deley MC, Zohar S. Dose-finding approach for dose escalation with overdose control considering incomplete observations. Stat Med. 2011;30(13):1584-94.

22. Paoletti X, Kramar A. A comparison of model choices for the continual reassessment method in phase I cancer trials. Stat Med. 2009;28(24):3012-28.

23. Lee SM, Ying Kuen Cheung null. Model calibration in the continual reassessment method. Clin Trials. 2009;6(3):227-38.

24. O'Quigley J. Continual reassessment designs with early termination. Biostatistics. 2002;3(1):87-99.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions