Hello bugs
I am Graduate student in Biostatistics.
I run two code for estimate the parameters in survival models but BUGS shows this
"Sorry something went wrong in procedure stack. Value in module GraphStack"
codes are attached.
In FG code this error is in the process of updating and in FS code this error is in the process of load inits.
Kind Regards
Bonakchi.
[log in to unmask]
# FS model
model {
C<-10000
for(i in 1 : N){
log(theta1[i]) <- B1[1]*Age[i]+B1[2]*Gender[i]+B1[3]*Treatment.course[i]+B1[4]*WBC[i]+B1[5]*plt[i]+B1[6]*LDH[i]+B1[7]*The.cell.line[i]+B1[8]*CNS[i]+B1[9]*hemorraahagic[i]+B1[10]*mediastinal.mass[i]+B1[11]*reumatoid.signs[i]+B1[12]*tumor.lysis.synd[i]+B1[13]*hepatosplenomegalia[i]+B1[14]*Lymphadenopathy[i]+B1[15]*Radiotherapy[i]
log(theta2[i])<- B2[1]*Age[i]+B2[2]*Gender[i]+B2[3]*Treatment.course[i]+B2[4]*WBC[i]+B2[5]*plt[i]+B2[6]*LDH[i]+B2[7]*The.cell.line[i]+B2[8]*CNS[i]+B2[9]*hemorraahagic[i]+B2[10]*mediastinal.mass[i]+B2[11]*reumatoid.signs[i]+B2[12]*tumor.lysis.synd[i]+B2[13]*hepatosplenomegalia[i]+B2[14]*Lymphadenopathy[i]+B2[15]*Radiotherapy[i]
zeros [i] <- 0
zeros[i] ~ dpois (zero.mean[i])
zero.mean[i] <- -l[i]+C
l[i]<-equals(cause[i],1)*L1[i]+
equals(cause[i],2)*L2[i]+
equals(cause[i],0)*L0[i]
L1[i]<-log(theta1[i])+logh1time[i]-theta1[i]*H1time[i]
logh1time[i]<-log(lambda1[1])*equals(D1[i],1)+
log(lambda1[2])*equals(D1[i],2)+
log(lambda1[3])*equals(D1[i],3)+
log(lambda1[4])*equals(D1[i],4)
H1time[i]<- G1[i]+G2[i]
G1[i]<- lambda1[1]*time[i]*equals(D1[i],1)+
(lambda1[1]*141+lambda1[2]*(time[i]-141))*equals(D1[i],2)+
((lambda1[1]+lambda1[2])*141+lambda1[3]*(time[i]-282))*equals(D1[i],3)
G2[i]<- ((lambda1[1]+lambda1[2]+lambda1[3])*141+lambda1[4]*(1-exp(-time[i]+423)))*equals(D1[i],4)
L2[i]<- log(theta2[i])+logh2time[i]-theta2[i]*H2time[i]-((lambda1[1]+lambda1[2]+lambda1[3])*141+lambda1[4])*theta1[i]
logh2time[i]<-log(lambda2[1])*equals(D2[i],1)+
log(lambda2[2])*equals(D2[i],2)+
log(lambda2[3])*equals(D2[i],3)+
log(lambda2[4])*equals(D2[i],4)
H2time[i]<- F1[i]+F2[i]
F1[i]<- lambda2[1]*time[i]*equals(D2[i],1)+
(lambda2[1]*92+lambda2[2]*(time[i]-92))*equals(D2[i],2)+
((lambda2[1]+lambda2[2])*92+lambda2[3]*(time[i]-184))*equals(D2[i],3)
F2[i]<- ((lambda2[1]+lambda2[2]+lambda2[3])*92+lambda2[4]*(time[i]-276))*equals(D2[i],4)
L0[i]<-log(L00[i])
L00[i]<-exp(-H1time[i]*theta1[i])-(1-exp(-H2time[i]*theta2[i]))*(exp(-((lambda1[1]+lambda1[2]+lambda1[3])*141+lambda1[4])*theta1[i]))
}
for (j in 1:3)
{
lambda1[j] ~ dgamma(.01, .01)
}
lambda1[4] ~ dgamma(2,2)
for (r in 1:4)
{
lambda2[r] ~ dgamma(.01, .01)
}
for(q in 1:15)
{
B1[q] ~ dnorm(0.0, 0.0001)
B2[q] ~ dnorm(0.0, 0.0001)
}
}
#initial
list(B1=c(0,.1,.15,.17,.2,.24,.3,.33,.36,.4,.55,.6,.62,.7,.73), B2=c(.1,.14,.2,.23,.26,.3,.34,.41,.45,.46,.5,.55,.6,.67,.71), lambda1 = c(.1,.1,.1,.1),lambda2 = c(.1,.1,.1,.1))
#Data
list(N = 200,D1=c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3,
3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4),
D2=c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3,
3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,
3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,
3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,
4, 4, 4, 4, 4, 4, 4, 4, 4),
Age=c(11, 9, 11, 5, 13, 4, 11, 2, 9, 16, 11, 9, 2, 13, 2, 4, 3, 1,
7, 3, 6, 13, 2, 2, 2, 8, 0.37, 12, 1, 10, 13, 2, 2, 4, 13, 4,
2, 5, 2.7, 2.3, 8, 13, 9, 2, 3.5, 12.5, 4, 3, 7, 3, 11, 13, 7,
4, 5, 5, 1, 4, 10, 3, 11, 5, 0.4, 3, 4, 2, 6, 2.5, 6, 2, 12,
3, 12, 6, 3, 3, 2, 11, 13, 2.5, 6, 7, 3, 3, 9, 0.83, 4, 2, 7,
2, 11, 2, 5, 4, 3, 4, 2, 2, 14, 2, 1, 3, 9, 2, 1, 3, 2, 4, 8,
5, 2, 3, 5, 3, 4, 2, 6, 9, 2, 1, 9, 4, 3, 7, 4, 12, 11, 2, 3,
13, 2, 10, 5, 8, 11, 8, 2, 9, 9, 3, 3, 6, 5, 7, 4, 5, 3, 4, 2,
2, 3, 9, 5, 4, 14, 2, 8, 2, 9, 3, 5, 4, 9, 4, 5, 9, 1, 1, 3,
2, 14, 5, 4, 3, 7, 1, 8, 13, 3, 4, 10, 3, 7, 7, 4, 6, 12, 4,
4, 3, 5, 5, 3, 3, 9, 7, 5, 3, 14, 8),
Gender=c(1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1, 0, 1,
0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1,
1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0,
0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1,
0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1,
0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1,
0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0,
0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0,
1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1,
0, 1, 1, 0, 0, 0, 1, 0, 1),
Treatment.course=c(1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1,
1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0),
WBC=c(3, 1, 1, 1, 1, 2, 1, 1, 2, 1, 4, 4, 1, 1, 3,
1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 4, 2, 4, 2, 1,
1, 1, 4, 1, 3, 1, 1, 1, 2, 1, 4, 1, 2, 1, 3, 2,
1, 1, 1, 3, 1, 2, 1, 1, 1, 4, 1, 1, 1, 4, 2, 3,
1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 1,
2, 1, 4, 2, 4, 1, 2, 1, 3, 1, 2, 1, 2, 1, 1, 2,
1, 1, 1, 4, 4, 3, 2, 1, 2, 3, 2, 2, 1, 1, 1, 1,
1, 2, 2, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 0, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1,
1, 1, 3, 2, 1, 1, 2, 1, 1, 3, 1, 1, 1, 2, 2, 1,
4, 1, 2, 1, 1, 1, 2, 1, 1, 2, 1, 4, 1, 1, 3, 1,
2, 2, 1, 1, 1, 1, 1, 1, 4),
plt=c(21000, 19000, 82000, 44000, 298000, 49000, 38000, 55000,
15000, 56000, 26000, 350000, 12000, 297000, 1000, 90000,
50000, 3000, 10000, 39000, 10000, 30000, 72000, 109000,
389000, 17000, 94000, 31000, 27000, 175000, 148000, 57000,
102000, 76000, 300000, 29000, 18000, 11000, 17000, 46000,
11000, 65000, 402000, 23000, 40000, 178000, 63000, 20000,
58000, 86000, 40000, 183000, 19000, 167000, 10000, 28000,
216000, 211000, 32000, 548000, 13000, 4000, 1600, 11000,
25000, 17000, 7000, 22000, 45000, 93000, 481000, 97000,
76000, 47000, 3000, 303000, 42000, 347000, 99000, 52000,
41000, 45000, 173000, 17000, 156000, 415000, 17000, 70000,
80000, 45000, 422000, 330000, 79000, 423000, 7000, 400000,
3000, 31000, 16000, 16000, 58000, 17000, 123000, 41000,
16000, 8000, 116000, 9000, 140000, 15000, 36000, 36000,
40000, 29000, 38000, 3000, 12000, 70000, 200000, 28000,
6000, 52000, 31000, 35000, 8000, 93000, 136000, 45000,
81000, 0, 330000, 38000, 93000, 30000, 147000, 29000,
190000, 27000, 175000, 122000, 79000, 20000, 35000, 29900,
186000, 32000, 65000, 56000, 211000, 359000, 129000, 229000,
200000, 290000, 127000, 126000, 40000, 55000, 250000,
292000, 1000, 25000, 66000, 97000, 38000, 257000, 384000,
178000, 93000, 148000, 14000, 127000, 10000, 28000, 38000,
31000, 10000, 67000, 12000, 9000, 173000, 10000, 7000,
32000, 17000, 30000, 49000, 96000, 22000, 25000, 59000,
57000, 57000, 10000, 95000, 161000, 32000, 18000, 112000,
30000),
LDH=c(236, 0, 604, 869, 823, 1549, 0, 403, 0, 496, 744,
3194, 697, 0, 1213, 1030, 0, 649, 346, 0, 1596, 240,
0, 0, 0, 2151, 480, 3699, 2731, 2532, 7900, 30, 300,
5666, 1577, 867, 733, 900, 533, 600, 0, 12402, 422,
1487, 879, 4231, 1918, 225, 661, 214, 0, 716, 1054,
0, 0, 2372, 0, 0, 0, 0, 0, 964, 1961, 0, 1554, 427,
848, 0, 0, 872, 335, 807, 1304, 665, 508, 0, 1451,
0, 862, 1078, 0, 0, 351, 1341, 0, 824, 0, 0, 617,
1245, 1256, 0, 0, 748, 0, 362, 410, 539, 642, 0, 1038,
615, 647000, 627, 1560, 1104, 500, 660, 1859, 566, 0,
0, 1399, 439, 1145, 512, 961, 1329, 0, 972, 0, 1419,
0, 471, 488, 0, 648, 6500, 1264, 0, 620, 408, 1342,
427, 719, 0, 840, 472, 282, 744, 8740, 484, 622, 1758,
492, 1052, 0, 642, 1656, 0, 1036, 538, 450, 0, 462,
1209, 0, 707, 476, 766, 602, 1987, 605, 3764, 1182,
516, 178, 384, 0, 618, 0, 2061, 886, 1289, 647, 1725,
0, 4599, 922, 442, 418, 1726, 410, 859, 1102, 1650,
1748, 435, 723, 1988, 1516, 1112, 0, 638, 856, 0, 490,
1019, 638, 8330),
The.cell.line=c(1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2,
1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,
1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,
1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1,
2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1,
1, 1, 0, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1,
1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1),
CNS=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0),
hemorraahagic=c(0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0),
mediastinal.mass=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0),
reumatoid.signs=c(0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1,
0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1,
1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0,
0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1,
0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1,
0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1,
0, 1, 1, 1, 1, 1, 0, 0, 1),
tumor.lysis.synd=c(0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0),hepatosplenomegalia=c(0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1,
0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1,
0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0,
0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1,
1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0,
0, 1, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1,
1, 0, 1, 1, 0, 1, 0, 1, 1),
Lymphadenopathy=c(0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0,
1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1,
0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1,
1, 0, 0, 1, 0, 0, 0, 1, 0),
Radiotherapy=c(0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1,
1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1,
1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 1, 1,
1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1,
1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0,
0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1,
1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1,
0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1,
0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1,
0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 0, 1),
time=c(22, 24, 26, 31, 42, 44, 48, 74, 81, 92, 98, 98,
138, 138, 140, 141, 142, 146, 149, 171, 180, 200, 201,
219, 226, 234, 271, 279, 287, 287, 288, 373, 373, 381,
395, 397, 410, 420, 420, 422, 445, 457, 461, 462, 468,
472, 486, 486, 507, 511, 520, 520, 544, 566, 604, 615,
639, 647, 650, 683, 688, 705, 725, 749, 752, 765, 766,
769, 793, 793, 913, 915, 930, 930, 937, 944, 947, 961,
976, 994, 997, 1012, 1013, 1023, 1025, 1035, 1053, 1066,
1067, 1093, 1101, 1134, 1159, 1185, 1187, 1210, 1223,
1224, 1330, 1330, 1351, 1394, 1401, 1402, 1411, 1412,
1417, 1422, 1464, 1467, 1479, 1493, 1502, 1534, 1548,
1559, 1592, 1605, 1637, 1649, 1659, 1666, 1666, 1673,
1676, 1696, 1712, 1716, 1724, 1745, 1759, 1763, 1794,
1796, 1806, 1826, 1828, 1830, 1831, 1844, 1853, 1858,
1859, 1887, 1901, 1908, 1927, 1930, 1938, 1938, 1956,
1957, 1988, 2007, 2015, 2019, 2031, 2075, 2091, 2103,
2106, 2107, 2139, 2166, 2173, 2199, 2200, 2200, 2209,
2210, 2217, 2220, 2306, 2310, 2328, 2358, 2359, 2390,
2402, 2435, 2440, 2448, 2457, 2520, 2612, 2623, 2641,
2642, 2711, 2748, 2777, 2791, 2839, 2844, 2861, 2878,
2897, 2904, 2925, 2932),
cause=c(2, 0, 2, 0, 0, 0, 2, 0, 1, 2, 0, 1, 0, 2, 1,
0, 0, 2, 0, 2, 1, 0, 0, 2, 2, 1, 2, 0, 1, 1, 1,
2, 2, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 0, 2, 1, 0, 0, 0, 1, 2, 1, 0, 0,
0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0,
1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0)
)
........................................................................................................................................................................
# FG model
model
{
# Set up data
for(i in 1 : N) {
for(j in 1 : T) {
# risk set = 1 if obs.t >= t
Y[i, j] <- step(obs.t[i] - t[j] + eps)
# counting process jump = 1 if obs.t in [ t[j], t[j+1] )
# i.e. if t[j] <= obs.t < t[j+1]
dN[i, j] <- Y[i, j ] *step(t[j+1] - obs.t[i] - eps)*(1-equals(fail[i],0))+Y[i, j ] *step(obs.t[i] -t[j] +eps)*equals(fail[i],2)
}
}
# Model
for(j in 1 : T)
# T is the number of observed failure time
{
for(i in 1 : N)
{
dN[i, j] ~ dpois(Idt[i, j])
#1 for cause 1 and 2 for cause 2
Idt[i, j] <- equals(fail[i],1)*Y[i, j] *exp(b[1]*Age[i]+b[2]*Gender[i]+b[3]*Treatment.course[i]+b[4]*WBC[i]+
b[5]*plt[i]+b[6]*LDH[i]+b[7]*The.cell.line[i]+b[8]*CNS[i]+b[9]*hemorraahagic[i]+b[10]*mediastinal.mass[i]+b[11]*reumatoid.signs[i]+b[12]*tumor.lysis.synd[i]+b[13]*hepatosplenomegalia[i]+b[14]*Lymphadenopathy[i]+b[15]*Radiotherapy[i]) * dL0[j]
}
dL0[j] ~ dgamma(mu[j], c)
mu[j] <- dL0.star[j] * c # prior mean hazard
}
for(k in 1 : p)
{
b[k] ~ dnorm(0.0, 0.0001)
}
c <- 2 r <-1
for (j in 1 : T) {
dL0.star[j] <- r * (t[j+1]-t[j])
}
}
#INITS
list( b=c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1),
dL0 = c(2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,2,2,2,2,2,2)
)
#DATA
list(N = 200,T=35,eps=0.00001,p=15,
Age=c(11, 9, 11, 5, 13, 4, 11, 2, 9, 16, 11, 9, 2, 13, 2, 4, 3, 1, 7, 3, 6, 13, 2, 2, 2, 8, 0.37, 12, 1, 10, 13, 2, 2, 4, 13, 4,
2, 5, 2.7, 2.3, 8, 13, 9, 2, 3.5, 12.5, 4, 3, 7, 3, 11, 13, 7,
4, 5, 5, 1, 4, 10, 3, 11, 5, 0.4, 3, 4, 2, 6, 2.5, 6, 2, 12,
3, 12, 6, 3, 3, 2, 11, 13, 2.5, 6, 7, 3, 3, 9, 0.83, 4, 2, 7,
2, 11, 2, 5, 4, 3, 4, 2, 2, 14, 2, 1, 3, 9, 2, 1, 3, 2, 4, 8,
5, 2, 3, 5, 3, 4, 2, 6, 9, 2, 1, 9, 4, 3, 7, 4, 12, 11, 2, 3,
13, 2, 10, 5, 8, 11, 8, 2, 9, 9, 3, 3, 6, 5, 7, 4, 5, 3, 4, 2,
2, 3, 9, 5, 4, 14, 2, 8, 2, 9, 3, 5, 4, 9, 4, 5, 9, 1, 1, 3,
2, 14, 5, 4, 3, 7, 1, 8, 13, 3, 4, 10, 3, 7, 7, 4, 6, 12, 4,
4, 3, 5, 5, 3, 3, 9, 7, 5, 3, 14, 8),
Gender=c(1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1, 0, 1,
0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1,
1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0,
0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1,
0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1,
0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1,
0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0,
0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0,
1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1,
0, 1, 1, 0, 0, 0, 1, 0, 1),
Treatment.course=c(1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1,
1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0),
WBC=c(3, 1, 1, 1, 1, 2, 1, 1, 2, 1, 4, 4, 1, 1, 3,
1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 4, 2, 4, 2, 1,
1, 1, 4, 1, 3, 1, 1, 1, 2, 1, 4, 1, 2, 1, 3, 2,
1, 1, 1, 3, 1, 2, 1, 1, 1, 4, 1, 1, 1, 4, 2, 3,
1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 1,
2, 1, 4, 2, 4, 1, 2, 1, 3, 1, 2, 1, 2, 1, 1, 2,
1, 1, 1, 4, 4, 3, 2, 1, 2, 3, 2, 2, 1, 1, 1, 1,
1, 2, 2, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 0, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1,
1, 1, 3, 2, 1, 1, 2, 1, 1, 3, 1, 1, 1, 2, 2, 1,
4, 1, 2, 1, 1, 1, 2, 1, 1, 2, 1, 4, 1, 1, 3, 1,
2, 2, 1, 1, 1, 1, 1, 1, 4),
plt=c(21000, 19000, 82000, 44000, 298000, 49000, 38000, 55000,
15000, 56000, 26000, 350000, 12000, 297000, 1000, 90000,
50000, 3000, 10000, 39000, 10000, 30000, 72000, 109000,
389000, 17000, 94000, 31000, 27000, 175000, 148000, 57000,
102000, 76000, 300000, 29000, 18000, 11000, 17000, 46000,
11000, 65000, 402000, 23000, 40000, 178000, 63000, 20000,
58000, 86000, 40000, 183000, 19000, 167000, 10000, 28000,
216000, 211000, 32000, 548000, 13000, 4000, 1600, 11000,
25000, 17000, 7000, 22000, 45000, 93000, 481000, 97000,
76000, 47000, 3000, 303000, 42000, 347000, 99000, 52000,
41000, 45000, 173000, 17000, 156000, 415000, 17000, 70000,
80000, 45000, 422000, 330000, 79000, 423000, 7000, 400000,
3000, 31000, 16000, 16000, 58000, 17000, 123000, 41000,
16000, 8000, 116000, 9000, 140000, 15000, 36000, 36000,
40000, 29000, 38000, 3000, 12000, 70000, 200000, 28000,
6000, 52000, 31000, 35000, 8000, 93000, 136000, 45000,
81000, 0, 330000, 38000, 93000, 30000, 147000, 29000,
190000, 27000, 175000, 122000, 79000, 20000, 35000, 29900,
186000, 32000, 65000, 56000, 211000, 359000, 129000, 229000,
200000, 290000, 127000, 126000, 40000, 55000, 250000,
292000, 1000, 25000, 66000, 97000, 38000, 257000, 384000,
178000, 93000, 148000, 14000, 127000, 10000, 28000, 38000,
31000, 10000, 67000, 12000, 9000, 173000, 10000, 7000,
32000, 17000, 30000, 49000, 96000, 22000, 25000, 59000,
57000, 57000, 10000, 95000, 161000, 32000, 18000, 112000,
30000),
LDH=c(236, 0, 604, 869, 823, 1549, 0, 403, 0, 496, 744,
3194, 697, 0, 1213, 1030, 0, 649, 346, 0, 1596, 240,
0, 0, 0, 2151, 480, 3699, 2731, 2532, 7900, 30, 300,
5666, 1577, 867, 733, 900, 533, 600, 0, 12402, 422,
1487, 879, 4231, 1918, 225, 661, 214, 0, 716, 1054,
0, 0, 2372, 0, 0, 0, 0, 0, 964, 1961, 0, 1554, 427,
848, 0, 0, 872, 335, 807, 1304, 665, 508, 0, 1451,
0, 862, 1078, 0, 0, 351, 1341, 0, 824, 0, 0, 617,
1245, 1256, 0, 0, 748, 0, 362, 410, 539, 642, 0, 1038,
615, 647000, 627, 1560, 1104, 500, 660, 1859, 566, 0,
0, 1399, 439, 1145, 512, 961, 1329, 0, 972, 0, 1419,
0, 471, 488, 0, 648, 6500, 1264, 0, 620, 408, 1342,
427, 719, 0, 840, 472, 282, 744, 8740, 484, 622, 1758,
492, 1052, 0, 642, 1656, 0, 1036, 538, 450, 0, 462,
1209, 0, 707, 476, 766, 602, 1987, 605, 3764, 1182,
516, 178, 384, 0, 618, 0, 2061, 886, 1289, 647, 1725,
0, 4599, 922, 442, 418, 1726, 410, 859, 1102, 1650,
1748, 435, 723, 1988, 1516, 1112, 0, 638, 856, 0, 490,
1019, 638, 8330),
The.cell.line=c(1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2,
1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,
1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,
1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1,
2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1,
1, 1, 0, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1,
1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1),
CNS=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0),
hemorraahagic=c(0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0),
mediastinal.mass=c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0),
reumatoid.signs=c(0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1,
0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1,
1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0,
0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1,
0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1,
0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1,
0, 1, 1, 1, 1, 1, 0, 0, 1),
tumor.lysis.synd=c(0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0),hepatosplenomegalia=c(0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1,
0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1,
0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0,
0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1,
1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0,
0, 1, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1,
1, 0, 1, 1, 0, 1, 0, 1, 1),
Lymphadenopathy=c(0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0,
1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1,
0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1,
1, 0, 0, 1, 0, 0, 0, 1, 0),
Radiotherapy=c(0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1,
1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1,
1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 1, 1,
1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1,
1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0,
0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1,
1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1,
0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1,
0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1,
0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 0, 1),
obs.t=c(22, 24, 26, 31, 42, 44, 48, 74, 81, 92, 98, 98,
138, 138, 140, 141, 142, 146, 149, 171, 180, 200, 201,
219, 226, 234, 271, 279, 287, 287, 288, 373, 373, 381,
395, 397, 410, 420, 420, 422, 445, 457, 461, 462, 468,
472, 486, 486, 507, 511, 520, 520, 544, 566, 604, 615,
639, 647, 650, 683, 688, 705, 725, 749, 752, 765, 766,
769, 793, 793, 913, 915, 930, 930, 937, 944, 947, 961,
976, 994, 997, 1012, 1013, 1023, 1025, 1035, 1053, 1066,
1067, 1093, 1101, 1134, 1159, 1185, 1187, 1210, 1223,
1224, 1330, 1330, 1351, 1394, 1401, 1402, 1411, 1412,
1417, 1422, 1464, 1467, 1479, 1493, 1502, 1534, 1548,
1559, 1592, 1605, 1637, 1649, 1659, 1666, 1666, 1673,
1676, 1696, 1712, 1716, 1724, 1745, 1759, 1763, 1794,
1796, 1806, 1826, 1828, 1830, 1831, 1844, 1853, 1858,
1859, 1887, 1901, 1908, 1927, 1930, 1938, 1938, 1956,
1957, 1988, 2007, 2015, 2019, 2031, 2075, 2091, 2103,
2106, 2107, 2139, 2166, 2173, 2199, 2200, 2200, 2209,
2210, 2217, 2220, 2306, 2310, 2328, 2358, 2359, 2390,
2402, 2435, 2440, 2448, 2457, 2520, 2612, 2623, 2641,
2642, 2711, 2748, 2777, 2791, 2839, 2844, 2861, 2878,
2897, 2904, 2925, 2932),
t=c(22, 26, 48, 81, 92, 98, 138, 140, 146, 171, 180,
219, 226, 234, 271, 287, 288, 373, 381, 445, 520, 566,
604, 650, 683, 688, 752, 793, 913, 961, 1134, 1210,
1351, 1411, 1712, 2199),
fail=c(2, 0, 2, 0, 0, 0, 2, 0, 1, 2, 0, 1, 0, 2, 1,
0, 0, 2, 0, 2, 1, 0, 0, 2, 2, 1, 2, 0, 1, 1, 1,
2, 2, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 0, 2, 1, 0, 0, 0, 1, 2, 1, 0, 0,
0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0,
1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0)
)
-------------------------------------------------------------------
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
