[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