gammatable {CFO} | R Documentation |
Generating table of threshold \gamma_L
and \gamma_R
in the calibration-free odds (CFO) design
Description
Generate all the possible thresholds under different m_C
, m_L
and m_R
Usage
gammatable(npatient, target,
para.prior = list(alp.prior = target, bet.prior = 1 - target))
Arguments
npatient |
the numbers of patients involved in the trial. |
target |
the target DLT rate. |
para.prior |
the prior parameters for a beta distribution, where set as |
Value
The gammatable()
function returns a list object comprising the following elements:
gammatb.left: the table of threshold
\gamma_L
under differentm_L
andm_C
wherem_C
andm_L
represent the number of patients at current dose level and left dose level.gammatb.right: the table of threshold
\gamma_R
under differentm_R
andm_C
wherem_C
andm_R
represent the number of patients at current dose level and right dose level.
Note
This function generate two matrices. gammatb.left
contains the threshold \gamma_L
,
and gammatb.right
contains the threhold \gamma_R
. For matrix gammatb.left
, the row index represent the number of patients
at left dose level, and the column index represent the number of patients at current dose level. For matrix gammatb.right
, the row index represent the number of patients
at right dose level, and the column index represent the number of patients at current dose level.
For example, if you want to get the threshold \gamma_L
in the case of m_C = 12, m_L = 13
, you can reach it by result$gammatb.left[13,12]
Author(s)
Jialu Fang, Ninghao Zhang, Wenliang Wang, and Guosheng Yin
References
Jin H, Yin G (2022). CFO: Calibration-free odds design for phase I/II clinical trials. Statistical Methods in Medical Research, 31(6), 1051-1066.
Examples
npatient <- 3; target <- 0.3
para.prior = list(alp.prior = target, bet.prior = 1 - target)
result <- gammatable(npatient, target, para.prior)
plot(result)
#This example may cost you a long time to run
npatient <- 30; target <- 0.3
para.prior = list(alp.prior = target, bet.prior = 1 - target)
result <- gammatable(npatient, target, para.prior)
plot(result)