#\dontrun{
data("GasolineYield", package = "zoib")
GasolineYield$batch <- as.factor(GasolineYield$batch)
# fixed effects zoib with batch treated as a 10-level qualitative variable
eg.fixed <- zoib(yield ~ temp + as.factor(batch)| 1, data=GasolineYield,
joint = FALSE, random = 0, EUID = 1:nrow(GasolineYield),
zero.inflation = FALSE, one.inflation = FALSE,
n.iter = 100, n.thin = 5)
post.sample<- eg.fixed$oripara
traceplot(post.sample)
autocorr.plot(post.sample)
gelman.diag(post.sample)
summary(post.sample)
#}Run the code above in your browser using DataLab