dabcd_max_power can be used for doubly adaptive biased coin design with maximal power
strategy for binary outcomes, targeting generalized Neyman allocation and generalized RSIHR allocation. The return
of this function is a vector of allocation probabilities to each arm, with the pre-specified number of participants in the trial.
dabcd_max_power(NN, Ntotal1, armn, BB, type, dabcd = FALSE, gamma = 2)A vector of allocation probabilities to each arm.
a vector representing the number of participants with success results for each arm estimated from the current data.
a vector representing the total number of participants for each arm estimated from the current data.
number of total arms in the trial.
the minimal allocation probability for each arm, which is within the range of \([0,1/armn]\).
allocation type, with choices from 'RSIHR' and 'Neyman'.
an indicator of whether to apply Hu & Zhang's formula (Hu2004RARtrials), with choices from FALSE and TRUE. TRUE represents allocation probabilities calculated using Hu & Zhang's formula; FALSE represents allocation probabilities calculated before applying Hu & Zhang's formula. Default value is set to FALSE.
tuning parameter in Hu & Zhang's formula (Hu2004RARtrials). When dabcd=FALSE, this parameter does not need
to be specified. Default value is set to 2.
Chuyao Xu, Thomas Lumley, Alain Vandal
The function simulates allocation probabilities for doubly adaptive biased coin design with maximal power strategy targeting
generalized Neyman allocation with 2-5 arms which is provided in Tymofyeyev2007RARtrials or
generalized RSIHR allocation with 2-3 arms which is provided in Jeon2010RARtrials, with modifications for typos
in Sabo2016RARtrials. All of those methods are not smoothed. The output of this function is based on Hu \& Zhang's formula Hu2004RARtrials.
With more than two armd the one-sided nominal level of each test is alphaa divided by arm*(arm-1)/2; a Bonferroni correction.
Hu2004RARtrials
Tymofyeyev2007RARtrials
Jeon2010RARtrials
Sabo2016RARtrials
dabcd_max_power(NN=c(54,67,85,63,70),Ntotal1=c(100,88,90,94,102),armn=5,BB=0.2, type='Neyman')
dabcd_max_power(NN=c(54,67,85,63),Ntotal1=c(100,88,90,94),armn=4,BB=0.2, type='Neyman')
Run the code above in your browser using DataLab