Coerce a BAS list object into a matrix
Models, coefficients, and standard errors in objects of class 'bas' are represented as a list of lists to reduce storage by omitting the zero entries. These functions coerce the list object to a matrix and fill in the zeros to facilitate other computations.
1 2 3 4 5 6 7
a 'bas' object
name of bas list to coerce
a vector of indices use to extract a subset
the total number of predictors,
list2matrix.bas(x, which) is equivalent to
list2matrix.which(x), however, the latter uses sapply rather than
which.matrix both coerce
x$which into a matrix.
a matrix representation of
x$what, with number of rows equal
to the length of which.models or total number of models and number of
Merlise Clyde firstname.lastname@example.org
1 2 3 4 5 6 7 8 9 10
## Not run: library(MASS) data(UScrime) UScrime[,-2] = log(UScrime[,-2]) crime.bic = bas.lm(y ~ ., data=UScrime, n.models=2^15, prior="BIC", initprobs= "eplogp") coef = list2matrix.bas(crime.bic, "ols") # extract all ols coefficients se = list2matrix.bas(crime.bic, "ols.se") models = list2matrix.which(crime.bic) #matrix of model indicators models = which.matrix(crime.bic$which, crime.bic$n.vars) #matrix of model indicators ## End(Not run)
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.