rm(list = ls())
setwd("c:\\E-drive\\Books\\LOGLIN3\\BAYES\\")
load("post-samp.Rda")
rowMeans(beta)
apply(beta, 1, mean) # Redundant command
sdv=c(sd(beta[1,]),sd(beta[2,]))
sdv
#apply(beta, 1, sd) # Redundant command
apply(beta, 1, quantile,
    probs=c(0.025,.25,0.5,.75,0.975))