# NOT RUN {
nDat = c(25, 25, 25, 25) # total number of patients
xDat = c(2, 3, 8, 6) # number of responses
alpha <- 1e-20
d0 <- 0.0
alpha1 = 50
beta1 = 10
tau2 <- 0.1
phi1 <- 0.1
deltaT <- 0.2
thetaT <- 0.60
res <- BCHM(nDat = nDat,
xDat = xDat,
alpha = alpha,
d0 = d0,
alpha1 = alpha1,
beta1 = beta1,
tau2 = tau2,
phi1 = phi1,
deltaT = deltaT,
thetaT = thetaT,
burnIn = 100,
MCIter = 200,
MCNum = 1000,
seed = 1000
)
print(res$SMatrix)
print(res$Result)
col <- res$Result[,4]
BCHMplot_post_value(res, col, HPD = 0.8)
# }
Run the code above in your browser using DataLab