mydata <- data_dose
BJSM_dose_result <- BJSM_binary(
data = data_dose, prior_dist = c("beta", "gamma"),
pi_prior = c(3, 17), normal.par = c(0.2, 100), beta_prior = c(2, 2),
n_MCMC_chain = 2, n.adapt = 100, MCMC_SAMPLE = 2000, ci = 0.95
)
Run the code above in your browser using DataLab