MCBUpperLimits: Simultaneous Upper Credible Intervals

In SMARTbayesR: Bayesian Set of Best Dynamic Treatment Regimes and Sample Size in SMARTs for Binary Outcomes

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")

SMARTbayesR documentation built on Oct. 1, 2021, 1:06 a.m.