getBootstrapQuantiles: getBootstrapQuantiles

In BayesianMCPMod: Simulate, Evaluate, and Analyze Dose Finding Trials with Bayesian MCPMod

getBootstrapQuantiles

Description

Calculates quantiles from bootstrapped dose predictions. Can be used to derive credible intervals to assess the uncertainty for the model fit.

Usage

getBootstrapQuantiles(bs_samples, quantiles)

Arguments

bs_samples	An object of class bootstrappedSample as created by getBootstrapSamples
quantiles	A vector of quantiles that should be evaluated

Value

A data frame with entries doses, models, and quantiles

Examples

posterior_list <- list(Ctrl = RBesT::mixnorm(comp1 = c(w = 1, m = 0, s = 1), sigma = 2),
                       DG_1 = RBesT::mixnorm(comp1 = c(w = 1, m = 3, s = 1.2), sigma = 2),
                       DG_2 = RBesT::mixnorm(comp1 = c(w = 1, m = 4, s = 1.5), sigma = 2) ,  
                       DG_3 = RBesT::mixnorm(comp1 = c(w = 1, m = 6, s = 1.2), sigma = 2) ,
                       DG_4 = RBesT::mixnorm(comp1 = c(w = 1, m = 6.5, s = 1.1), sigma = 2))
models         <- c("exponential", "linear")
dose_levels    <- c(0, 1, 2, 4, 8)
fit            <- getModelFits(models      = models,
                               posterior   = posterior_list,
                               dose_levels = dose_levels,
                               simple      = TRUE)
                               
bs_samples     <- getBootstrapSamples(model_fits = fit,
                                      n_samples  = 10, # speeding up example run time
                                      doses      = c(0, 6, 8))
                           
getBootstrapQuantiles(bs_samples = bs_samples,
                      quantiles  = c(0.025, 0.5, 0.975))

BayesianMCPMod documentation built on May 29, 2024, 9:14 a.m.