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