Institut für Mathematik
A novel Bayesian approach for variable selection in linear regression models
A novel Bayesian approach to the problem of variable selection in multiple linear regression models is proposed. In particular, a hierarchical setting which allows for direct speciﬁcation of a priori beliefs about the number of nonzero regression coefﬁcients as well as a speciﬁcation of beliefs that given coefﬁcients are nonzero is presented. This is done by introducing a new prior for a random set which holds the indices of the predictors with nonzero regression coefﬁcients. To guarantee numerical stability, a g-prior with an additional ridge parameter is adopted for the unknown regression coefﬁcients. In order to simulate from the joint posterior distribution an intelligent random walk Metropolis-Hastings algorithm which is able to switch between different models is proposed. For the model transitions a novel proposal, which prefers to add a priori or empirically important predictors to the model and further tries to remove less important ones, is used. Testing the algorithm on real and simulated data illustrates that it performs at least on par and often even better than other well-established methods. Finally, it is proven that under some nominal assumptions, the presented approach is consistent in terms of model selection.
Senka Haznadar (senka [dot] haznadar [at] aau [dot] at)