Test a cer design for early rejection of hypotheses and calculate the CER for adaptions

Usage

cer_final_test(design, p_values)

Arguments

design: A cer_design object
p_values: A list of p-values for the hypotheses

Value

a cer_design object, which now also includes the rejection status of the hypotheses after the final test

Examples

as <- function(x,t) 2-2*pnorm(qnorm(1-x/2)/sqrt(t))
design <- cer_design(
 correlation=rbind(H1=c(1, NA),
                   H2=c(NA, 1)),
 weights=c(2/3, 1/3),
 alpha=0.05,
 test_m=rbind(c(0, 1),
              c(1, 0)),
 alpha_spending_f=as,
 t=0.5)

design <- cer_interim_test(design, c(0.001, 0.02))

design <- cer_drop_hypotheses(design, c(TRUE, FALSE))
design <- cer_adapt_bounds(design)

design <- cer_final_test(design, c(NA, 0.01))
design$rej
#> [1] TRUE TRUE