normalBACT: Normal distribution for Bayesian Adaptive trials

In thevaachandereng/BACT: Simulation and Analysis of Adaptive Bayesian Clinical Trials

Description

Simulation of continuous (normally distributed) data for Bayesian adaptive trials with various inputs to control for power, sample size, type I error rate, etc.

Usage

normalBACT(
  mu_treatment,
  sd_treatment,
  mu_control = NULL,
  sd_control = NULL,
  mu0_treatment = NULL,
  sd0_treatment = NULL,
  N0_treatment = NULL,
  mu0_control = NULL,
  sd0_control = NULL,
  N0_control = NULL,
  N_total,
  lambda = 0.3,
  lambda_time = NULL,
  interim_look = NULL,
  EndofStudy,
  block = 2,
  rand_ratio = c(1, 1),
  discount_function = "identity",
  alternative = "greater",
  prop_loss_to_followup = 0.15,
  h0 = 0,
  futility_prob = 0.05,
  expected_success_prob = 0.9,
  prob_ha = 0.95,
  N_impute = 10,
  number_mcmc = 10000,
  alpha_max = 1,
  fix_alpha = FALSE,
  weibull_scale = 0.135,
  weibull_shape = 3,
  method = "fixed"
)

Arguments

mu_treatment	scalar. Mean outcome in the treatment arm.
sd_treatment	scalar. Standard deviation of outcome in the treatment.
mu_control	scalar. Mean outcome in the control arm.
sd_control	scalar. Standard deviation of outcome in the control arm. arm.
mu0_treatment	scalar. Mean of the historical treatment group.
sd0_treatment	scalar. Standard deviation of the historical treatment group.
N0_treatment	scalar. Number of observations of the historical treatment group.
mu0_control	scalar. Mean of the historical control group.
sd0_control	scalar. Standard deviation of the historical control group.
N0_control	scalar. Number of observations of the historical control group.
N_total	scalar. Total sample size.
lambda	vector. Enrollment rates across simulated enrollment times. See enrollment for more details.
lambda_time	vector. Enrollment time(s) at which the enrollment rates change. Must be same length as lambda. See enrollment for more details.
interim_look	vector. Sample size for each interim look. Note: the maximum sample size should not be included.
EndofStudy	scalar. Length of the study.
block	scalar. Block size for generating the randomization schedule.
rand_ratio	vector. Randomization allocation for the ratio of control to treatment. Integer values mapping the size of the block. See randomization for more details.
discount_function	character. If incorporating historical data, specify the discount function. Currently supports the Weibull function (discount_function = "weibull"), the scaled-Weibull function (discount_function = "scaledweibull"), and the identity function (discount_function = "identity"). The scaled-Weibull discount function scales the output of the Weibull CDF to have a maximum value of 1. The identity discount function uses the posterior probability directly as the discount weight. Default value is "identity". See bdpnormal for more details.
alternative	character. The string specifying the alternative hypothesis, must be one of "greater" (default), "less" or "two.sided".
prop_loss_to_followup	scalar. Overall proportion of subjects lost to follow-up.
h0	scalar. Threshold for comparing two mean values. Default is h0 = 0.
futility_prob	scalar. Probability of stopping early for futility.
expected_success_prob	scalar. Probability of stopping early for success.
prob_ha	scalar. Probability of alternative hypothesis.
N_impute	scalar. Number of imputations for Monte Carlo simulation of missing data.
number_mcmc	scalar. Number of Markov Chain Monte Carlo draws in sampling posterior.
alpha_max	scalar. Maximum weight the discount function can apply. Default is 1. For a two-arm trial, users may specify a vector of two values where the first value is used to weight the historical treatment group and the second value is used to weight the historical control group.
fix_alpha	logical. Fix alpha at alpha_max? Default value is FALSE.
weibull_scale	scalar. Scale parameter of the Weibull discount function used to compute alpha, the weight parameter of the historical data. Default value is 0.135. For a two-arm trial, users may specify a vector of two values where the first value is used to estimate the weight of the historical treatment group and the second value is used to estimate the weight of the historical control group. Not used when discount_function = "identity".
weibull_shape	scalar. Shape parameter of the Weibull discount function used to compute alpha, the weight parameter of the historical data. Default value is 3. For a two-arm trial, users may specify a vector of two values where the first value is used to estimate the weight of the historical treatment group and the second value is used to estimate the weight of the historical control group. Not used when discount_function = "identity".
method	character. Analysis method with respect to estimation of the weight parameter alpha. Default method "mc" estimates alpha for each Monte Carlo iteration. Alternate value "fixed" estimates alpha once and holds it fixed throughout the analysis. See the the bdpsurvival vignette vignette("bdpsurvival-vignette", package="bayesDP") for more details.

Value

A list of output for a single trial simulation:

mu_treatment

scalar. The input parameter of mean value of the outcome in the treatment group.

p_control

scalar. The input parameter of mean value of the outcome in the control group.

sd_treatment

scalar. The input parameter of standard deviation of the outcome in the control group.

sd_control

scalar. The input parameter of standard deviation of the outcome in the control group.

prob_of_accepting_alternative

scalar. The input parameter of probability threshold of accepting the alternative.

margin

scalar. The margin input value of difference between mean estimate of treatment and mean estimate of the control.

alternative

character. The input parameter of alternative hypothesis.

interim_look

vector. The sample size for each interim look.

N_treatment

scalar. The number of patients enrolled in the experimental group for each simulation.

N_control

scalar. The number of patients enrolled in the control group for each simulation.

N_enrolled

vector. The number of patients enrolled in the trial (sum of control and experimental group for each simulation).

N_complete

scalar. The number of patients who completed the trial and had no loss to follow-up.

post_prob_accept_alternative

vector. The final probability of accepting the alternative hypothesis after the analysis is done.

est_final

scalar. The final estimate of the difference in posterior estimate of treatment and posterior estimate of the control group.

stop_futility

scalar. Did the trial stop for futility during imputation of patient who had loss to follow up? 1 for yes and 0 for no.

stop_expected_success

scalar. Did the trial stop for early success during imputation of patient who had loss to follow up? 1 for yes and 0 for no.

est_interim

scalar. The interim estimate of the difference in posterior estimate of treatment and posterior estimate of the control group.

thevaachandereng/BACT documentation built on July 24, 2020, 2:35 a.m.