# NOT RUN {
testBF.prior(100)
library(MASS)
data(Pima.tr)
# use g = n
bas.glm(type ~ ., data=Pima.tr, family=binomial(),
betaprior=testBF.prior(nrow(Pima.tr)),
modelprior=uniform(), method="BAS")
# }
Run the code above in your browser using DataLab