#\examples<>

#<ss.cons.2norm.avg.hpdlimits>
accepted.diff <- 0.01
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.avg.hpdlimits(accepted.diff, prior1, prior2)
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), known.sd=c(1, 1.2))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), known.sd=c(1.1, 1.4))
z <- ss.cons.2norm.avg.hpdlimits(accepted.diff, prior1, prior2)

#<ss.cons.2norm.avg.hpdlimits.mymarg>
accepted.diff <- 0.01
clinical.prior <- list(mu0=c(10, 10), n0=c(10, 9),  prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.avg.hpdlimits.mymarg(accepted.diff, prior1, prior2, clinical.prior)

#<ss.cons.2norm.avg.hpdlimits.bothmarg>
accepted.diff <- 0.01
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.avg.hpdlimits.bothmarg(accepted.diff, prior1, prior2)
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), known.sd=c(1, 1.2))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), known.sd=c(1.1, 1.4))
z <- ss.cons.2norm.avg.hpdlimits.bothmarg(accepted.diff, prior1, prior2)

#<ss.cons.2norm.prob.hpdlimits>
accepted.diff <- 0.01
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.prob.hpdlimits(accepted.diff, prior1, prior2, prob=0.75)

#<ss.cons.2norm.prob.hpdlimits.mymarg>
accepted.diff <- 0.01
clinical.prior <- list(mu0=c(10, 10), n0=c(10, 9),  prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.prob.hpdlimits.mymarg(accepted.diff, prior1, prior2, clinical.prior, prob=0.75)

#<ss.cons.2norm.prob.hpdlimits.bothmarg>
accepted.diff <- 0.01
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.prob.hpdlimits.bothmarg(accepted.diff, prior1, prior2, prob=0.75)

#<ss.cons.2norm.avg.cdf>
accepted.cdf.diff <- 0.01
cdf.points <- c(-1, 1)
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.avg.cdf(cdf.points, accepted.cdf.diff, prior1, prior2)

#<ss.cons.2norm.avg.cdf.mymarg>
accepted.cdf.diff <- 0.01
cdf.points <- c(-1, 1)
clinical.prior <- list(mu0=c(10, 10), n0=c(10, 9),  prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.avg.cdf.mymarg(cdf.points, accepted.cdf.diff, prior1, prior2, clinical.prior)

#<ss.cons.2norm.avg.cdf.bothmarg>
accepted.cdf.diff <- 0.01
cdf.points <- c(-1, 1)
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.avg.cdf.bothmarg(cdf.points, accepted.cdf.diff, prior1, prior2)

#<ss.cons.2norm.prob.cdf>
accepted.cdf.diff <- 0.01
cdf.points <- c(-1, 1)
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.prob.cdf(cdf.points, accepted.cdf.diff, prior1, prior2, prob=0.85)

#<ss.cons.2norm.prob.cdf.mymarg>
accepted.cdf.diff <- 0.01
cdf.points <- c(-1, 1)
clinical.prior <- list(mu0=c(10, 10), n0=c(10, 9),  prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.prob.cdf.mymarg(cdf.points, accepted.cdf.diff, prior1, prior2, clinical.prior, prob=0.85)

#<ss.cons.2norm.prob.cdf.bothmarg>
accepted.cdf.diff <- 0.01
cdf.points <- c(-1, 1)
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
z <- ss.cons.2norm.prob.cdf.bothmarg(cdf.points, accepted.cdf.diff, prior1, prior2)

#<ss.cons.2norm.avg.q>
accepted.diff <- 0.01
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
quantiles <- c(0.2, 0.75)
z <- ss.cons.2norm.avg.q(quantiles, accepted.diff, prior1, prior2)

#<ss.cons.2norm.avg.q.mymarg>
accepted.diff <- 0.01
clinical.prior <- list(mu0=c(10, 10), n0=c(10, 9),  prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
quantiles <- c(0.2, 0.75)
z <- ss.cons.2norm.avg.q.mymarg(quantiles, accepted.diff, prior1, prior2, clinical.prior)

#<ss.cons.2norm.avg.q.bothmarg>
accepted.diff <- 0.01
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
quantiles <- c(0.2, 0.75)
z <- ss.cons.2norm.avg.q.bothmarg(quantiles, accepted.diff, prior1, prior2)

#<ss.cons.2norm.prob.q>
accepted.diff <- 0.01
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
quantiles <- c(0.2, 0.75)
z <- ss.cons.2norm.prob.q(quantiles, accepted.diff, prior1, prior2, prob=0.77)

#<ss.cons.2norm.prob.q.mymarg>
accepted.diff <- 0.01
clinical.prior <- list(mu0=c(10, 10), n0=c(10, 9),  prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
quantiles <- c(0.2, 0.75)
z <- ss.cons.2norm.prob.q.mymarg(quantiles, accepted.diff, prior1, prior2, clinical.prior, prob=0.77)

#<ss.cons.2norm.prob.q.bothmarg>
accepted.diff <- 0.01
prior1 <- list(mu0=c(15, 20), n0=c(25, 25), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
prior2 <- list(mu0=c(5, 6),   n0=c(25, 30), prec.shape=c(23, 25), prec.rate=c(2100, 2000))
quantiles <- c(0.2, 0.75)
z <- ss.cons.2norm.prob.q.bothmarg(quantiles, accepted.diff, prior1, prior2, prob=0.77)
