On 4 April 2012 05:35, Jack Tanner <i...@hotmail.com> wrote: > (Using BRugs 0.7-5, R 2.14.2 32-bit on 64-bit Windows 7, OpenBUGS 3.2.1) > > 1. BRugs crashes R for me as follows. Sorry about the lack of detail; > please let > me know if / how to supply a more useful bug report on this issue. > > fit <- BRugsFit(...) > # BRugs and OpenBUGS runs fine, the parameter estimates are reasonable > # across 3 chains > samplesBgr("beta") # crash > samplesBgr("beta", plot=FALSE) # also crash > > Have you plotted your histories? I haven't used samplesBgr() much, so I don't know how stable it is (although I do know it's slow).
> 2. I'd like to run a bunch of models via BRugs, saving the state after > each one, > like so: > > for (m in models) { > fit <- BRugsFit(paste(m$modelname, ".bug"), ...) > save(fit, paste(m$modelname, ".RData)) > } > > After I load the saved fitted objects in a new R process, BRugs functions > like > samplesStats() no longer work. Is it possible to somehow load the saved > objects > and use samplesStats and other BRugs functions? > > samplesStats() calls internal OpenBUGS functions (not R functions), so that would mean saving the whole BUGS run (like externalise in OpenBUGS itself. From you code fit$Stats should give you the same as sampleStats('*'): if you want more use BRugsFit(..., coda=T) and work with the coda object it produces (check the documentation for BRugsFit and coda). Bob -- Bob O'Hara Biodiversity and Climate Research Centre Senckenberganlage 25 D-60325 Frankfurt am Main, Germany Tel: +49 69 798 40216 Mobile: +49 1515 888 5440 WWW: http://www.bik-f.de/root/index.php?page_id=219 Blog: http://blogs.nature.com/boboh Journal of Negative Results - EEB: www.jnr-eeb.org [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.