MABOUST {MABOUST}

Conduct the MABOUST Trial design.

Description

Performs posterior sampling for the MABOUST design and determines whether the trial should continue and what treatment(s) are optimal.

Usage

MABOUST(
  Y,
  T1,
  X,
  ACTIVE,
  FUTILITY,
  nTreat,
  nCat,
  UT,
  DeltaVEC,
  gamma,
  PSPIKE,
  ADJ,
  B
)

Arguments

Y

Ordinal Outcome Vector, labeled 1,...,J

T1

Treatment Indicator, labeled 1,...,K.

X

Matrix of patient covariates.

ACTIVE

Binary indicator of active treatments. This vector must be length K, and have a 1 for each entry corresponding to an active treatment and 0 otherwise.

FUTILITY

Binary indicator of whether a futility decision will be allowed.

nTreat

Number of treatments in consideration, i.e. K.

nCat

Number of ordinal outcome categories, i.e. J.

UT

Vector of numerical utility scores to give outcomes 1,...,J.

DeltaVEC

Vector of \Delta values to test.

gamma

Length 3 vector of cutoff parameters.

PSPIKE

Prior probability of a pairwise null. PSPIKE=1 means no clustering is possible.

ADJ

Integer for whether or not we should adjust for covariates.

B

Number of MCMC iterations to perform.

Value

The set of active treatments to continue, an optimal treatment, or a set of equally optimal treatments. Also reports posterior mean utilities and ordinal outcome probabilities as well as pairwise comparisons of utility similarity, when appropriate.

References

[1] Chapple and Clement (2020), MABOUST: A Multi-Armed Bayesian Ordinal Outcome Utility-Based Sequential Trial. Submitted.

Examples

##Clinical Parameters
nCat = 6
nTreat = 3
UT = c(0,10,20,80,90,100)
DeltaVEC  = c(5,10)
###Which treatments are active?
ACTIVE = c(1,0,1) ###Treatments 1, 3 are active
FUTILITY = 1 ###Futility look is allowed.
###Design parameters
gamma= c(.5, .05, .05)
PSPIKE = .9
set.seed(1)
##Generate Random Data
n=300
Y=sample(1:nCat,n,replace=TRUE)
T1 = sample(1:nTreat,n,replace=TRUE)
X=matrix(rnorm(n*2),ncol=2)
###Number of iterations
B=100
PSPIKE = .9
ADJ = 1
MABOUST(Y, T1, X, ACTIVE, FUTILITY, nTreat, nCat, UT, DeltaVEC, gamma, PSPIKE, ADJ,B )