I am trying to estimate the following model that is based on the eyes
mixture model in the winbugs manual. All the independent variables are
standardized with x2 to x6 being factor scores. The model compiles but
I get a trap66 error. I have cut and paste the trap message after the model.
I have not been able to get much insight from the past
answers to this problem posted in the winbugs archives. I would appreciate any feedback in making this message go away.
Thanks for all your help,
Raghu
_______________________________________________
The model is as below (I have left out the data part):
model
{
for( i in 1 : N ) {
y[i] ~ dbin(p[i],100)
logit(p[i]) <- beta[i]
beta[i] ~ dnorm(mu[i], tau)
x1[i] <- (x0[i] - mean(x0[])) / sd(x0[])
mu[i] <- alpha1[T[i]] *x1[i] + alpha2[T[i]] * x2[i] + alpha3[T[i]] * x3[i] +
alpha4[T[i]] * x4[i] + alpha5[T[i]] * x5[i]
+ alpha6[T[i]] * x6[i] + alpha7[T[i]]
T[i] ~ dcat(P[])
}
P[1:2] ~ ddirch(alpha[])
theta1 ~ dnorm(0.0,1.0E-6)I(0,)
theta2 ~ dnorm(0.0,1.0E-6)I(0,)
theta3 ~ dnorm(
0.0,1.0E-6)I(0,)
theta4 ~ dnorm(0.0,1.0E-6)I(0,)
theta5 ~ dnorm(0.0,1.0E-6)I(0,)
theta6 ~ dnorm(0.0,1.0E-6)I(0,)
theta7 ~ dnorm(0.0,1.0E-6)I(0,)
alpha1[1] ~ dnorm(
0.0,1.0E-6)
alpha1[2] <- alpha1[1] + theta1
alpha2[1] ~ dnorm(0.0,1.0E-6)
alpha2[2] <- alpha2[1] + theta2
alpha3[1] ~ dnorm(0.0,1.0E-6)
alpha3[2] <- alpha3[1] + theta3
alpha4[1] ~ dnorm(0.0,1.0E-6)
alpha4[2] <- alpha4[1] + theta4
alpha5[1] ~ dnorm(0.0,1.0E-6)
alpha5[2] <- alpha5[1] + theta5
alpha6[1] ~ dnorm(
0.0,1.0E-6)
alpha6[2] <- alpha6[1] + theta6
alpha7[1] ~ dnorm(0.0,1.0E-6)
alpha7[2] <- alpha7[1] + theta7
tau ~ dgamma(0.001,0.01)
sigma <- 1 / sqrt(tau)
}
# Initial Values
list(alpha1=c(-3,NA), alpha2=c(1,NA), alpha3=c(5,NA), alpha4=c(-2,NA), alpha5=c(4,NA), alpha6=c(-3,NA), alpha7=c(0, NA),
theta1=2, theta2=5, theta3=7, theta4=1, theta5=4, theta6=2, theta7=3, tau=.1,
T = c(1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 2, 2,
1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 2, 1, 2, 1,
1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 2, 2, 1, 1, 1, 1, 2,
1, 1, 1, 2, 2, 1, 2, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1,
1, 1, 2, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 2, 2, 2,
2, 2, 1, 1, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 2, 1,
1, 1, 2, 2, 1, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 1,
1, 2, 2, 1, 2, 2, 2, 1, 2, 2, 1, 2, 1, 2, 1, 1, 1,
1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 2, 2, 2, 1, 2, 1, 2,
1, 2, 2, 2, 2, 2, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 2,
2, 2, 2, 2, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 2, 2,
2, 1, 2, 1, 1, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1,
1, 1, 2, 1, 1, 2, 2, 1, 2, 2, 1, 2, 1, 1, 1, 2, 2,
2, 2, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 1, 1,
1, 2, 2, 2, 2, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 2 ))
___________Trap Message_______________
TRAP 66 (postcondition violated)
UpdaterRejection.Updater.MCMC [000005B5H]
.deriv REAL -0.4106504744994321
.e REAL 0.0
.f
UpdaterRejection.Factory [101EA480H]
.i INTEGER 1
.k INTEGER 2
.lambdaL REAL 0.3619576527234962
.lambdaR REAL 2.435160950974093
.left REAL 3.525264210791281
.leftStar REAL 1.526054617532264
.linear REAL -1.505024195595193
.logFleft REAL -2.721441549298178
.logFmode REAL 39.06067862813028
.logFright REAL 34.56539095090484
.logQ REAL 1.0
.mode REAL 30.44507053697182
.mu REAL -1.505024195595193
.oldValue REAL 30.66604750192777
.overRelax BOOLEAN FALSE
.pL REAL 0.07248862937658747
.pM REAL 0.4897631891870716
.pR REAL 0.4377481814363409
.prec REAL -0.08431690610713367
.prior GraphStochastic.Node [102AC050H]
.rand REAL 4.67240400587114
.res INTEGER 0
.right REAL 35.31539243863755
.rightStar REAL -6.076427113308684
.s REAL 4.993302478421561
.sigma REAL 4.870321901665728
.tau REAL
2.0
.u REAL 0.4226899135032156
.updater UpdaterRejection.Updater [10618A40H]
.v REAL -0.2120758516044264
.values POINTER NIL
.x REAL 4.67240400587114
.y REAL
0.1493762485447229
BugsUpdaters.Updater.MCMC [000005B8H]
.chain INTEGER 0
.depth INTEGER 3
.i INTEGER 241
.name BugsNames.Name [10121F20H]
.ok BOOLEAN TRUE
.overRelax BOOLEAN FALSE
.res INTEGER 0
.size INTEGER 253
.string ARRAY 120 OF CHAR "" ...
.updater BugsUpdaters.Updater [10604930H]
BugsSampler.MCMC
[0000036DH]
.chain INTEGER 0
.cursor BugsSampler.List [106199C0H]
.depth INTEGER 3
.depth0 INTEGER 1
.depth1 INTEGER 4
.offset INTEGER 2287500
.ok BOOLEAN TRUE
.overRelax BOOLEAN FALSE
.updater BugsUpdaters.Updater [10604930H]
BugsCmds.Action.Step [00000B1FH]
.a BugsCmds.Action [101DE4D0H]
.chain INTEGER 0
.i INTEGER 4
.j INTEGER 0
.numChains INTEGER 1
BugsCmds.Action.Do [00000C58H]
.a BugsCmds.Action [101DE4D0H]
.elapsedTime LONGINT 032A21F103290384H
.profileList ARRAY 2024 OF CHAR "" ...
.res INTEGER 22041183
.s ARRAY 120 OF CHAR 3X ...
Services.Exec [00000136H]
.a Services.Action [101DE4D0H]
.t POINTER [03290384H]
Services.IterateOverActions [000002F4H]
.p Services.Action [101DE4D0H]
.t POINTER NIL
.time LONGINT 204620548
Services.StdHook.Step [0000034DH]
.h Services.StdHook [1013E380H]
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 525432
Kernel.Try [00003A61H]
.a INTEGER 525432
.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 4243644
.s ARRAY 256 OF SHORTCHAR "lú"" ...
.w INTEGER 2147348480
.wParam INTEGER 1
.wnd INTEGER 525432
<system> (pc=77D48733H, fp=0022FB84H)
<system> (pc=77D48815H, fp=0022FBECH)
<system> (pc=77D489CCH, fp=0022FC4CH)
<system> (pc=77D496C6H, fp=0022FC5CH)
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
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail