pipe.design: Dual-Agent Dose Escalation for Phase I Trials using the PIPE...
In pipe.design: Dual-Agent Dose Escalation for Phase I Trials using the PIPE Design
Description
Usage
Arguments
Details
Value
Note
Author(s)
References
See Also
Examples
View source: R/pipe.design_0.5.1.R
Description
Implements the Product of Independent beta Probabilities dose Escalation (PIPE) design for dual-agent Phase I trials as described in Mander A.P. and Sweeting M.J (2015) Statistics in Medicine. <doi: 10.1002/sim.6434>
Usage
1
2
3
4
5
Arguments
N
Sample size for the trial. If not specified this is assumed to be one greater than the sample size of the patients recruited thus far (i.e. the number of rows of data plus one)
S
Number of simulations to conduct. If set to 1 then a single trial is run and invokes different output from print.pipe and plot.pipe
c
Cohort size. Must be divisible by N
theta
The target toxicity probability
pi
A matrix with rows denoting levels of drug A and columns denoting levels of drug B. Each element gives the true probability of the outcome (dose-limiting toxicity) for that dose combination. If omitted then the trial will recommend only for the next dose combination
prior.med
A matrix with rows denoting levels of drug A and columns denoting levels of drug B. Each element gives the prior median probability of the outcome (dose-limiting toxicity) for that dose combination. Should be specified in combination with prior.ss
prior.ss
A matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the prior sample size for that dose combination. Should be specified in combination with prior.med
strategy
A string detailing the dose escalation strategy. Options are "ss" (sample size) and "ss-random" (weighted randomisation by sample size). See details below
admis
A string detailing how to choose admissible doses around the maximum tolerated contour (MTC). Options are "adjacent" and "closest". See details below
constraint
A string detailing any dose-skipping constraint to be applied. Options are "neighbouring" or "no.dose.skip" (see details below). In addition, to prevent diagonal escalation the option "-nodiag" can be appended, i.e. "neighbouring-nodiag" or "no.dose.skip-nodiag". Defaults to no constraint (constraint="none").
epsilon
Any number between 0 and 1 giving a safety threshold. Doses whose (weighted) posterior probability(dose>MTC)>epsilon will not be experimented on. See details below
mode
A string indicating the running mode of the program. Options are "sim" (the default, used for simulations), "nodlt" (every patient recruited is assumed to be DLT free), or "alldlt" (every patient recruited is assumed to suffer a DLT).
data
An optional named data frame giving information about dose and toxicity from previously recruited patients. If missing then it is assumed that no data have thus far been collected. Contains the following variables:
patient Recruited patient numbers, 1,...,n
doseA Dose levels of recruited patients for drug A
doseB Dose levels of recruited patients for drug B
tox An indicator variable for each patient (1=toxicity, 0=no toxicity)
a
(Optional). A matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the shape1 parameter (a) from a prior Beta(a,b) distribution. Can be used as an alternative to specifying prior.med and prior.ss and must be used in conjunction with b
b
(Optional). A matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the shape2 parameter (b) from a prior Beta(a,b) distribution. Can be used as an alternative to specifying prior.med and prior.ss and must be used in conjunction with a
alternate
Logical variable (defaults to FALSE). Should the design always deescalate if above the MTC and escalate if below (subject to admissibility, safety and dose skipping constraints)?
uppertox.constraint
Any number between 0 and 1 giving the upper toxicity safety constraint. No dose combination that lies above the most likely contour for this constraint may be dosed. Defaults to NULL
stop
(Optional). Value of a stopping threshold. An alternative stopping rule to a fixed sample size specification (through N) or termination of a trial due to safety constraints (epsilon or uppertox.constraint). The trial is stopped if the posterior probability of being > theta at the lowest dose combination is greater than stop
non.admissible
(Optional) A matrix with rows denoting levels of drug A and columns denoting levels of drug B. Each element is logical indicating whether the dose combination is non-admissible (TRUE) or not. Non-admissible doses can never be experimented on in during the trial.
seed
(Optional) Starting seed to allow replication of the trial or simulation.
Details
This function implements a Product of Independent Beta Probabilities (PIPE) Phase I dose-escalation design for a dual-agent trial. An independent Beta prior is placed on the probability of dose limiting toxicity (DLT) for each dose combination. Simple conjugate Bayesian analysis is used to obtain the posterior at each dose combination. This model is considered a 'working model' only as it does not account for monotonically increasing toxicity by dose.
The dose combination assigned to the next cohort of patients is, however, obtained by calculating a maximum tolerated contour (MTC) that bisects the dose space and respects monotonicity. Admissible dose combinations considered for the next cohort are those that are either ‘adjacent’ or ‘closest’ to the MTC (see Mander A.P. and Sweeting M.J., Statistics in Medicine 2015) and can be specified using the admis argument. If there is more than one dose combination then either the one with the smallest sample size (least experimented on) is chosen, by specifying strategy="ss", or admissible doses are chosen at random, weighted by the inverse of the sample size (strategy="ss-random").
Admissible doses can be constrained to be within a set of neighbouring doses to the current dose combination (constraint="neighbouring") or within the set of doses consisting of up to one dose level above any previously experimented drug A and drug B level (constraint="no.dose.skip"). To additionally prevent diagonal escalation (e.g. where both drugs are escalated together), use either (constraint="neighbouring-nodiag") or (constraint="no.dose.skip-nodiag").
In addition, a safety constraint can be added, to avoid potential overdosing. Dose combinations whose expected probability of being above the MTC is greater than a threshold epsilon are excluded from the admissible doses for the next cohort. The expected probability is calculated by averaging over all possible MTCs, weighted by their posterior probabilities. An alternative safety constraint can be used by specifying uppertox.constraint=c, where c is the upper toxicity level from which a contour is derived (e.g. c=0.6 calculates the most likely contour for a 60% chance of DLT). No doses above this contour will be used.
The program will either run as a simulation (the default), or by using a deterministic set of outcomes (all DLTs / no DLTs) specified using mode="alldlt" and mode="nodlt", respectively.
If the user specifies a data-frame using the data argument, then every simulated trial will begin with the specified doses and DLT outcomes. This option is useful for assess operating characteristics mid-way through the trial, or just to obtain the dose for the next cohort.
A trial of size N is generally run, except if no doses are deemed safe according to one of the two safety constraints described above, at which point the trial will stop. Users may additionally specify a stopping rule using stop=s, where the trial is stopped if the posterior probability that the lowest dose combination is above the target toxicity probability is greater than s. Note, this stopping rule relies only on the prior and data collected for the first dose combination and does not borrow strength from other doses. Hence it should be used with caution.
One can force deescalation if above the MTC and escalation if below (subject to meeting other specified constraints) by specifying alternate=TRUE. This constraint ensures the design follows the properties of 'coherence' (Cheung YK. Biometrika 2005).
If a single trial is specified by S=1 then by additionally specifying densities=TRUE the returned object will include calculations of the probability of toxicity densities for each dose-combination after each cohort (with cohort 0 representing the prior). The densities are calculated over a grid of probabilities α_1,...,α_K through weighting by the contour probabilities for each α_k:
p(α_k< π_{i,j} ≤ α_{k+1}) = ∑_{C_s \in \mathcal{C}} C_s[i,j] P_{C_s}(α_{k+1})- ∑_{C_s \in \mathcal{C}} C_s[i,j] P_{C_s}(α_{k})
Further details regarding the notation can be found in Mander A.P. and Sweeting M.J. (2015) Statistics in Medicine. <doi: 10.1002/sim.6434>
Value
pipe.design returns an object of class "pipe" or "pipe.sim"; the latter occuring when more than one simulation has been conducted.
The function print (i.e. print.pipe or print.pipe.sim) can be used to obtain summary information about the design used, the data observed, the maximum tolerated contour estimate, and the next recommended dose level.
An object of class "pipe" (when N=1) or "pipe.sim" (when N>1) is a list with the following components:
r.simA list of matrices showing the number of DLTs for each dose combination in each simulation
n.simA list of matrices showing the number of patients dosed at each combination in each simulation
rec.i.simAn S x (N/c) matrix with elements representing the level of drug A for each cohort
rec.j.simAn S x (N/c) matrix with elements representing the level of drug B for each cohort
expA matrix showing the experimentation proportions for each dose combination averaged across simulations
recA matrix showing the proportion of times each dose combination is recommended, where the denominator is the total number of recommended phase II doses over the simulations
dltsA vector giving the average number of DLTs occuring in each simulated trial
mat.listA list of binary matrices showing the most likely MTC before each cohort is recruited. Output only if S=1
uppermat.listA list of binary matrices showing the most likely contour for uppertox.constraint before each cohort is recruited. Output only if S=1 and uppertox.constraint is used
uppermat2.listA list of binary matrices showing the doses that fail the constraint imposed by epsilon. Output only if S=1 and epsilon is used
r.listA list of matrices showing the cumulative number of DLTs for each dose combination before each cohort is recruited. Output only if S=1
n.listA list of matrices showing the cumulative number of patients at each dose combination before each cohort is recruited. Output only if S=1
n.rpIIA vector listing the number of doses recommended for Phase II from each simulated trial
no.not.treatedThe total number of patients not treated (in the trial) accumulated across the simulated trials due to a trial stopping early
piThe true probability of DLT for each dose combination, as specified in the input argument
thetaThe target toxicity probability, as specified in the input argument
rpII.listA list giving the recommended Phase II dose combinations for each simulation
aA matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the shape1 parameter (a) from a prior Beta(a,b) distribution
bA matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the shape2 parameter (b) from a prior Beta(a,b) distribution
plot.densityA list of objects of length c+1. The kth element is a ggplot object giving the toxicity density plots for each dose-combination after the (k-1)st cohort has been recruited
Note
The use of the uppertox.constraint safety constraint is not documented in Mander and Sweeting, Statistics in Medicine 2015, and some preliminary investigations have shown that this constraint can cause rigidity in the escalation process, and therefore should be used with caution.
Author(s)
Michael Sweeting mjs212@medschl.cam.ac.uk (University of Cambridge, UK)
References
Mander A.P., Sweeting M.J. A product of independent beta probabilities dose escalation design for dual-agent phase I trials. Statistics in Medicine (2015) 34(8): 1261–1276. <doi: 10.1002/sim.6434>
Cheung Y.K. Coherence principles in dose-finding studies. Biometrika (2005) 92(4): 863–873.
See Also
print.pipe, print.pipe.sim, plot.pipe, plot.pipe.sim
Examples
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
## Reproducing Figure 5 from Mander and Sweeting, Statistics in Medicine 2015.
## The true probability of DLT for a 6x6 grid of dose combinations
## (Scenario 3 from Mander and Sweeting, Statistics in Medicine 2015)
pi <-matrix(c(0.02,0.10,0.20,0.30,0.35,0.45,0.06,0.14,0.24,0.34,0.39,0.49,0.12,0.20,
0.30,0.40,0.45,0.55,0.17,0.25,0.35,0.45,0.50,0.60,0.22,0.30,0.40,0.50,0.60,0.70,0.30,
0.38,0.48,0.58,0.68,0.78),nrow=6,ncol=6)
## Using a weak prior with prior sample size 1/36 for each dose combination and prior
## median of Scenario 1
prior.med<-matrix(c(0.02,0.03,0.06,0.10,0.18,0.23,0.03,0.05,0.09,0.13,0.21,0.27,0.06,
0.09,0.14,0.18,0.26,0.30,0.11,0.14,0.18,0.23,0.30,0.36,0.18,0.21,0.26,0.30,0.39,0.44,
0.23,0.27,0.3,0.36,0.44,0.49),nrow=6,ncol=6)
## Using a neighbouring escalation constraint
## Selecting the closest admissible doses
## Using a safety constraint with epsilon = 0.8
## Obtain dose recommendations for first cohort
cohort1<-pipe.design(N=2,S=1,c=2,theta=0.3,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
cohort1 ## Recommends starting at (1,1)
## Second cohort
## Assume no toxicities are seen in first cohort
data1<-data.frame(patient=1:2,doseA=rep(1,2),doseB=rep(1,2),tox=rep(0,2))
cohort2<-pipe.design(data=data1,S=1,c=2,theta=0.3,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
cohort2 ## Recommends dosing at (2,2)
## Third cohort
## Assume one toxicity is seen in second cohort
data2<-rbind(data1,data.frame(patient=3:4,doseA=rep(2,2),doseB=rep(2,2),tox=c(1,0)))
cohort3<-pipe.design(data=data2,S=1,c=2,theta=0.3,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
cohort3 ## Recommends dosing at (1,3)
## Simulating a single trial with sample size of 40, cohort size 2
## and target toxicity of 0.3, using true probabilities of toxicity pi
## Not run:
set.seed(700)
example<-pipe.design(N=40,S=1,c=2,theta=0.3,pi=pi,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
example
plot(example)
## End(Not run)
## Not run:
## Run a Shiny GUI version of pipe.design
runShinyPIPE()
## End(Not run)
## Not run:
## Conducting a simulation study (1000 trials) to investigate
## operating characteristics for Scenario 3
## (as reported in Table III in Mander and Sweeting, Statistics in Medicine 2015)
set.seed(262)
scen3<-pipe.design(N=40,S=1000,c=2,theta=0.3,pi=pi,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
scen3
## End(Not run)
pipe.design documentation built on May 2, 2019, 12:39 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Arguments Details Value Note Author(s) References See Also Examples
View source: R/pipe.design_0.5.1.R
Description
Implements the Product of Independent beta Probabilities dose Escalation (PIPE) design for dual-agent Phase I trials as described in Mander A.P. and Sweeting M.J (2015) Statistics in Medicine. <doi: 10.1002/sim.6434>
Usage
1 2 3 4 5 |
Arguments
N |
Sample size for the trial. If not specified this is assumed to be one greater than the sample size of the patients recruited thus far (i.e. the number of rows of |
S |
Number of simulations to conduct. If set to 1 then a single trial is run and invokes different output from |
c |
Cohort size. Must be divisible by |
theta |
The target toxicity probability |
pi |
A matrix with rows denoting levels of drug A and columns denoting levels of drug B. Each element gives the true probability of the outcome (dose-limiting toxicity) for that dose combination. If omitted then the trial will recommend only for the next dose combination |
prior.med |
A matrix with rows denoting levels of drug A and columns denoting levels of drug B. Each element gives the prior median probability of the outcome (dose-limiting toxicity) for that dose combination. Should be specified in combination with |
prior.ss |
A matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the prior sample size for that dose combination. Should be specified in combination with |
strategy |
A string detailing the dose escalation strategy. Options are "ss" (sample size) and "ss-random" (weighted randomisation by sample size). See details below |
admis |
A string detailing how to choose admissible doses around the maximum tolerated contour (MTC). Options are "adjacent" and "closest". See details below |
constraint |
A string detailing any dose-skipping constraint to be applied. Options are "neighbouring" or "no.dose.skip" (see details below). In addition, to prevent diagonal escalation the option "-nodiag" can be appended, i.e. "neighbouring-nodiag" or "no.dose.skip-nodiag". Defaults to no constraint (constraint="none"). |
epsilon |
Any number between 0 and 1 giving a safety threshold. Doses whose (weighted) posterior probability(dose>MTC)>epsilon will not be experimented on. See details below |
mode |
A string indicating the running mode of the program. Options are "sim" (the default, used for simulations), "nodlt" (every patient recruited is assumed to be DLT free), or "alldlt" (every patient recruited is assumed to suffer a DLT). |
data |
An optional named data frame giving information about dose and toxicity from previously recruited patients. If missing then it is assumed that no data have thus far been collected. Contains the following variables:
|
a |
(Optional). A matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the |
b |
(Optional). A matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the |
alternate |
Logical variable (defaults to FALSE). Should the design always deescalate if above the MTC and escalate if below (subject to admissibility, safety and dose skipping constraints)? |
uppertox.constraint |
Any number between 0 and 1 giving the upper toxicity safety constraint. No dose combination that lies above the most likely contour for this constraint may be dosed. Defaults to NULL |
stop |
(Optional). Value of a stopping threshold. An alternative stopping rule to a fixed sample size specification (through |
non.admissible |
(Optional) A matrix with rows denoting levels of drug A and columns denoting levels of drug B. Each element is logical indicating whether the dose combination is non-admissible (TRUE) or not. Non-admissible doses can never be experimented on in during the trial. |
seed |
(Optional) Starting seed to allow replication of the trial or simulation. |
Details
This function implements a Product of Independent Beta Probabilities (PIPE) Phase I dose-escalation design for a dual-agent trial. An independent Beta prior is placed on the probability of dose limiting toxicity (DLT) for each dose combination. Simple conjugate Bayesian analysis is used to obtain the posterior at each dose combination. This model is considered a 'working model' only as it does not account for monotonically increasing toxicity by dose.
The dose combination assigned to the next cohort of patients is, however, obtained by calculating a maximum tolerated contour (MTC) that bisects the dose space and respects monotonicity. Admissible dose combinations considered for the next cohort are those that are either ‘adjacent’ or ‘closest’ to the MTC (see Mander A.P. and Sweeting M.J., Statistics in Medicine 2015) and can be specified using the admis argument. If there is more than one dose combination then either the one with the smallest sample size (least experimented on) is chosen, by specifying strategy="ss", or admissible doses are chosen at random, weighted by the inverse of the sample size (strategy="ss-random").
Admissible doses can be constrained to be within a set of neighbouring doses to the current dose combination (constraint="neighbouring") or within the set of doses consisting of up to one dose level above any previously experimented drug A and drug B level (constraint="no.dose.skip"). To additionally prevent diagonal escalation (e.g. where both drugs are escalated together), use either (constraint="neighbouring-nodiag") or (constraint="no.dose.skip-nodiag").
In addition, a safety constraint can be added, to avoid potential overdosing. Dose combinations whose expected probability of being above the MTC is greater than a threshold epsilon are excluded from the admissible doses for the next cohort. The expected probability is calculated by averaging over all possible MTCs, weighted by their posterior probabilities. An alternative safety constraint can be used by specifying uppertox.constraint=c, where c is the upper toxicity level from which a contour is derived (e.g. c=0.6 calculates the most likely contour for a 60% chance of DLT). No doses above this contour will be used.
The program will either run as a simulation (the default), or by using a deterministic set of outcomes (all DLTs / no DLTs) specified using mode="alldlt" and mode="nodlt", respectively.
If the user specifies a data-frame using the data argument, then every simulated trial will begin with the specified doses and DLT outcomes. This option is useful for assess operating characteristics mid-way through the trial, or just to obtain the dose for the next cohort.
A trial of size N is generally run, except if no doses are deemed safe according to one of the two safety constraints described above, at which point the trial will stop. Users may additionally specify a stopping rule using stop=s, where the trial is stopped if the posterior probability that the lowest dose combination is above the target toxicity probability is greater than s. Note, this stopping rule relies only on the prior and data collected for the first dose combination and does not borrow strength from other doses. Hence it should be used with caution.
One can force deescalation if above the MTC and escalation if below (subject to meeting other specified constraints) by specifying alternate=TRUE. This constraint ensures the design follows the properties of 'coherence' (Cheung YK. Biometrika 2005).
If a single trial is specified by S=1 then by additionally specifying densities=TRUE the returned object will include calculations of the probability of toxicity densities for each dose-combination after each cohort (with cohort 0 representing the prior). The densities are calculated over a grid of probabilities α_1,...,α_K through weighting by the contour probabilities for each α_k:
p(α_k< π_{i,j} ≤ α_{k+1}) = ∑_{C_s \in \mathcal{C}} C_s[i,j] P_{C_s}(α_{k+1})- ∑_{C_s \in \mathcal{C}} C_s[i,j] P_{C_s}(α_{k})
Further details regarding the notation can be found in Mander A.P. and Sweeting M.J. (2015) Statistics in Medicine. <doi: 10.1002/sim.6434>
Value
pipe.design returns an object of class "pipe" or "pipe.sim"; the latter occuring when more than one simulation has been conducted.
The function print (i.e. print.pipe or print.pipe.sim) can be used to obtain summary information about the design used, the data observed, the maximum tolerated contour estimate, and the next recommended dose level.
An object of class "pipe" (when N=1) or "pipe.sim" (when N>1) is a list with the following components:
r.simA list of matrices showing the number of DLTs for each dose combination in each simulation
n.simA list of matrices showing the number of patients dosed at each combination in each simulation
rec.i.simAn S x (N/c) matrix with elements representing the level of drug A for each cohort
rec.j.simAn S x (N/c) matrix with elements representing the level of drug B for each cohort
expA matrix showing the experimentation proportions for each dose combination averaged across simulations
recA matrix showing the proportion of times each dose combination is recommended, where the denominator is the total number of recommended phase II doses over the simulations
dltsA vector giving the average number of DLTs occuring in each simulated trial
mat.listA list of binary matrices showing the most likely MTC before each cohort is recruited. Output only if
S=1uppermat.listA list of binary matrices showing the most likely contour for
uppertox.constraintbefore each cohort is recruited. Output only ifS=1anduppertox.constraintis useduppermat2.listA list of binary matrices showing the doses that fail the constraint imposed by
epsilon. Output only ifS=1andepsilonis usedr.listA list of matrices showing the cumulative number of DLTs for each dose combination before each cohort is recruited. Output only if
S=1n.listA list of matrices showing the cumulative number of patients at each dose combination before each cohort is recruited. Output only if
S=1n.rpIIA vector listing the number of doses recommended for Phase II from each simulated trial
no.not.treatedThe total number of patients not treated (in the trial) accumulated across the simulated trials due to a trial stopping early
piThe true probability of DLT for each dose combination, as specified in the input argument
thetaThe target toxicity probability, as specified in the input argument
rpII.listA list giving the recommended Phase II dose combinations for each simulation
aA matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the
shape1parameter (a) from a prior Beta(a,b) distributionbA matrix with rows denoting levels of drug A and columns levels of drug B. Each element gives the
shape2parameter (b) from a prior Beta(a,b) distributionplot.densityA list of objects of length
c+1. The kth element is a ggplot object giving the toxicity density plots for each dose-combination after the (k-1)st cohort has been recruited
Note
The use of the uppertox.constraint safety constraint is not documented in Mander and Sweeting, Statistics in Medicine 2015, and some preliminary investigations have shown that this constraint can cause rigidity in the escalation process, and therefore should be used with caution.
Author(s)
Michael Sweeting mjs212@medschl.cam.ac.uk (University of Cambridge, UK)
References
Mander A.P., Sweeting M.J. A product of independent beta probabilities dose escalation design for dual-agent phase I trials. Statistics in Medicine (2015) 34(8): 1261–1276. <doi: 10.1002/sim.6434>
Cheung Y.K. Coherence principles in dose-finding studies. Biometrika (2005) 92(4): 863–873.
See Also
print.pipe, print.pipe.sim, plot.pipe, plot.pipe.sim
Examples
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 | ## Reproducing Figure 5 from Mander and Sweeting, Statistics in Medicine 2015.
## The true probability of DLT for a 6x6 grid of dose combinations
## (Scenario 3 from Mander and Sweeting, Statistics in Medicine 2015)
pi <-matrix(c(0.02,0.10,0.20,0.30,0.35,0.45,0.06,0.14,0.24,0.34,0.39,0.49,0.12,0.20,
0.30,0.40,0.45,0.55,0.17,0.25,0.35,0.45,0.50,0.60,0.22,0.30,0.40,0.50,0.60,0.70,0.30,
0.38,0.48,0.58,0.68,0.78),nrow=6,ncol=6)
## Using a weak prior with prior sample size 1/36 for each dose combination and prior
## median of Scenario 1
prior.med<-matrix(c(0.02,0.03,0.06,0.10,0.18,0.23,0.03,0.05,0.09,0.13,0.21,0.27,0.06,
0.09,0.14,0.18,0.26,0.30,0.11,0.14,0.18,0.23,0.30,0.36,0.18,0.21,0.26,0.30,0.39,0.44,
0.23,0.27,0.3,0.36,0.44,0.49),nrow=6,ncol=6)
## Using a neighbouring escalation constraint
## Selecting the closest admissible doses
## Using a safety constraint with epsilon = 0.8
## Obtain dose recommendations for first cohort
cohort1<-pipe.design(N=2,S=1,c=2,theta=0.3,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
cohort1 ## Recommends starting at (1,1)
## Second cohort
## Assume no toxicities are seen in first cohort
data1<-data.frame(patient=1:2,doseA=rep(1,2),doseB=rep(1,2),tox=rep(0,2))
cohort2<-pipe.design(data=data1,S=1,c=2,theta=0.3,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
cohort2 ## Recommends dosing at (2,2)
## Third cohort
## Assume one toxicity is seen in second cohort
data2<-rbind(data1,data.frame(patient=3:4,doseA=rep(2,2),doseB=rep(2,2),tox=c(1,0)))
cohort3<-pipe.design(data=data2,S=1,c=2,theta=0.3,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
cohort3 ## Recommends dosing at (1,3)
## Simulating a single trial with sample size of 40, cohort size 2
## and target toxicity of 0.3, using true probabilities of toxicity pi
## Not run:
set.seed(700)
example<-pipe.design(N=40,S=1,c=2,theta=0.3,pi=pi,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
example
plot(example)
## End(Not run)
## Not run:
## Run a Shiny GUI version of pipe.design
runShinyPIPE()
## End(Not run)
## Not run:
## Conducting a simulation study (1000 trials) to investigate
## operating characteristics for Scenario 3
## (as reported in Table III in Mander and Sweeting, Statistics in Medicine 2015)
set.seed(262)
scen3<-pipe.design(N=40,S=1000,c=2,theta=0.3,pi=pi,prior.med=prior.med,
prior.ss=matrix(1/36,ncol=6,nrow=6),strategy="ss",constraint="neighbouring",
epsilon=0.8,admis="closest",alternate=FALSE)
scen3
## End(Not run)
|
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.