# for Windows first install the BAS package by
# 1) download windows zip file from http://www.stat.duke.edu/~clyde/BAS/
# 2) In R, use the Package menu to install from a zip file on your computer

# library(BAS, lib.loc="~/Rlib64")

library(BAS)
library(MASS)
data(UScrime)  # loads in US Crime data
summary(UScrime)

UScrime[,-2] = log(UScrime[,-2])
help(BAS)
help(bas.lm)
crime.blm =  bas.lm(y ~ ., data=UScrime, n.models=2^15,
                     prior="g-prior",alpha=length(UScrime[,2]),
                     initprobs= "eplogp") 

summary(crime.blm)
plot(crime.blm)
image(crime.blm, subset=-1)
image(crime.blm, subset=-1, top=2^15)
image(crime.blm, subset=-1, top=50)
# takes a while to run:
 crime.coef = coefficients(crime.blm)
 crime.coef
 par(mfrow=c(3,2))
 plot(crime.coef, ask=FALSE)