## ----packages, message=FALSE-----------------------------------------------
# load CNPBayes
suppressMessages(library(CNPBayes))
# load packages for manipulating and visualizing data
suppressMessages(library(dplyr))
suppressMessages(library(tidyr))
suppressMessages(library(ggplot2))

## ----post-sim--------------------------------------------------------------
set.seed(1)
N <- 200
n <- 81
lrr <- c(rnorm(100, -0.5, sd=0.1), rnorm(100, 0, sd=0.1))
mp <- McmcParams(iter=600, burnin=10, nStarts=4)
mlist <- gibbs(model="SBP", mp=mp, dat=lrr, k_range=c(2, 2))

## ----chains----------------------------------------------------------------
model <- mlist[[1]]
figs <- ggChains(model)
figs[[1]]
figs[[2]]

## ----marginal-lik, eval=FALSE----------------------------------------------
#  m.lik <- round(sapply(mlist, marginal_lik), 1)
#  m.lik