if (FALSE) {
# CRM example
fit <- stan_crm(skeleton = c(0.1, 0.2, 0.35, 0.6), target = 0.2,
model = 'empiric', beta_sd = sqrt(1.34), seed = 123,
doses = c(1, 1, 2, 2, 2),
tox = c(0, 0, 0, 0, 0),
weights = c(1, 1, 0.9, 0.1, 0.1))
l <- weights_at_dose(fit)
length(l) # 4
l[[1]] # c(1, 1)
l[[2]] # c(0.9, 0.1, 0.1)
l[[3]] # c()
weights_at_dose(fit, dose = 2) # c(0.9, 0.1, 0.1)
}
Run the code above in your browser using DataLab