plots the posterior density function of coefficients in a BMA regression
pltdensbma_g(res)
* res = the results typed list of a Bayesian Model Averaging regression
* nothing (printed on a graphic window)
load(GROCERDIR+'/data/crime.dat') // take log of data listn = ['crime';'m';'ed';'po1';'po2';'lf';'mf';'pop';'nw';'u1';'u2';'gdp';'ineq';'prob';'time'] ; for i =1:size(listn,1) execstr('l'+listn(i)+' = log('+listn(i)+')') end // BMA estimation without printing the results rbma = bma_g('lcrime',30000,'lm','so','led','lpo1','lpo2','llf','lmf','lpop','lnw','lu1','lu2',... 'lgdp','lineq','lprob','ltime','burnin=10000','mcmc=''mc3''','noprint') // now plot the results pltdensbma_g(rbma) | ![]() | ![]() |