library("DBR")
library("coda")
data("pain")
est <- dbr(
interference ~ severity + age
, pain
, control = dbr.control(
nsmp = 50
, nburnin = 25
)
)
coda_wrapper(est, autocorr, lag = c(0, 5))
Run the code above in your browser using DataLab