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 [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