setwd("c:\\E-drive\\Books\\LOGLIN3\\BAYES\\") load("trauma_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.05,0.95))