gammatable | R Documentation |
\gamma_L
and \gamma_R
in the calibration-free odds (CFO) designGenerate all the possible thresholds under different m_C
, m_L
and m_R
gammatable(npatient, target,
para.prior = list(alp.prior = target, bet.prior = 1 - target))
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 |
The gammatable()
function returns a list object comprising the following elements:
gammatb.left: the table of threshold \gamma_L
under different m_L
and m_C
where m_C
and m_L
represent the number of patients at current dose level and left dose level.
gammatb.right: the table of threshold \gamma_R
under different m_R
and m_C
where m_C
and m_R
represent the number of patients at current dose level and right dose level.
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]
Jialu Fang, Ninghao Zhang, Wenliang Wang, and Guosheng Yin
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.
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.