# NOT RUN {
# }
# 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.ebg = update(crime.bic, newprior="EB-global")
crime.zs = update(crime.bic, newprior="ZS-null")
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab