This function coerces the list object of models to a matrix and fill in the zeros to facilitate other computations.
which.matrix(which, n.vars)
a 'bas' model object x$which
the total number of predictors, x$n.vars
a matrix representation of x$which
, with number of rows equal
to the length of which.models or total number of models and number of
columns x$n.vars
which.matrix
coerces
x$which
into a matrix.
Other as.matrix methods:
list2matrix.bas()
,
list2matrix.which()
# NOT RUN {
data(Hald)
Hald.bic <- bas.lm(Y ~ ., data=Hald, prior="BIC", initprobs="eplogp")
# matrix of model indicators
models <- which.matrix(Hald.bic$which, Hald.bic$n.vars)
# }
Run the code above in your browser using DataLab