Learn R Programming

RARtrials (version 0.0.2)

dabcd_max_power: Allocation Probabilities Using Doubly Adaptive Biased Coin Design with Maximal Power Strategy for Binary Endpoint

Description

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.

Usage

dabcd_max_power(NN, Ntotal1, armn, BB, type, dabcd = FALSE, gamma = 2)

Value

A vector of allocation probabilities to each arm.

Arguments

NN

a vector representing the number of participants with success results for each arm estimated from the current data.

Ntotal1

a vector representing the total number of participants for each arm estimated from the current data.

armn

number of total arms in the trial.

BB

the minimal allocation probability for each arm, which is within the range of \([0,1/armn]\).

type

allocation type, with choices from 'RSIHR' and 'Neyman'.

dabcd

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.

gamma

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.

Author

Chuyao Xu, Thomas Lumley, Alain Vandal

Details

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.

References

Hu2004RARtrials

Tymofyeyev2007RARtrials

Jeon2010RARtrials

Sabo2016RARtrials

Examples

Run this code
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