data("AlcoholUse", package = "zoib")
AlcoholUse$Grade = as.factor(AlcoholUse$Grade)
post.obj <- zoib(Percentage ~ Grade+Days+Gender|1|Grade+Days+Gender|1,
data = AlcoholUse, random = 1, EUID= AlcoholUse$County,
zero.inflation = TRUE, one.inflation = FALSE, joint = FALSE,
n.iter=1000, n.thin=5)
post.sample <- post.obj$oripara
post.sample.c1<- post.sample[[1]][11:200,]
post.sample.c2<- post.sample[[2]][11:200,]
post.sample <- mcmc.list(as.mcmc(post.sample.c1),as.mcmc(post.sample.c2))
summary(post.sample)Run the code above in your browser using DataLab