x <- conjugate(
s1 = rnorm(10, 50, 10), s2 = rnorm(10, 60, 12), method = "t",
priors = list(list(mu = 40, sd = 10), list(mu = 45, sd = 8)),
rope_range = c(-5, 8), rope_ci = 0.89,
cred.int.level = 0.89, hypothesis = "equal"
)
plot(x)
Run the code above in your browser using DataLab