library("DBR")
data("pain")
est <- dbr(
interference ~ severity + age
, pain
, control = dbr.control(
nsmp = 50
, nburnin = 25
)
)
print(est)
plot(est)
coef(est, prob = c(0.05, 0.5, 0.95))
summary(est, context = pain[2, ])
Run the code above in your browser using DataLab