I am trying to use Bayesian model averaging for variable selection with a large number of variables. In R, the BMS package allows to apply the method, with the option of using MCMC sampler (Metropolis Hastings algorithm) when the number of covariates is large.
Here is a sample code:
data(datafls) fls1 = bms(datafls, burn = 50000, iter=100000, g = "BRIC", mprior = "uniform", nmodel = 2000, mcmc="bd", user.int=F) result = coef(fls1)
However, if you run the same code twice, the results (i.e. the posterior probabilities) would be completely different. Does anyone know how to tune the code, so that the results are consistent for every run?