source("binomialOneArm.R")


testrunNull <- binomialOneArm(nreps = 10000, nmax = 19, futilitybound = 0.05, checkpoints = c(10, 13, 16),
           postcut = 0.95, betapriorparam1 = 0.15, betapriorparam2 = 0.85, p_true = 0.15, p_0 = 0.15, seed = 42)
testrunAlt1 <- binomialOneArm(nreps = 10000, nmax = 19, futilitybound = 0.05, checkpoints = c(10, 13, 16),
           postcut = 0.95, betapriorparam1 = 0.15, betapriorparam2 = 0.85, p_true = 0.35, p_0 = 0.15, seed = 42)
testrunAlt2 <- binomialOneArm(nreps = 10000, nmax = 19, futilitybound = 0.05, checkpoints = c(10, 13, 16),
           postcut = 0.95, betapriorparam1 = 0.15, betapriorparam2 = 0.85, p_true = 0.45, p_0 = 0.15, seed = 42)
testrunNull$output
testrunAlt1$output
testrunAlt2$output





testrunNull <- binomialOneArm(nreps = 10000, nmax = 19, futilitybound = 0.05, checkpoints = c(10, 13, 16),
           postcut = 0.97, betapriorparam1 = 0.15, betapriorparam2 = 0.85, p_true = 0.15, p_0 = 0.15, seed = 42)
testrunAlt1 <- binomialOneArm(nreps = 10000, nmax = 19, futilitybound = 0.05, checkpoints = c(10, 13, 16),
           postcut = 0.97, betapriorparam1 = 0.15, betapriorparam2 = 0.85, p_true = 0.35, p_0 = 0.15, seed = 42)
testrunAlt2 <- binomialOneArm(nreps = 10000, nmax = 19, futilitybound = 0.05, checkpoints = c(10, 13, 16),
           postcut = 0.97, betapriorparam1 = 0.15, betapriorparam2 = 0.85, p_true = 0.45, p_0 = 0.15, seed = 42)
testrunNull$output
testrunAlt1$output
testrunAlt2$output
 


testrunNull <- binomialOneArm(nreps = 10000, nmax = 19, futilitybound = 0.15, checkpoints = c(10, 13, 16),
           postcut = 0.95, betapriorparam1 = 0.15, betapriorparam2 = 0.85, p_true = 0.15, p_0 = 0.15, seed = 42)
testrunAlt1 <- binomialOneArm(nreps = 10000, nmax = 19, futilitybound = 0.15, checkpoints = c(10, 13, 16),
           postcut = 0.95, betapriorparam1 = 0.15, betapriorparam2 = 0.85, p_true = 0.35, p_0 = 0.15, seed = 42)
testrunAlt2 <- binomialOneArm(nreps = 10000, nmax = 19, futilitybound = 0.15, checkpoints = c(10, 13, 16),
           postcut = 0.95, betapriorparam1 = 0.15, betapriorparam2 = 0.85, p_true = 0.45, p_0 = 0.15, seed = 42)
testrunNull$output
testrunAlt1$output
testrunAlt2$output