
Exact repeated confidence intervals for vaccine or prevention efficacy
Source:R/ciBinomialExact.R
repeatedCIBinomialExact.RdInverts exact binomial efficacy tests at each completed analysis to obtain repeated confidence intervals for vaccine or prevention efficacy.
Usage
repeatedCIBinomialExact(
gsD,
n.I = NULL,
x = NULL,
conf.level = 0.95,
tol = 1e-08,
maxiter = 100
)Arguments
- gsD
A
gsSurvobject with non-bindingtest.type1, 4, 6, or 8.- n.I
Increasing integer total event counts at completed analyses. If
NULL, planned integer event counts fromtoInteger(gsD)are used.- x
Integer experimental-arm event counts at the analyses in
n.I.- conf.level
Two-sided confidence level.
- tol
Absolute tolerance for bisection on the conditional binomial event-probability scale.
- maxiter
Maximum bisection iterations for each confidence limit.
Value
A data frame with one row per completed analysis containing the observed counts, efficacy estimate, repeated confidence limits, confidence level, and one-sided tail level.
Details
A two-sided interval with confidence level 1 - alpha uses the same
spending function, spending times, and count-path ordering in both
directions, with alpha / 2 in each tail. The lower efficacy limit inverts
the usual lower event-count efficacy test. The upper efficacy limit inverts
its mirror image after exchanging experimental- and control-arm event
counts. This mirrored test is not the design's futility boundary.
Non-binding futility and harm are ignored in both directions. Coverage is generally conservative because the exact rejection regions are discrete. Spending time remains relative to the planned final event count.
This follows the repeated-confidence-interval construction of Jennison and Turnbull (1984), using exact Bernoulli ordering as in Coe and Tamhane (1993).
References
Jennison, C. and Turnbull, B. W. (1984). Repeated confidence intervals for group sequential clinical trials. Controlled Clinical Trials, 5, 33–45.
Coe, P. R. and Tamhane, A. C. (1993). Exact repeated confidence intervals for Bernoulli parameters in a group sequential clinical trial. Controlled Clinical Trials, 14, 19–29.
Examples
design <- gsSurv(
k = 3, test.type = 4, timing = c(.45, .7), ratio = 3,
hr = .3, hr0 = .7
)
counts <- toBinomialExact(design)$n.I
# \donttest{
repeatedCIBinomialExact(design, counts, x = c(12, 23, 38))
#> Analysis n.I x estimate conf.low conf.high conf.level tail_alpha
#> 1 1 31 12 0.7894737 0.3792904 0.9345737 0.95 0.025
#> 2 2 48 23 0.6933333 0.3477728 0.8562642 0.95 0.025
#> 3 3 68 38 0.5777778 0.2902891 0.7471530 0.95 0.025
# }