diff --git a/R/bootstrapping.R b/R/bootstrapping.R index b73bd4c..db97781 100644 --- a/R/bootstrapping.R +++ b/R/bootstrapping.R @@ -9,6 +9,7 @@ #' @return tbd #' @export getBootsrapBands <- function ( + model_fits, n_samples = 1e3, alpha = c(0.05, 0.5), @@ -16,9 +17,10 @@ getBootsrapBands <- function ( dose_seq = NULL ) { + mu_hat_samples <- sapply(attr(model_fits, "posterior"), RBesT::rmix, n = n_samples) - sd_hat <- summary.postList(model_fits)[, 2] + sd_hat <- summary.postList(attr(model_fits, "posterior"))[, 2] dose_levels <- model_fits[[1]]$dose_levels model_names <- names(model_fits)