Learn R Programming

bpp (version 1.0.4)

bpp_binary: Bayesian Predictive Power (BPP) for Binary Endpoint

Description

Compute BPP for a binary endpoint.

Usage

bpp_binary(prior = c("normal", "flat"), successdelta, pi1, n1,
                       pi2, n2, priormean, ...)

Arguments

prior

Prior density on effect sizes.

successdelta

The proportion difference that defines success at the final analysis. We assume that higher proportions are better, e.g. as for response in oncology. Typically chosen to be the minimal detectable difference, i.e. the critical on the scale of the effect size of interest corresponding to the significance level at the final analysis.

pi1

Assumed response proportion in intervention arm. Used to compute standard error at final analysis.

n1

Sample size in intervention arm. Used to compute standard error at final analysis.

pi2

Assumed response proportion in control arm. Used to compute standard error at final analysis.

n2

Sample size in control arm. Used to compute standard error at final analysis.

priormean

Prior mean.

...

Further arguments specific to the chosen prior (see bpp_binary for examples).

Value

A real number, the bpp.

References

Rufibach, K., Jordan, P., Abt, M. (2016a). Sequentially Updating the Likelihood of Success of a Phase 3 Pivotal Time-to-Event Trial based on Interim Analyses or External Information. J. Biopharm. Stat., 26(2), 191--201.

Rufibach, K., Burger, H.U., Abt, M. (2016b). Bayesian Predictive Power: Choice of Prior and some Recommendations for its Use as Probability of Success in Drug Development. Pharm. Stat., 15, 438--446.

Examples

Run this code
# NOT RUN {
# simple example with response proportions
pi1 <- 0.6              # proportion in intervention arm
pi2 <- 0.45             # proportion in control arm
n1 <- 174
n2 <- 174

# MDD at final analysis - proportion difference that corresponds to "success"
mdd <- 0.1

# prior if normal
pi20 <- 0.44
pi10 <- 0.64
n0 <- 50
priormean <- pi10 - pi20
sd0 <- sqrt(pi20 * (1 - pi20) / (n0 / 2) + pi10 * (1 - pi10) / (n0 / 2))

bpp0 <- bpp_binary(prior = "normal", successdelta = mdd, pi1 = pi1, n1 = n1,
           pi2 = pi2, n2 = n2, priormean = priormean, priorsigma = sd0)
bpp0

# prior if flat
width1 <- 0.5
height1 <- 1.5

bpp0_1 <- bpp_binary(prior = "flat", successdelta = mdd, pi1 = pi1, n1 = n1,
                    pi2 = pi2, n2 = n2, priormean = priormean, 
                    width = width1, height = height1)
bpp0_1
# }

Run the code above in your browser using DataLab