get.Q.bwer | R Documentation |
Get Efficacy Cutoff based on Basket-Wise Error Rate (BWER) Control
get.Q.bwer(object, alpha = 0.1, digits = 3, Qclust = NULL)
object |
returned by post.infer. |
alpha |
basket-wise type I error control. |
digits |
number of digits in the cutoffs. |
Qclust |
|
It returns the efficacy cutoffs.
N <- rbind(
c(10, 25),
c(10, 25),
c(10, 25),
c(10, 25),
c(10, 25)) # interim sample size and total sample size for each indication
scenarios <- rbind( c(0.15, 0.15, 0.15, 0.15, 0.15), c(0.3, 0.3, 0.3, 0.3, 0.3) )
res <- generate.data(N = N, ORRs = scenarios, ntrial = 1000, seed = 343809)
post <- post.infer(res, pnull = rep(0.15,5), stopbounds = cbind(c(1,1,1,1,1)),
ModelFit = "localPP", method = "PEB", a = 2, delta = 0.3)
(Q <- get.Q.bwer(post, alpha = 0.1, digits = 3, Qclust = rep(1, 5)))
Qmat <- array(NA, dim = dim(post$postprob))
for(i in 1:5) Qmat[,,i] <- Q[i]
apply(post$postprob>Qmat, c(1,3), mean)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.