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