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))
total_weight_at_dose(fit) # c(2, 1.1, 0, 0)
total_weight_at_dose(fit, dose = 2) # 1.1
}
Run the code above in your browser using DataLab