model { for (i in 1:1334) { LOGCFU[i] ~ dt(MEAN[i], INVSIGSQ[KVTRTN[i]], V[KVTRTN[i]]) LOGCFUC[i] ~ dunif(0, 1) X[i] <- density(LOGCFU[i], DATA[i]) } for (i in 1335:1360) { LOGCFUC[i] ~ dt(MEAN[i], INVSIGSQ[KVTRTN[i]], V[KVTRTN[i]])C(, C[i]) X[i] <- cumulative(LOGCFUC[i], C[i]) } for (i in 1:1360) { MEAN[i] <- SMU[NSUBJID[i], 1] - SMU[NSUBJID[i], 2]*TIME[i] - SMU[NSUBJID[i], 3]*SGAMMA[NSUBJID[i]]* log((exp((TIME[i] - SKAPPA[NSUBJID[i]])/SGAMMA[NSUBJID[i]]) + exp(-(TIME[i] - SKAPPA[NSUBJID[i]])/ SGAMMA[NSUBJID[i]]))/(exp((SKAPPA[NSUBJID[i]])/SGAMMA[NSUBJID[i]]) + exp(-(SKAPPA[NSUBJID[i]])/ SGAMMA[NSUBJID[i]]))) PPO[i] <- X[i] } for (i in 1:85) { SALPHA[i] <- SMU[i, 1] SBETA1[i] <- SMU[i, 2] SBETA2[i] <- SMU[i, 3] SLAMBDA1[i] <- SBETA1[i] - SBETA2[i] SLAMBDA2[i] <- SBETA1[i] + SBETA2[i] SMU[i, 1:3] ~ dmnorm(MMU[OKVTRTN[1 + 16*(i - 1)], 1:3], MOMGINV[OKVTRTN[1 + 16*(i - 1)], 1:3, 1:3]) SKAPPA[i] ~ dnorm(MKAPPA[OKVTRTN[1 + 16*(i - 1)]], KAPINVSQ[OKVTRTN[1 + 16*(i - 1)]])T(2, 11) SGAMMA[i] ~ dnorm(MGAMMA[OKVTRTN[1 + 16*(i - 1)]], GAMINVSQ[OKVTRTN[1 + 16*(i - 1)]])T(0.1, 2) SEBA002[i] <- SBETA1[i] + SBETA2[i]*SGAMMA[i]*(log(exp((2 - SKAPPA[i])/SGAMMA[i]) + exp(-(2 - SKAPPA[i])/SGAMMA[i])) - log(exp((0 - SKAPPA[i])/SGAMMA[i]) + exp(-(0 - SKAPPA[i])/SGAMMA[i])))/2 SEBA007[i] <- SBETA1[i] + SBETA2[i]*SGAMMA[i]*(log(exp((7 - SKAPPA[i])/SGAMMA[i]) + exp(-(7 - SKAPPA[i])/SGAMMA[i])) - log(exp((0 - SKAPPA[i])/SGAMMA[i]) + exp(-(0 - SKAPPA[i])/SGAMMA[i])))/7 SEBA014[i] <- SBETA1[i] + SBETA2[i]*SGAMMA[i]*(log(exp((14 - SKAPPA[i])/SGAMMA[i]) + exp(-(14 - SKAPPA[i])/SGAMMA[i])) - log(exp((0 - SKAPPA[i])/SGAMMA[i]) + exp(-(0 - SKAPPA[i])/SGAMMA[i])))/14 SEBA214[i] <- SBETA1[i] + SBETA2[i]*SGAMMA[i]*(log(exp((14 - SKAPPA[i])/SGAMMA[i]) + exp(-(14 - SKAPPA[i])/SGAMMA[i])) - log(exp((2 - SKAPPA[i])/SGAMMA[i]) + exp(-(2 - SKAPPA[i])/SGAMMA[i])))/12 SEBA714[i] <- SBETA1[i] + SBETA2[i]*SGAMMA[i]*(log(exp((14 - SKAPPA[i])/SGAMMA[i]) + exp(-(14 - SKAPPA[i])/SGAMMA[i])) - log(exp((7 - SKAPPA[i])/SGAMMA[i]) + exp(-(7 - SKAPPA[i])/SGAMMA[i])))/7 } for (i in 1:6) { V[i] ~ dunif(2, 100) MALPHA[i] <- MMU[i, 1] MBETA1[i] <- MMU[i, 2] MBETA2[i] <- MMU[i, 3] MLAMBDA1[i] <- MBETA1[i] - MBETA2[i] MLAMBDA2[i] <- MBETA1[i] + MBETA2[i] MMU[i, 1:3] ~ dmnorm(D[1:3], IDENX[1:3, 1:3]) MKAPPA[i] ~ dunif(2, 11) MGAMMA[i] ~ dunif(0.1, 2) INVSIGSQ[i] ~ dgamma(0.0001, 0.0001) SIGSQ[i] <- 1/INVSIGSQ[i] KAPSIGSQ[i] ~ dunif(0.01, 30) KAPINVSQ[i] <- 1/KAPSIGSQ[i] GAMSIGSQ[i] ~ dunif(0.01, 5) GAMINVSQ[i] <- 1/GAMSIGSQ[i] SMUTILDA[i, 1:3] ~ dmnorm(MMU[i, 1:3], MOMGINV[i, 1:3, 1:3]) SALPHATILDA[i] <- SMUTILDA[i, 1] SBETA1TILDA[i] <- SMUTILDA[i, 2] SBETA2TILDA[i] <- SMUTILDA[i, 3] MEBA002[i] <- MBETA1[i] + MBETA2[i]*MGAMMA[i]*(log(exp((2 - MKAPPA[i])/MGAMMA[i]) + exp(-(2 - MKAPPA[i])/MGAMMA[i])) - log(exp((0 - MKAPPA[i])/MGAMMA[i]) + exp(-(0 - MKAPPA[i])/MGAMMA[i])))/2 MEBA007[i] <- MBETA1[i] + MBETA2[i]*MGAMMA[i]*(log(exp((7 - MKAPPA[i])/MGAMMA[i]) + exp(-(7 - MKAPPA[i])/MGAMMA[i])) - log(exp((0 - MKAPPA[i])/MGAMMA[i]) + exp(-(0 - MKAPPA[i])/MGAMMA[i])))/7 MEBA014[i] <- MBETA1[i] + MBETA2[i]*MGAMMA[i]*(log(exp((14 - MKAPPA[i])/MGAMMA[i]) + exp(-(14 - MKAPPA[i])/MGAMMA[i])) - log(exp((0 - MKAPPA[i])/MGAMMA[i]) + exp(-(0 - MKAPPA[i])/MGAMMA[i])))/14 MEBA214[i] <- MBETA1[i] + MBETA2[i]*MGAMMA[i]*(log(exp((14 - MKAPPA[i])/MGAMMA[i]) + exp(-(14 - MKAPPA[i])/MGAMMA[i])) - log(exp((2 - MKAPPA[i])/MGAMMA[i]) + exp(-(2 - MKAPPA[i])/MGAMMA[i])))/12 MEBA714[i] <- MBETA1[i] + MBETA2[i]*MGAMMA[i]*(log(exp((14 - MKAPPA[i])/MGAMMA[i]) + exp(-(14 - MKAPPA[i])/MGAMMA[i])) - log(exp((7 - MKAPPA[i])/MGAMMA[i]) + exp(-(7 - MKAPPA[i])/MGAMMA[i])))/7 for (j in 0:14) { MPLOT[i, j + 1] <- MMU[i, 1] - MMU[i, 2]*j - MMU[i, 3]*MGAMMA[i]*log((exp((j - MKAPPA[i])/MGAMMA[i]) + exp(-(j - MKAPPA[i])/MGAMMA[i]))/ (exp((MKAPPA[i])/MGAMMA[i]) + exp(-(MKAPPA[i])/MGAMMA[i]))) } } for (i in 1:6) { MEBAD002[i] <- MEBA002[i] - MEBA002[6] MEBAD007[i] <- MEBA007[i] - MEBA007[6] MEBAD014[i] <- MEBA014[i] - MEBA014[6] MEBAD214[i] <- MEBA214[i] - MEBA214[6] MEBAD714[i] <- MEBA714[i] - MEBA714[6] } D[1] <- 0 D[2] <- 0 D[3] <- 0 IDENX[1, 1] <- 0.0001 IDENX[1, 2] <- 0 IDENX[1, 3] <- 0 IDENX[2, 1] <- IDENX[1, 2] IDENX[2, 2] <- 0.0001 IDENX[2, 3] <- 0 IDENX[3, 1] <- IDENX[1, 3] IDENX[3, 2] <- IDENX[2, 3] IDENX[3, 3] <- 0.0001 for (i in 1:6) { MOMGINV[i, 1:3, 1:3] ~ dwish(IDEN[i, 1:3, 1:3], 3) MOMEGA[i, 1:3, 1:3] <- inverse(MOMGINV[i, 1:3, 1:3]) ALPSIGSQ[i] <- MOMEGA[i, 1, 1] BT1SIGSQ[i] <- MOMEGA[i, 2, 2] BT2SIGSQ[i] <- MOMEGA[i, 3, 3] ALPBT1SIGSQ[i] <- MOMEGA[i, 1, 2] ALPBT2SIGSQ[i] <- MOMEGA[i, 1, 3] BT1BT2SIGSQ[i] <- MOMEGA[i, 2, 3] } IDEN[1, 1, 1] <- 1.9845083337 IDEN[1, 1, 2] <- 0.1424260547 IDEN[1, 1, 3] <- -0.255846266 IDEN[1, 2, 1] <- IDEN[1, 1, 2] IDEN[1, 2, 2] <- 0.0216738387 IDEN[1, 2, 3] <- -0.001714245 IDEN[1, 3, 1] <- IDEN[1, 1, 3] IDEN[1, 3, 2] <- IDEN[1, 2, 3] IDEN[1, 3, 3] <- 0.1007109971 IDEN[2, 1, 1] <- 3.2813361357 IDEN[2, 1, 2] <- 0.23404549 IDEN[2, 1, 3] <- -0.42176524 IDEN[2, 2, 1] <- IDEN[2, 1, 2] IDEN[2, 2, 2] <- 0.0360726551 IDEN[2, 2, 3] <- -0.001910575 IDEN[2, 3, 1] <- IDEN[2, 1, 3] IDEN[2, 3, 2] <- IDEN[2, 2, 3] IDEN[2, 3, 3] <- 0.161357375 IDEN[3, 1, 1] <- 1.9485281239 IDEN[3, 1, 2] <- 0.1391632094 IDEN[3, 1, 3] <- -0.247040917 IDEN[3, 2, 1] <- IDEN[3, 1, 2] IDEN[3, 2, 2] <- 0.0213492645 IDEN[3, 2, 3] <- -0.000989594 IDEN[3, 3, 1] <- IDEN[3, 1, 3] IDEN[3, 3, 2] <- IDEN[3, 2, 3] IDEN[3, 3, 3] <- 0.0942288949 IDEN[4, 1, 1] <- 2.6021610817 IDEN[4, 1, 2] <- 0.1865544366 IDEN[4, 1, 3] <- -0.337417074 IDEN[4, 2, 1] <- IDEN[4, 1, 2] IDEN[4, 2, 2] <- 0.0291376766 IDEN[4, 2, 3] <- -0.001787706 IDEN[4, 3, 1] <- IDEN[4, 1, 3] IDEN[4, 3, 2] <- IDEN[4, 2, 3] IDEN[4, 3, 3] <- 0.1312407105 IDEN[5, 1, 1] <- 7.9831351344 IDEN[5, 1, 2] <- 0.5673254998 IDEN[5, 1, 3] <- -1.050960876 IDEN[5, 2, 1] <- IDEN[5, 1, 2] IDEN[5, 2, 2] <- 0.0944775205 IDEN[5, 2, 3] <- 0.0019173075 IDEN[5, 3, 1] <- IDEN[5, 1, 3] IDEN[5, 3, 2] <- IDEN[5, 2, 3] IDEN[5, 3, 3] <- 0.4189948514 IDEN[6, 1, 1] <- 3.2598540528 IDEN[6, 1, 2] <- 0.2344166633 IDEN[6, 1, 3] <- -0.41072744 IDEN[6, 2, 1] <- IDEN[6, 1, 2] IDEN[6, 2, 2] <- 0.0354100862 IDEN[6, 2, 3] <- -0.002830085 IDEN[6, 3, 1] <- IDEN[6, 1, 3] IDEN[6, 3, 2] <- IDEN[6, 2, 3] IDEN[6, 3, 3] <- 0.1559887863 }