[R] error calling Winbugs using R2WinBugs to run a multi-level model
Uwe Ligges
ligges at statistik.tu-dortmund.de
Wed Jun 6 19:17:10 CEST 2012
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.isojunno at 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 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.
>>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.
More information about the R-help
mailing list