setwd("c:\\E-drive\\Books\\LOGLIN3\\BAYES\\")
load("post-samp.Rda")
# load("post-samp-18.Rda")
#create a placeholder matrix for the output
out=rep(1,15)
out=matrix(out,ncol=3)

af=c(.9,.75,.5,.25,.1)
for(k in 1:5){
lda=(log(af[k]/(1-af[k]))-beta[1,])/beta[2,]
out[k,]=quantile(lda,c(.05,.5,.95))
}
out