data("BiRepeated", package = "zoib")
post.obj <- zoib(y1|y2 ~ x|1|x, data= BiRepeated,
random=1, EUID= BiRepeated$id,
zero.inflation = FALSE, one.inflation = FALSE,
prior.Sigma = "UN.unif", n.iter=4000, n.thin=20)
post.sample <- post.obj$oripara
post.sample.c1 <- post.sample[[1]][51:200,]
post.sample.c2 <- post.sample[[2]][51:200,]
post.sample <- mcmc.list(mcmc(post.sample.c1),mcmc(post.sample.c2))
summary(post.sample)Run the code above in your browser using DataLab