On 06.06.2012 16:51, ilai wrote:
Untested because I don't have (use) winbugs and you didn't provide dat*.
But consider
a<- 4 ; f<- 6
list('a','f')
list(a,f)
list(a=a,f=f)
My guess is you wanted sp.data to be a named list, not a list of names...
That's also OK, from ?bugs.data:
data: either a named list (names corresponding to variable names in
the ‘model.file’) of the data for the ‘WinBUGS’ model, _or_
(which is not recommended and unsafe) a vector or list of the
names of the data objects used by the model
To really now what is going on, I'd need the model file and the data. My
suspicion is that the matrices have to be transposed.
Best,
Uwe Ligges
HTH
On Wed, Jun 6, 2012 at 4:12 AM, Saana Isojunno<
saana.isoju...@googlemail.com> wrote:
Dear all,
I'm calling Winbugs (1.4.3) through R2WinBugs (2.1-18 coda_0.14-7) to
fit a switching random walk model, but come up with an instant trap
with the log only displaying 'check('.
I will paste the trap with session info below; I'd be very grateful
for any ideas.
Couple of leads:
1. I presume the problem relates to the r package itself or the way I
call bugs(), because I can use the same text files specifying the
model and data directly in Winbugs and it runs fine (i.e syntax ok,
compilation ok, updates slow but no traps).
2. The problem occurs in r only when I try to fit the model to
multiple individuals, i.e. the data contains a matrix of step lengths
(rows) and individuals (columns) instead of a vector for just one
individual. I get the same error message regardless of the number of
data rows in each column (I even tried just one).
The model loops over the path of each animal, estimating a hidden
movement state and their parameters. For 4 individuals with 100 data
points each the data looks something like this:
dat1 : num 100
dat2 : int 4
dat3 : num [1:4] 8 4 2 5
dat4 : num [1:100, 1:4] 1 1 1 1 1 2 2 2 2 2 ...
dat5 : num [1:100, 1:4] 2 2 2 2 2 1 2 2 2 2 ...
dat6 : num [1:100, 1:4] 16 34.3 33.5 27.9 14.9 ...
dat7 : num [1:100, 1:4] 0.357 0.474 0.487 0.495 0.524 ...
dat8 : num [1:50, 1:4] 36.4 294.5 24.4 21.1 422.8 ...
This is how I've called WinBugs in r:
# write data to text file
sp.data = list("dat1","dat2","dat3","dat4","dat5","dat6","dat7","dat8")
bugs.data(sp.data, digits=5, data.file="dir1\\data1.txt")
# test the model runs
fit =
bugs(data=paste("C:\\Users\\User1\\Documents\\dir1\\data1.txt",dataFile,sep=""),
inits=NULL, parameters.to.save=list('par1','par2','par3'),
model.file=modelFile,
debug=TRUE, n.chains=3, n.iter=20, n.burnin=3, n.thin=1,
digits=4)
## The trap
incompatible copy
BugsScript.Action.Do [00000436H]
.a BugsScript.Action [025B6790H]
.argNum INTEGER 0
.bugsCommands ARRAY 240 OF CHAR 7877X, 75A5X, 0B17X, 3701X
...
.p ARRAY 3, 120 OF CHAR Elements
.s BugsScanners.Scanner Fields
.scriptCommand ARRAY 240 OF CHAR "#Bugs:check" ...
.vectorName BOOLEAN FALSE
Services.Exec [00000136H]
.a Services.Action [025B6790H]
.t POINTER [64E10170H]
Services.IterateOverActions [000002F4H]
.p Services.Action [025B6790H]
.t POINTER NIL
.time LONGINT 4375656
Services.StdHook.Step [0000034DH]
.h Services.StdHook [0248E380H]
HostWindows.Idle [00004A86H]
.focus BOOLEAN FALSE
.tick Controllers.TickMsg Fields
.w HostWindows.Window NIL
HostMenus.TimerTick [00003422H]
.lParam INTEGER 0
.ops Controllers.PollOpsMsg Fields
.wParam INTEGER 1
.wnd INTEGER 1311298
Kernel.Try [00003A61H]
.a INTEGER 1311298
.b INTEGER 1
.c INTEGER 0
.h PROCEDURE HostMenus.TimerTick
HostMenus.ApplWinHandler [00003841H]
.Proc PROCEDURE NIL
.hit BOOLEAN FALSE
.lParam INTEGER 0
.message INTEGER 275
.res INTEGER 1664639202
.s ARRAY 256 OF SHORTCHAR "" ...
.w INTEGER 1970768325
.wParam INTEGER 1
.wnd INTEGER 1311298
<system> (pc=75778816H, fp=0027FB38H)
<system> (pc=7577898DH, fp=0027FBB0H)
<system> (pc=75778AB8H, fp=0027FC14H)
<system> (pc=757790E2H, fp=0027FC24H)
HostMenus.Loop [00003BDEH]
.done BOOLEAN FALSE
.f SET {0..5}
.n INTEGER 0
.res INTEGER 0
.w HostWindows.Window NIL
Kernel.Start [00002B8CH]
.code PROCEDURE HostMenus.Loop
## my current R session specs:
R version 2.15.0 (2012-03-30)
Platform: x86_64-pc-mingw32/x64 (64-bit)
locale:
[1] LC_COLLATE=English_United Kingdom.1252
[2] LC_CTYPE=English_United Kingdom.1252
[3] LC_MONETARY=English_United Kingdom.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United Kingdom.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] R2WinBUGS_2.1-18 coda_0.14-7 lattice_0.20-6
loaded via a namespace (and not attached):
[1] grid_2.15.0 tools_2.15.0
## also tried these specs:
R version 2.13.2 (2011-09-30)
Platform: x86_64-pc-mingw32/x64 (64-bit)
locale:
[1] LC_COLLATE=English_United Kingdom.1252
[2] LC_CTYPE=English_United Kingdom.1252
[3] LC_MONETARY=English_United Kingdom.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United Kingdom.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] R2WinBUGS_2.1-18 coda_0.14-4 lattice_0.19-33
loaded via a namespace (and not attached):
[1] grid_2.13.2 tools_2.13.2
Thank you!
Saana Isojunno
PhD candidate
Scottish Oceans Institute
University of St Andrews
Scotland
______________________________________________
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.
[[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.
______________________________________________
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.
