I've run into persistent problems with OpenBUGS crashing when using
BRugs .53 and .71, and am hoping someone has suggestions. There is
obviously something unusual going on in my environment, but I'm at a
loss as to where to begin to try to solve it.
In a nutshell, what happens is that, as soon as I call "modelCheck()"
in BRugs, it gets an error or crashes ... but only some of the time
(90%< p< 100%). Following are details:
1. OpenBUGS 3.0.3 + BRugs 0.531:
It works occasionally, but approximately 90% of the time, I get the
following error from modelCheck():
Error in handleRes(res[[3]]) :
An OpenBUGS module or procedure was called that did not exist.
The specific code seems not to matter, but here is an example (model
taken from the OpenBUGS tutorial):
modelString =
" model
{
for (i in 1:N) {
r[i] ~ dbin(p[i], n[i])
b[i] ~ dnorm(0, tau)
logit(p[i])<- alpha0 + alpha1 * x1[i] + alpha2 * x2[i]
+ alpha12 * x1[i] * x2[i] + b[i]
}
alpha0 ~ dnorm(0, 1.0E-6)
alpha1 ~ dnorm(0, 1.0E-6)
alpha2 ~ dnorm(0, 1.0E-6)
alpha12 ~ dnorm(0, 1.0E-6)
tau ~ dgamma(0.001, 0.001)
sigma<- 1 / sqrt(tau)
}
"
print(modelString)
writeLines(modelString,con="model3.txt")
modelCheck( "model3.txt" )
Which (usually) produces:
modelCheck( "model3.txt" )
Error in handleRes(res[[3]]) :
An OpenBUGS module or procedure was called that did not exist.
I've copied at the end of this message an example from a single R
session that shows how it may work sometimes but not others.
2. OpenBUGS 3.2.1 + BRugs 0.71:
As above, the model occasionally works, but mostly it crashes R on
the modelCheck() line with the error "R for Windows GUI front-end has
stopped working".
I've tried the following combinations to try to get it to work:
A. Win7 32-bit + R 2.13 + BRugs 0.531 from standard CRAN repository
(installed from R)
B. Win7 32-bit + R 2.13 + BRugs 0.71 + OpenBUGS 3.2.1 (package and
EXE from OpenBUGS site)
C. Win7 64-bit [different machine] + R 2.13 (32-bit) + BRugs 0.531
D. Win7 64-bit + R 2.13 (32-bit) + BRugs 0.71 + OpenBugs 3.2.1
E. Win7 32-bit + *R 2.12.2* + BRugs 0.531 from standard CRAN
repository (installed from R)
F. Win7 64-bit + R 2.12.2 (32-bit) + BRugs 0.531
G. Win7 64-bit + R 2.12.2 (32-bit) + BRugs 0.71 + OpenBugs 3.2.1
H. ... and various combinations of removing all R and openbugs
versions, clean installs, and trying again
In all cases, it occasionally works -- especially on first re-install
of the package -- but then generally gives the "handleRes" error with
BRugs 0.531 or crashes with BRugs 0.71.
Thanks in advance for any suggestions!
=== example from a a single R session ===
R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i386-pc-mingw32/i386 (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
Loading required package: proto
Loading required package: grid
Loading required package: reshape
Loading required package: plyr
Attaching package: 'reshape'
The following object(s) are masked from 'package:plyr':
rename, round_any
[Previously saved workspace restored]
install.packages("BRugs")
Installing package(s) into 'C:/Users/cchap/Documents/R/win-library/2.13'
(as 'lib' is unspecified)
--- Please select a CRAN mirror for use in this session ---
also installing the dependency 'coda'
trying URL
'http://cran.fhcrc.org/bin/windows/contrib/2.13/coda_0.14-4.zip'
Content type 'application/zip' length 200622 bytes (195 Kb)
opened URL
downloaded 195 Kb
trying URL
'http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.13/BRugs_0.5-3.1.zip'
Content type 'application/zip' length 5613828 bytes (5.4 Mb)
opened URL
downloaded 5.4 Mb
package 'coda' successfully unpacked and MD5 sums checked
package 'BRugs' successfully unpacked and MD5 sums checked
The downloaded packages are in
C:\Users\cchap\AppData\Local\Temp\RtmpoNcfyf\downloaded_packages
library(BRugs)
Loading required package: coda
Loading required package: lattice
Welcome to BRugs running on OpenBUGS version 3.0.3
# Specify the model in BUGS language, but save it as a string in R:
modelString = "
+ model {
+ for ( i in 1:nFlips ) {
+ y[i] ~ dbern( theta )
+ }
+ theta ~ dunif( priorA , priorB )
+ priorA<- 0.6
+ priorB<- 1.0
+ }
+ " # close quote to end modelString
# Write the modelString to a file, using R commands:
writeLines(modelString,con="model3.txt")
# Use BRugs to send the model.txt file to BUGS, which checks the
model syntax:
modelCheck( "model3.txt" )
model is syntactically correct
[ . run the model successfully!! .]
[ . now change just the prior definition and try to re-run .]
modelString = "
+ model {
+ for ( i in 1:nFlips ) {
+ y[i] ~ dbern( theta )
+ }
+ theta ~ dunif( priorA , priorB )
+ priorA<- 0.4
+ priorB<- 0.6
+ }
+ " # close quote to end modelString
# Write the modelString to a file, using R commands:
writeLines(modelString,con="model3.txt")
# Use BRugs to send the model.txt file to BUGS, which checks the
model syntax:
modelCheck( "model3.txt" )
Error in handleRes(res[[3]]) :
An OpenBUGS module or procedure was called that did not exist.
[[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.
