Dear Bugs users,
I got a trap 66: postcondition violated error when trying to run a
binary logistic model. I searched the archive and found that many
people met this problem before. But I have not been able to get much
insight from the past answers to this problem posted in the winbugs
archives. Are there any persons who could tell me how to deal with this
situation? Is there anything wrong with my initial value? I would
appreciate any feedback in making this message go away.
Cheers
Tianji Cai
model {
for( i in 1 : 200 ) {
mu[i] ~ dnorm (0,tau)
for( j in 1 : 2 ) {
y[i,j] ~ dbin(k[i,j],1)
logit(k[i,j]) <-beta0 + mu[i] + inprod(beta[],x[,i,j])
}
}
beta0 ~ dnorm(0,1.0E-6)
for( p in 1 : 3 ) {
beta[p] ~ dnorm(0,1.0E-6)
}
tau~ dgamma(0.001,0.001)
sigma <- 1 / sqrt(tau)
}
list(beta0=0.01,beta1=0.01,beta2=0.01,beta3=0.01,mu=c(-0.06233,-0.1318,-0.04821,-0.07616,-0.06129,-0.04842,-0.0943,-0.05525,-0.1066,-0.05969,-0.06451,-0.06853,-0.04246,-0.04521,-0.08278,-0.07013,-0.1095,-0.09044,-0.06405,0.2209,-0.07029,-0.06236,0.2368,-0.05497,-0.08377,0.07456,-0.06142,-0.07017,-0.07087,-0.06462,0.5588,0.1259,-0.07882,0.2588,-0.05963,-0.06762,0.5444,-0.07938,-0.1219,-0.1197,-0.06593,0.2593,-0.0874,0.08601,0.2065,-0.07259,-0.06837,-0.09773,-0.06828,-0.06993,0.06575,-0.08374,0.2642,-0.06587,-0.04652,0.2622,0.2529,-0.05982,-0.08469,-0.07373,-0.1117,0.2379,-0.08004,-0.07063,-0.07908,-0.1358,-0.07947,-0.05022,-0.06236,-0.1264,-0.1245,-0.04786,0.2125,-0.09951,-0.03685,-0.067,-0.07344,-0.04183,0.1894,-0.05353,-0.1117,0.2662,-0.05344,-0.07508,-0.07351,-0.1082,0.2146,-0.1043,0.1819,-0.06674,-0.1453,-0.0942,-0.05416,-0.06924,-0.07701,-0.08942,0.2061,-0.07565,-0.05257,-0.1023,-0.09429,-0.08815,-0.1434,-0.06793,-0.06887,-0.09467,-0.06233,-0.05364,-0.1192,-0.1068,0.2523!
,-0.1209,0.2251,-0.1128,-0.05631,-0.08959,-0.05609,-0.07341,0.6015,-0.07288,-0.06327,-0.1234,-0.07751,-0.05912,0.5657,0.163,-0.08884,-0.1316,0.1979,-0.08886,-0.07585,0.2166,0.1971,0.2196,0.2354,0.2678,0.2139,0.2416,0.293,0.2637,-0.06403,0.2527,0.5038,-0.0981,-0.08673,-0.05482,0.1922,-0.08256,0.214,-0.1177,-0.05974,-0.06255,-0.09704,-0.09788,-0.08042,-0.09557,-0.03519,-0.1144,0.2783,-0.1604,-0.09812,-0.1071,0.2519,0.212,-0.06767,-0.06402,-0.1373,-0.08689,-0.07571,-0.09666,-0.03853,-0.05912,-0.0605,-0.07536,0.1738,-0.09463,-0.1269,-0.07748,-0.05664,0.2441,-0.06842,0.2259,-0.0686,-0.1007,0.5051,-0.1419,-0.09372,-0.08739,0.2117,-0.06861,-0.0751,-0.08019,0.2018,-0.09615,0.2194,-0.1353,-0.08311,-0.1564,0.259,-0.161), tau=
0.01)
TRAP 66 (postcondition violated)
UpdaterRejection.Updater.MCMC [000005B5H] .deriv REAL -3.710017026928494
.e REAL 0.0286304125979853
.f UpdaterRejection.Factory [01152170H] .i INTEGER 0
.k INTEGER 9999
.lambdaL REAL 0.5022615450579188
.lambdaR REAL 0.2695405419278885
.left REAL -5.20626732895774
.leftStar REAL 4.410488149153482
.linear REAL -3.987651160131523E+292
.logFleft REAL -10.56516899601791
.logFmode REAL -1772.983705814284
.logFright REAL -inf
.logQ REAL 5.369276064085929E-315
.mode REAL 4.783846633731022
.mu REAL 5.26517591373825E-315
.oldValue REAL 4.510830896808214
.overRelax BOOLEAN FALSE
.pL REAL 0.02516457633762698
.pM REAL 0.3845547811469755
.pR REAL 0.5902806425153976
.prec REAL -6.908626796291777
.prior GraphStochastic.Node [01102360H] .rand REAL 6.565894725501128
.res INTEGER 0
.right REAL 5.317451235720053
.rightStar REAL -inf
.s REAL 19.9590701754402
.sigma REAL 0.5336046019890319
.tau REAL -4.05041566895954E+50
.u REAL 0.6792391132001016
.updater UpdaterRejection.Updater [01151420H] .v REAL -0.3791545466815478
.values POINTER NIL
.x REAL 6.565894725501128
.y REAL 0.6615808222729623
BugsUpdaters.Updater.MCMC [000005B8H] .chain INTEGER 0
.depth INTEGER 1
.i INTEGER 0
.name BugsNames.Name [01015B60H] .ok BOOLEAN TRUE
.overRelax BOOLEAN FALSE
.res INTEGER 0
.size INTEGER 1
.string ARRAY 120 OF CHAR "" ...
.updater BugsUpdaters.Updater [01150A60H] BugsSampler.MCMC
[0000036DH] .chain INTEGER 0
.cursor BugsSampler.List [011628C0H] .depth INTEGER 1
.depth0 INTEGER 1
.depth1 INTEGER 3
.offset INTEGER 2287504
.ok BOOLEAN TRUE
.overRelax BOOLEAN FALSE
.updater BugsUpdaters.Updater [01150A60H] BugsCmds.Action.Step
[00000B1FH] .a BugsCmds.Action [011D72D0H] .chain INTEGER 0
.i INTEGER 48
.j INTEGER 0
.numChains INTEGER 1
BugsCmds.Action.Do [00000C58H] .a BugsCmds.Action [011D72D0H]
.elapsedTime LONGINT 628021F1627F0384H
.profileList ARRAY 2024 OF CHAR "" ...
.res INTEGER 1642156639
.s ARRAY 120 OF CHAR "" ...
Services.Exec [00000136H] .a Services.Action [011D72D0H]
.t POINTER [627F0384H]
Services.IterateOverActions [000002F4H]
.p Services.Action [011D72D0H] .t POINTER NIL
.time LONGINT 9992687
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 526808
Kernel.Try [00003A61H] .a INTEGER 526808
.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 14024704
.s ARRAY 256 OF SHORTCHAR "?", 9X, 0AX ...
.w INTEGER 4243644
.wParam INTEGER 1
.wnd INTEGER 526808
<system> (pc=77D13A67H, fp=0022FB88H)
<system> (pc=77D13B36H, fp=0022FBF0H)
<system> (pc=77D13D90H, fp=0022FC50H)
<system> (pc=77D1438BH, 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
|