Dear Bugs Users
I'm a new WinBugs user and I'm working with binary data.
Firstly I´ve took the example of WinBugs about “Beetles: choice of link
function” and the correspondent model program is:
model
{
for( i in 1 : N ) {
r[i] ~ dbin(p[i],n[i])
logit(p[i]) <- alpha.star + beta * (x[i] - mean(x[]))
#p[i]<- phi(alpha.star + beta * (x[i] - mean(x[])))
#cloglog(p[i]) <- alpha.star + beta * (x[i] - mean(x[]))
rhat[i] <- n[i] * p[i]
}
alpha <- alpha.star - beta * mean(x[])
beta ~ dnorm(0.0,0.001)
alpha.star ~ dnorm(0.0,0.001)
}
I´ve run these program using Bliss (1935) data and every thing was ok.
I need to run these same model program to another set data.
When I used logit link function it ran ok, but it doesn´t work as well as I
wish to probit link function.
The error message is following bellow. Please could you help me with these
question?
Best regards.
Simone
undefined real result
Math.Ln [0000018BH]
.x REAL -inf
UpdaterDFreeARS.SampleHull [00001198H]
.d REAL inf
.h UpdaterDFreeARS.Hull
.lBounded BOOLEAN FALSE
.uBounded BOOLEAN FALSE
.n INTEGER 23
.uMax REAL 4597705.059932021
.yMax REAL -1306.553150852046
.lowerBound REAL -inf
.upperBound REAL inf
.x ARRAY 51 OF REAL Elements
.y ARRAY 51 OF REAL Elements
.z ARRAY 101 OF REAL Elements
.U ARRAY 101 OF REAL Elements
.expU ARRAY 101 OF REAL Elements
.dz ARRAY 101 OF REAL Elements
.I ARRAY 101 OF REAL Elements
.i INTEGER 1
.lower REAL -0.003608119902629419
.rand REAL 2.289103025445015E-7
.sample REAL 2.708217690480629E-315
.squeeze BOOLEAN Undefined191
.upper REAL 8.15781491100198E-69
UpdaterDFreeARS.Ars [0000173FH]
.calls INTEGER 499
.h UpdaterDFreeARS.Hull Fields
.logRand REAL -0.004850367932147044
.lower REAL -0.003608119902629419
.node GraphStochastic.Node NIL
.res INTEGER 0
.sample REAL 2.708217690480629E-315
.squeeze BOOLEAN Undefined191
.updater UpdaterDFreeARS.Updater [010E7A90H]
.upper REAL 8.15781491100198E-69
UpdaterDFreeARS.Updater.MCMC [00002F00H]
.h UpdaterDFreeARS.Hull Fields
.i INTEGER 835468831
.k INTEGER 835387952
.oldValue REAL 0.0
.overRelax BOOLEAN FALSE
.prior GraphStochastic.Node [010DBAB0H]
.rand REAL -1.204323445990148
.res INTEGER 0
.updater UpdaterDFreeARS.Updater [010E7A90H]
.values POINTER NIL
BugsUpdaters.Updater.MCMC [000005B8H]
.chain INTEGER 0
.depth INTEGER 1
.i INTEGER 0
.name BugsNames.Name [010D78A0H]
.ok BOOLEAN TRUE
.overRelax BOOLEAN FALSE
.res INTEGER 0
.size INTEGER 1
.string ARRAY 120 OF CHAR "" ...
.updater BugsUpdaters.Updater [010E7A40H]
BugsSampler.MCMC [0000036DH]
.chain INTEGER 0
.cursor BugsSampler.List [010E7B40H]
.depth INTEGER 1
.depth0 INTEGER 1
.depth1 INTEGER 2
.offset INTEGER 2287504
.ok BOOLEAN TRUE
.overRelax BOOLEAN FALSE
.updater BugsUpdaters.Updater [010E7A40H]
BugsCmds.Action.Step [00000B29H]
.a BugsCmds.Action [01071620H]
.chain INTEGER 0
.i INTEGER 1
.j INTEGER 0
.numChains INTEGER 1
BugsCmds.Action.Do [00000C62H]
.a BugsCmds.Action [01071620H]
.elapsedTime LONGINT 331221FB33110384H
.profileList ARRAY 2024 OF CHAR "d" ...
.res INTEGER 2291454
.s ARRAY 120 OF CHAR 0A974X, 77D2X, 83B0X, 97X, 0A915X ...
Services.Exec [00000136H]
.a Services.Action [01071620H]
.t POINTER [33110384H]
Services.IterateOverActions [000002F4H]
.p Services.Action [01071620H]
.t POINTER NIL
.time LONGINT 1061093
Services.StdHook.Step [0000034DH]
.h Services.StdHook [0101E380H]
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 131572
Kernel.Try [00003A61H]
.a INTEGER 131572
.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 13369344
.s ARRAY 256 OF SHORTCHAR "pú"" ...
.w INTEGER 4243644
.wParam INTEGER 1
.wnd INTEGER 131572
<system> (pc=77D286CAH, fp=0022FB88H)
<system> (pc=77D2879EH, fp=0022FBF0H)
<system> (pc=77D28A30H, fp=0022FC50H)
<system> (pc=77D2910FH, fp=0022FCA0H)
Kernel.Start [00002B8CH]
.code PROCEDURE HostMenus.Loop
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail [log in to unmask]
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
|