rm(list = ls())
setwd("c:\\E-drive\\Books\\LOGLIN3\\BAYES\\")
load("post-samp-thin.Rda")
post_slope=density(beta[2,])
plot(post_slope,main="")