# NOT RUN {
pi1 <- pi_MCSim_beta(M = 5000, seed = 10, a = 1, b = 1)
pi2 <- pi_MCSim_beta(M = 5000, seed = 10, a = 2, b = 18)
MFM_optimal_n(cvec = c(1, 1), pivec = cbind(pi1, pi2), Rvec = c(0.8, 0.7), thres_CR = 0.05)
# }
Run the code above in your browser using DataLab