Given output from quantile regression stochastic search variable selection, this function returns a table of the 'best' models together with their associated empirical posterior probability.
topmodels(qrssvs, nmodels = 5, abbreviate = FALSE, minlength = 3)
An object of class qrssvs
. Typically this will be the
gamma
component of the list returned by SSVSquantreg
.
The number of models to tabulate.
Logical: should the names of the predictors be abbreviated?
If abbreviate
is set to TRUE
, the minimum
length of the abbreviations.
A table with the models and their associated posterior probability. The models are arranged in descending order of probability.
# NOT RUN {
# }
# NOT RUN {
set.seed(1)
epsilon<-rnorm(100)
set.seed(2)
x<-matrix(rnorm(1000),100,10)
y<-x[,1]+x[,10]+epsilon
qrssvs<-SSVSquantreg(y~x)
topmodels(qrssvs$gamma)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab