Retrieving estimation results from "bayes [...] saving(xyz)" command

[Stefan Lang]

Dear statalisters,

I was running a number of bayesian regressions when my pc crashed.
Using the "saving" suffix, I have now a number of datasets, containing the results from my 50000 draws.

Is there any way I can retrieve the resulting estimates from those datasets?



Just to illustrate what I want - the structure of my do file is the following:

bayes, [...] saving(xyz, replace): [...]
estimates store abc

Usually, after all the estimating is finished, I can access my estimates for post-estimation by typing

estimates use abc

Obviously, this doesn't work now. So is there a way to access my estimates from the xyz - datasets?

I think I have researched the manual quite thoroughly, but probably I'm missing something.

Thanks a lot in advance!

Best,
Stefan

[Nikolay Balov (StataCorp)]

Hi Stefan,

You can try the following workaround. First, copy your simulation dataset to a safe place other than your current directory, rerun your model with small number of iterations, say 1, keeping all other options the same, then replace the newly created dataset xyz.dta with the saved one, created before the crash. At the end you should have an estimation results working with your saved simulation data.

cap mkdir temp
cp xyz.dta temp/xyz.dta, replace

bayes, ... mcmcs(1) saving(xyz, replace) : ...
estimates store abc

cp temp/xyz.dta xyz.dta, replace

bayesstats summary

[Stefan Lang]

Hej Nikolay,

at first: sorry for my delayed answer!
Thank you very much for your advice!
Your advice worked out pretty well!

Best,
Stefan