# NOT RUN {
# standard deviation of measurments in one group
stDev <- 24
# number of patients at interim
n1_int <- 46
n2_int <- 46
# number of patients at final analysis
n1 <- 92
n2 <- 92
# MDD at final analysis (corresponds to delta = 10 for 80% power)
mdd <- 7.023506
# efficacy boundary
effi <- 15
# futility boundary --> chosen informally
futi <- 0
# prior
priormean <- 12.3
sd0 <- 4.2
# flat prior
width1 <- 25
height1 <- 0.02
thetas <- seq(-0.65, 0.3, by = 0.01)
bpp_1interim_continuous(prior = "normal", successmean = mdd, stDev = stDev,
n1 = c(n1_int, n1), n2 = c(n2_int, n2),
IntEffBoundary = effi, IntFutBoundary = futi, IntFix = 1,
priormean = priormean, propA = 0.5, thetas,
priorsigma = sd0)[[1]]
bpp_1interim_continuous(prior = "flat", successmean = mdd, stDev = stDev,
n1 = c(n1_int, n1), n2 = c(n2_int, n2), IntEffBoundary = effi,
IntFutBoundary = futi, IntFix = 1,
priormean = 12.3, propA = 0.5, thetas = thetas,
width = width1, height = height1)[[1]]
# }
Run the code above in your browser using DataLab