# NOT RUN {
# standard deviation of measurments in one group
stDev <- 24
# sample size at final analysis
n1 <- 92
n2 <- 92
# MDD at final analysis (corresponds to delta = 10 for 80% power)
mdd <- 7.023506
# prior
priormean <- 12.3
# standard error for prior, based on Phase 2 data
sig1 <- 26.1
n1p <- 25
sig2 <- 33.6
n2p <- 25
sd0 <- sqrt(sig1 ^ 2 / n1p + sig2 ^ 2 / n2p)
# flat prior
width1 <- 25
height1 <- 0.02
# bpps
bpp_continuous(prior = "normal", successmean = mdd, stDev = stDev,
n1 = n1, n2 = n2, priormean = priormean, priorsigma = sd0)
bpp_continuous(prior = "flat", successmean = mdd, stDev = stDev,
n1 = n1, n2 = n2, priormean = priormean,
width = width1, height = height1)
# }
Run the code above in your browser using DataLab