# NOT RUN {
#simulate data
ngibbs<- 1000
y<- (-1000 * 501:1500)/(-500 + 501:1500) + rnorm(ngibbs, 0, 0.1)
dat<- matrix(c(1, y), 1, 1001)
dat<- data.frame(dat)
names(dat)[1]<- "id"
#run function
traceplot(data = dat, ngibbs = ngibbs, type = "LML")
# }
Run the code above in your browser using DataLab