## Not run:
# library(MASS)
# data(UScrime)
# UScrime[,-2] = log(UScrime[,-2])
# crime.bic = bas.lm(y ~ ., data=UScrime, n.models=2^15, prior="BIC",initprobs= "eplogp")
# crime.zs = update(crime.bic, newprior="ZS-null")
# ## End(Not run)Run the code above in your browser using DataLab