MCBUpperLimits {SMARTbayesR}

Simultaneous Upper Credible Intervals

Description

Compute simultaneous upper credible intervals from draws of embedded dynamic treatment regime probabilities given by the function argument “thetadraws”.

Usage

MCBUpperLimits(thetadraws, alpha = 0.05, design = "design-1", type = "log-OR")

Arguments

thetadraws

draws of the embedded dynamic treatment regimes.

alpha

the probability of excluding the true best EDTR from the set of best.

design

specifies to which SMART design to apply function: either design-1, general, or design-3.

type

summary statistic: log-OR, log-RR, or RD

Value

Upper 1-alpha level simultaneous credible interval limits for the embedded dynamic treatment regimes.

Examples

dat <- SimDesign1(sample_size=250,
                  response_prob = c(0.5,0.9,0.3,0.7,0.5,0.8),
                  stage_one_trt_one_response_prob = 0.7,
                  stage_one_trt_two_response_prob = 0.4)

x <- PosteriorTrtSeqProb(niter = 1000, dat, design = "design-1")

thetadraws <- PosteriorEDTRProbs(x, design = "design-1")

MCBUpperLimits(thetadraws,
               alpha = 0.05,
               design = "design-1",
               type = "log-OR")