Nothing
## ---- echo=FALSE--------------------------------------------------------------
knitr::kable(data.frame(
i = 1:6,
Pretreated = c(F, F, F, T, T, T),
PDL1 = c('Low', 'Medium', 'High', 'Low', 'Medium', 'High'),
x1 = rep(c(0, 1), each = 3),
x2 = c(1, 0, 0, 1, 0, 0),
x3 = c(0, 1, 0, 0, 1, 0)
))
## ----message=FALSE------------------------------------------------------------
library(trialr)
peps2_sc <- function() peps2_get_data(num_patients = 60,
prob_eff = c(0.167, 0.192, 0.5, 0.091, 0.156, 0.439),
prob_tox = rep(0.1, 6),
eff_tox_or = rep(1, 6))
set.seed(123)
dat <- peps2_sc()
## -----------------------------------------------------------------------------
c(dat$alpha_mean, dat$alpha_sd)
## -----------------------------------------------------------------------------
knitr::kable(
head(with(dat, data.frame(eff, tox, x1, x2, x3)), 10)
)
## ---- results = "hide"--------------------------------------------------------
fit <- stan_peps2(dat$eff, dat$tox, dat$cohorts)
## ---- message = FALSE, fig.width = 7, fig.height = 7, fig.cap = "Posterior Prob(Efficacy) in the six PePS2 cohorts"----
rstan::plot(fit, pars = 'prob_eff')
## ---- eval=TRUE---------------------------------------------------------------
decision <- peps2_process(fit)
knitr::kable(
with(decision, data.frame(ProbEff, ProbAccEff, ProbTox, ProbAccTox, Accept)),
digits = 3
)
## ---- eval=FALSE--------------------------------------------------------------
# set.seed(123)
# run_sims <- function(num_sims = 10, sample_data_func = peps2_sc,
# summarise_func = peps2_process, ...) {
# sims <- list()
# for(i in 1:num_sims) {
# print(i)
# dat <- sample_data_func()
# fit <- stan_peps2(dat$eff, dat$tox, dat$cohorts, ...)
# sim <- summarise_func(fit)
# sims[[i]] <- sim
# }
# return(sims)
# }
#
# sims <- run_sims(num_sims = 10, sample_data_func = peps2_sc,
# summarise_func = peps2_process)
## ---- eval=FALSE--------------------------------------------------------------
# apply(sapply(sims, function(x) x$Accept), 1, mean)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.