[R] Error in bugs.run -- R2WinBUGS

Cleber Nogueira Borges klebyn at yahoo.com.br
Wed Jun 18 23:17:18 CEST 2008


Hi,

I tried to use MethComp library and this library make use of the WinBUGS 
by R2WinBuGUS,
but I get the follow error in bugs.run:

*Error in bugs.run(n.burnin, bugs.directory, WINE = WINE, useWINE = 
useWINE,  : *
  Look at the log file and
try again with 'debug=TRUE' to figure out what went wrong within Bugs.

Anyone can help-me, please?
Thanks

Cleber



library( MethComp )
library( R2WinBUGS )
library( Epi )
# source("c:/stat/r/bxc/library.sources/methcomp/r/tab.repl.r")
# source("c:/stat/r/bxc/library.sources/methcomp/r/MCx.r")
# The oximetry example
#-------------------------------------------------------------------------------
# First the full story
data( ox )
str( ox )
tab.repl( ox )
*options( bugs.directory="c:/R/WinBUGS14/")*
( ox.mi.ir <- MethComp( ox, n.iter=5000 ) )
Comparison of 2 methods, using 354 measurements
on 61 items, with up to 3 replicate measurements,
(with the values in the set: 1 2 3 )
( 2 * 61 * 3 = 366 ):

No. items with measurements on each method:

Simulation run of a model with
- method by item and item by replicate interaction:
- using 4 chains run for 5000 iterations
  (of which 2500 are burn-in),
- monitoring every 3 values of the chain:
- giving a posterior sample of 3333 observations.

*Error in bugs.run(n.burnin, bugs.directory, WINE = WINE, useWINE = 
useWINE,  : *
  Look at the log file and
try again with 'debug=TRUE' to figure out what went wrong within Bugs.
In addition: Warning message:
In pt(q, df, lower.tail, log.p) : NaNs produced
 >


 > sessionInfo()
R version 2.7.0 (2008-04-22)
i386-pc-mingw32

locale:
LC_COLLATE=Portuguese_Brazil.1252;LC_CTYPE=Portuguese_Brazil.1252;LC_MONETARY=Portuguese_Brazil.1252;LC_NUMERIC=C;LC_TIME=Portuguese_Brazil.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base    

other attached packages:
[1] Epi_1.0.8       R2WinBUGS_2.1-8 coda_0.13-2     lattice_0.17-8  
MethComp_0.2.0  nlme_3.1-89   

loaded via a namespace (and not attached):
[1] grid_2.7.0  tools_2.7.0
 >









		
_______________________________________________________ 

Experimente já e veja as novidades.



More information about the R-help mailing list