[R] Errors when saving output from WinBUGS to R

[Uwe Ligges]

On 03.10.2012 17:51, yiyiwang wrote:
> Dear all
>
> I used R2WinBUGS package's bugs() function to generate MCMC results. Then I
> tried to save the simulation draws in R, using read.bugs() function. Here is
> a simple test:
> ######################
> library(coda)
> library(R2WinBUGS)
> #fake some data to test
> beta0=1
> beta1=1.5
> beta2=-1
> beta3=2
> N=200
> x1=rnorm(N, mean=0,sd=1)
> x2=rnorm(N, mean=0,sd=1)
> x3=rnorm(N, mean=0,sd=1)
> lambda2= exp(beta0+ beta1*x1+ beta2 * x2 + beta3 * x3)
> y=rep(NA,N)
> for (i in 1: N) {
> y[i]=rpois(1,lambda2[i])
> }
> #generate inputs required by winbugs
> data <- list(Y20=y,x1=x1, x2=x2, x3=x3, N=N)
> inits<- function(){
> list(beta0=0, beta1=0, beta2=0, beta3=0)}
> parameters <- c("beta0", "beta1", "beta2", "beta3")
> #call winbugs using bugs() function
> model <- bugs(data, inits, parameters,
> model.file="C:/Jinmodel/test_codes/model_poisson.txt", n.chains=1,
> n.iter=30, n.burnin=0,n.thin=1, codaPkg=TRUE, DIC=TRUE, debug=TRUE,
> bugs.directory="C:/Users/Yiyi/Desktop/winbugs14/WinBUGS14")
> ###########################################
>
> I can get pretty good traceplots & parameter estimate outputs in WinBUGS.
> However, when I switched back to the R Console (by pressing ESC key)

.... which is the explanation. Just close WinBUGS regularly rather than 
interrupting the R function that executes WinBUGS and waits for its results.

You may also want to update to a recent version of R.

Best,
Uwe Ligges


 > , a
> warning message occurred:
> Warning message:
> running command '"C:/Users/Yiyi/Desktop/winbugs14/WinBUGS14/WinBUGS14.exe"
> /par "C:/Users/Yiyi/AppData/Local/Temp/Rtmpgdvr4Z/script.txt"' had status 93
>
> I ignored this message  & typed in the following code to save the outputs
> into R:
> ###################
> myres <- read.bugs(model)
>   ####################
> An error message occurred: Error in is.vector(X) : object 'model' not found
>
> I've been trying to solve this mystery for several weeks but couldn't make
> any breakthough.
>
> Any suggestions are highly appreciated!
>
> Best,
>
> Yiyi
>
> p.s. 0) the program is running in a Windows 7 Home Premium, 64-bit operating
> system machine.
>               R version:   2.13.2, WinBUGS 14
>         1) some old posts suggest running the software as the administrator
> (to kill the warning message). Can anyone shed more light on how to log in
> as the administrator? I am the sole user of my own laptop, not using the
> pc's at my office.  So I doubt administrator is the issue.
>         2) The second time I run the bugs() function, no warning message of
> status 93 emerged.
>         3) perhaps the warning message & the error message are somehow
> related??
>
>
>
> --
> View this message in context: http://r.789695.n4.nabble.com/Errors-when-saving-output-from-WinBUGS-to-R-tp4644899.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at 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.
>