data(Hald)
hald.bic <- bas.lm(Y ~ ., data=Hald, prior="BIC",
initprobs= "eplogp")
coef <- list2matrix.bas(hald.bic, "mle") # extract all coefficients
se <- list2matrix.bas(hald.bic, "mle.se")
models <- list2matrix.which(hald.bic) #matrix of model indicators
models <- which.matrix(hald.bic$which, hald.bic$n.vars) #matrix of model indicators
Run the code above in your browser using DataLab