Как получить апостериорные резюме для 500 образцов начальной загрузки одновременно, используя Openbugs? - PullRequest
0 голосов
/ 28 октября 2019

Я использовал пакет Openbugs для запуска байесовской модели лечения с логистическим распределением и получения апостериорных сводок. Я могу получить последующие сводки для одного набора данных успешно, но у меня есть вопрос здесь. Как запустить модель для 500 образцов начальной загрузки, извлеченных из реального набора данных одновременно? Я представил свои данные и модель ниже.

################################model
model { 
for (i in 1:252) {
Sstar[i] <- p[i]+(1-p[i])*S0[i]; fstar[i] <- (1-p[i])*f0[i]
S0[i] <- 1/(1+lam[i]*pow(t[i],kap))
log(f0[i]) <- log(kap)+log(lam[i])+(kap-1)*log(t[i])-2*log(1+lam[i]*pow(t[i],kap))
log(lam[i]) <-beta[1]+beta[2]*age[i]+beta[3]*sex[i]+beta[4]*bmi[i]+beta[5]*smoking[i]+beta[6]*Diabetes[i]+
beta[7]*hypertension[i]+beta[8]*diagnosis.age[i]+beta[9]*hour.w[i]+beta[10]*time.w[i]+beta[11]*ch.k[i]+beta[12]*job.less[i]+beta[13]*job.house[i]+beta[14]*job.farm[i]+beta[15]*job.retired[i]+beta[16]*job.oth[i]+beta[17]*blood.a[i]+
beta[18]*blood.b[i]+beta[19]*blood.ab[i]+beta[20]*marital_r[i]


p[i] <- 1/(1+exp(-(eta[1]+eta[2]*age[i]+eta[3]*sex[i]+eta[4]*bmi[i]+eta[5]*smoking[i]+eta[6]*Diabetes[i]+
eta[7]*hypertension[i]+eta[8]*diagnosis.age[i]+eta[9]*hour.w[i]+eta[10]*time.w[i]+eta[11]*ch.k[i]+eta[12]*job.less[i]+
eta[13]*job.house[i]+eta[14]*job.farm[i]+eta[15]*job.retired[i]+eta[16]*job.oth[i]+eta[17]*blood.a[i]+
eta[18]*blood.b[i]+eta[19]*blood.ab[i]+eta[20]*marital_r[i])))
# log-likelihood#
log(L[i]) <- d[i]*log(fstar[i])+(1-d[i])*log(Sstar[i]);logL[i] < - log(L[i])
zeros[i] < - 0
zeros[i] ~ dloglik(logL[i])}
Dv< - -2*sum(logL[])
# priors
for (j in 1:20) {eta[j]~dnorm( 0, 300)}
for (k in 1:20) {beta[k]~dnorm( 0,300) }
kap ~ dunif(0, 100)}

###############################data(dialysisthesis.use1.txt)
age[]   sex[]   bmi[]   smoking[]   Diabetes[]  hypertension[]  diagnosis.age[] hour.w[]    time.w[]    ch.k[]  job.less[]  job.house[] job.farm[]  job.retired[]   job.oth[]   blood.a[]   blood.b[]   blood.ab[]  marital_r[] t[] d[]
18  1   15.22   0   0   1   15  4.00    2   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    3   0
45  1   21.45   0   1   0   43  4.00    3   5.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    2   0
64  1   21.08   1   1   1   44  4.00    3   12.00   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    20  0
64  1   22.77   0   0   1   34  3.50    2   16.00   0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    30  0
60  1   29.07   0   1   1   52  4.00    3   4.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    1.00    8   0
59  1   24.16   0   1   0   49  4.00    3   6.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    10  0
58  1   25.53   0   1   0   36  4.00    2   6.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    22  0
55  1   18.51   0   0   0   45  4.00    3   6.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
56  1   22.49   1   1   0   46  4.00    1   5.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
30  1   19.53   0   0   1   20  4.00    3   14.00   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    10  0
83  1   22.49   1   1   1   73  4.00    3   12.00   0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    10  0
38  1   20.00   0   0   0   35  3.50    3   12.00   0.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    3   0
36  1   20.31   0   1   0   31  4.00    3   8.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    5   0
85  1   23.44   1   1   1   75  4.00    3   12.00   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
47  1   22.27   1   0   1   37  4.00    3   8.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    10  0
17  1   16.65   0   0   1   12  3.50    3   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    5   0
87  1   21.48   1   1   0   77  4.00    3   5.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    10  0
46  1   24.97   1   0   0   41  4.00    3   18.00   0.00    0.00    0.00    0.00    1.00    1.00    0.00    0.00    1.00    5   0
66  1   23.18   1   0   0   56  4.00    2   0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    10  0
56  1   20.76   0   1   0   46  4.00    3   5.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    10  0
53  1   22.04   1   0   0   43  4.00    2   5.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    10  0
72  1   22.49   1   1   0   62  4.00    3   5.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
27  1   17.78   0   0   0   24  4.00    3   12.00   1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    3   0
31  1   20.00   0   0   1   26  4.00    3   4.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    5   0
72  1   25.26   0   1   1   64  4.00    2   12.00   0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    8   0
24  1   17.82   0   0   0   14  3.50    3   12.00   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    10  0
28  1   23.05   0   0   0   18  4.00    3   5.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    10  0
51  1   23.15   0   1   0   38  4.00    3   7.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    13  0
71  1   22.66   1   1   0   61  3.50    3   5.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    1.00    10  0
69  1   25.61   0   1   0   58  4.00    3   5.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    11  0
18  1   18.88   0   0   0   1   4.00    2   10.00   1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    17  0
38  1   29.69   0   0   1   28  4.00    2   8.00    0.00    0.00    0.00    0.00    1.00    1.00    0.00    0.00    0.00    10  0
48  1   23.53   0   0   1   42  3.00    3   5.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    1.00    6   0
41  1   22.86   0   0   1   33  4.00    3   0.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    8   0
26  1   26.56   1   0   1   21  4.00    3   5.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    5   0
67  1   25.39   1   1   1   47  3.00    3   12.00   0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    20  0
45  1   20.76   1   1   0   35  4.00    3   4.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
79  1   25.71   1   0   1   59  4.00    3   4.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    20  0
10  1   17.36   0   0   0   10  4.00    2   0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1   0
25  1   17.86   1   0   0   15  4.00    2   12.00   1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
35  1   26.63   1   0   1   30  4.00    2   0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    5   0
17  1   17.86   0   0   0   14  4.00    2   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    3   0
50  1   22.04   1   0   0   40  4.00    2   0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
62  1   27.34   1   1   1   42  4.00    2   12.00   0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    20  0
45  1   24.97   1   1   1   25  4.00    2   10.00   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    20  0
50  1   24.97   1   0   1   40  3.00    3   2.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
70  1   21.30   1   1   1   64  4.00    2   12.00   1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    6   0
39  1   28.91   0   0   1   37  4.00    3   4.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    2   0
65  1   24.58   1   0   1   59  4.00    3   0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    6   0
56  1   20.44   0   0   0   33  4.00    3   5.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    23  0
57  1   20.34   0   0   0   55  4.00    3   18.00   0.00    0.00    0.00    0.00    1.00    1.00    0.00    0.00    1.00    2   0
58  1   24.84   0   1   1   48  4.00    3   0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    10  0
42  1   21.64   0   1   1   35  4.00    3   9.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    7   0
60  1   21.08   1   0   1   45  4.00    1   8.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    15  0
48  1   18.73   1   0   1   38  4.00    3   8.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
27  1   16.79   0   0   1   26  4.00    2   8.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1   0
44  1   29.04   1   1   1   32  4.00    2   0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    12  0
43  1   23.88   1   0   1   33  4.00    2   4.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
54  1   24.80   0   1   1   51  4.00    3   9.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    3   0
53  1   22.31   0   1   1   41  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    12  0
51  1   21.88   0   0   0   49  4.00    3   0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    2   0
35  1   24.91   0   0   1   25  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
63  1   27.18   0   1   0   53  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    10  0
50  1   25.51   1   1   0   40  4.00    3   8.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    10  0
75  1   23.44   1   1   1   65  4.00    3   16.00   1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
72  1   24.22   0   0   0   68  4.00    3   12.00   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    4   0
71  1   18.62   0   1   1   61  4.00    3   12.00   0.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    1.00    10  0
43  1   18.69   1   0   1   28  4.00    3   5.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    15  0
70  1   25.39   1   1   1   60  4.00    2   12.00   0.00    0.00    0.00    0.00    1.00    1.00    0.00    0.00    1.00    10  0
65  1   24.03   0   1   0   55  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
53  1   30.10   0   1   0   42  4.00    3   9.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    11  0
35  1   24.03   0   0   0   59  4.00    2   0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    24  0
68  1   25.28   1   1   1   56  4.00    2   5.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    12  0
25  1   24.34   0   0   1   24  4.00    3   4.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    1   0
53  1   35.59   1   0   1   43  4.00    3   8.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
57  1   27.43   1   0   1   49  4.00    3   8.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    8   0
27  1   19.72   0   0   1   22  4.00    3   12.00   0.00    0.00    0.00    0.00    1.00    1.00    0.00    0.00    1.00    5   0
48  1   24.97   1   0   1   43  4.00    3   12.00   0.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    1.00    5   0
75  1   30.92   0   1   1   55  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    20  0
68  1   22.66   1   1   1   57  4.00    3   12.00   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    11  0
75  1   23.44   0   0   0   70  4.00    3   12.00   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    5   0
40  1   24.44   0   1   0   35  4.00    3   0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    5   0
38  1   20.00   0   0   1   33  4.00    3   5.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    5   0
45  1   22.04   1   1   0   40  4.00    3   0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    5   0
36  1   18.73   0   1   0   31  4.00    3   0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    5   0
81  1   25.95   0   1   1   60  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    21  0
58  1   22.04   1   1   0   43  4.00    3   0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    15  0
71  1   24.22   0   1   0   61  4.00    3   12.00   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    10  0
63  1   22.49   1   1   1   43  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    20  0
53  1   24.61   1   0   1   48  4.00    3   8.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    1.00    5   0
64  1   21.80   0   0   1   54  4.00    3   0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
65  1   25.39   0   1   1   42  4.00    3   5.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    23  0
76  1   22.83   0   0   1   74  4.00    2   0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    2   0
64  1   28.06   0   1   0   54  4.00    3   0.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    10  0
78  1   19.83   0   0   1   48  4.00    3   12.00   0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    30  0
30  1   35.16   1   1   1   20  4.00    3   9.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    10  0
51  1   19.87   1   0   1   44  4.00    3   5.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    7   0
68  1   20.57   1   1   0   53  4.00    3   0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    15  0
49  1   27.34   1   1   1   48  4.00    3   0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    1   0
86  1   22.49   1   1   1   79  4.00    3   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    7   0
73  1   29.30   1   0   0   58  4.00    3   5.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    15  0
42  1   25.93   1   0   0   40  4.00    3   9.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    2   0
60  1   19.31   1   0   1   50  4.00    3   5.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
52  1   24.22   0   0   1   48  4.00    3   6.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    4   0
64  1   23.63   0   1   0   56  4.00    3   10.00   0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    8   0
63  1   24.73   0   1   1   53  4.00    3   5.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    1.00    10  0
51  1   22.04   1   1   0   48  3.50    3   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    3   0
55  1   25.95   0   1   1   20  4.00    3   0.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    35  0
74  1   20.00   0   0   1   73  4.00    3   0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    0.00    1.00    1   0
65  1   24.61   1   1   1   30  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    35  0
63  1   24.91   0   0   1   58  3.00    3   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    5   0
63  1   27.68   0   1   1   58  3.00    2   9.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    5   0
63  1   29.05   0   1   1   51  4.00    3   4.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    12  0
87  1   24.57   1   1   1   83  3.00    3   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    1.00    4   0
15  1   23.67   0   0   1   12  3.00    3   6.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    3   0
84  1   20.89   0   1   1   64  4.00    3   12.00   1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    20  0
57  1   16.65   1   1   1   45  3.00    2   6.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    1.00    12  0
78  1   20.76   1   0   0   73  4.00    2   5.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    1.00    5   1
62  1   22.04   0   1   0   52  4.00    3   0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    0.00    1.00    10  1
40  1   20.76   1   1   0   30  4.00    3   1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    10  1
17  1   20.00   0   0   0   12  3.00    3   6.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    5   1
87  1   21.60   1   1   1   72  3.00    2   1.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    15  1
60  1   20.76   1   1   0   40  4.00    3   5.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    20  1
11  1   20.83   0   0   0   8   4.00    3   5.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    3   1
14  1   13.89   0   0   0   12  3.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    2   1
45  1   24.21   0   1   1   37  4.00    3   0.00    0.00    0.00    0.00    0.00    1.00    1.00    0.00    0.00    1.00    8   1
91  1   17.90   0   0   1   86  4.00    3   9.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    5   1
78  1   17.04   1   0   0   48  3.00    3   9.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    30  0
60  1   18.73   1   1   1   40  4.00    3   5.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    20  1
57  1   25.71   0   1   1   37  3.00    3   16.00   0.00    0.00    0.00    1.00    0.00    1.00    0.00    0.00    1.00    20  1
30  1   17.78   1   1   1   25  3.00    2   5.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    5   1
53  1   22.04   1   0   1   43  4.00    3   12.00   1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  1
69  1   23.88   1   1   1   59  4.00    3   12.00   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    10  1
65  1   23.44   1   1   0   55  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  1
64  1   26.64   1   0   1   59  3.00    3   9.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    0.00    1.00    5   1
70  1   31.11   0   1   0   65  4.00    3   5.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    5   1
61  2   19.56   0   1   1   51  3.00    3   16.00   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
33  2   25.10   0   0   1   23  4.00    3   6.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
65  2   28.91   1   1   0   41  3.00    3   12.00   0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    24  0
58  2   16.02   1   1   1   40  3.00    3   9.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    18  0
53  2   21.50   0   1   0   48  3.00    3   0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    1.00    1.00    5   0
37  2   33.59   0   1   0   29  4.00    3   8.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    8   0
62  2   19.10   0   1   1   54  4.00    3   0.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    0.00    1.00    8   0
41  2   23.12   0   0   0   2   3.00    2   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    39  0
31  2   19.05   0   0   1   29  4.00    3   0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    2   0
70  2   25.00   0   1   1   50  4.00    3   12.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    20  0
45  2   17.78   0   0   1   40  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    5   0
21  2   17.78   0   0   1   16  4.00    2   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    5   0
55  2   22.04   0   1   0   45  3.00    3   5.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
44  2   21.48   0   0   0   39  3.00    2   8.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    5   0
45  2   18.73   0   0   1   35  3.00    2   10.00   0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
59  2   19.56   0   1   1   49  3.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
44  2   23.14   0   0   1   36  4.00    2   2.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    8   0
17  2   21.50   0   0   0   12  4.00    2   14.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    5   0
71  2   20.81   0   0   1   61  4.00    3   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
36  2   19.82   0   1   1   36  3.00    3   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    1   0
56  2   24.82   0   1   1   46  4.00    2   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
47  2   18.36   0   1   1   37  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    10  0
63  2   17.58   0   1   1   53  3.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    10  0
36  2   22.27   0   0   1   26  4.00    2   9.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
28  2   23.14   0   0   0   18  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    10  0
38  2   30.67   0   0   1   27  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    11  0
36  2   35.26   0   0   1   35  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    1   0
14  2   20.41   0   0   1   12  4.00    2   5.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    2   0
16  2   20.71   0   0   1   12  4.00    3   14.00   1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    4   0
34  2   16.46   0   0   0   24  4.00    2   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
70  2   25.39   1   1   0   45  4.00    2   12.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    25  0
50  2   23.78   1   1   0   30  3.00    2   4.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    20  0
49  2   22.22   1   1   0   39  3.00    2   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
48  2   26.67   1   0   0   38  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
64  2   28.54   0   1   1   34  4.00    3   5.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    30  0
35  2   23.46   0   0   0   30  4.00    2   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    5   0
12  2   20.83   0   0   0   1   3.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    11  0
20  2   17.75   0   0   0   15  4.00    3   8.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    5   0
56  2   21.64   1   0   0   52  4.00    3   6.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    4   0
58  2   24.03   1   1   1   45  3.00    2   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    13  0
65  2   26.14   0   0   1   50  4.00    2   9.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    15  0
59  2   20.42   1   1   1   52  4.00    2   12.00   0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    7   0
56  2   23.83   0   1   1   34  4.00    2   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    22  0
56  2   23.51   0   1   0   84  4.00    2   12.00   0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    28  0
64  2   25.80   0   0   1   62  3.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    2   0
74  2   27.27   0   0   1   51  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    23  0
62  2   23.94   0   0   1   34  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    28  0
72  2   22.83   1   0   1   55  4.00    3   12.00   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    17  0
58  2   19.98   0   0   1   48  3.50    3   12.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    1.00    10  0
71  2   26.30   0   1   1   54  3.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    17  0
26  2   18.73   0   0   0   22  4.00    3   10.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    4   0
44  2   26.37   0   0   1   24  3.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    20  0
46  2   23.44   1   0   0   36  4.00    3   2.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    10  0
54  2   16.85   1   0   1   53  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    1   0
32  2   15.82   0   0   1   22  4.00    2   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
71  2   19.56   1   1   0   46  4.00    2   5.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    25  1
71  2   27.73   1   1   1   56  4.00    2   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    15  0
50  2   23.71   1   1   1   34  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    16  0
32  2   24.61   0   0   1   28  4.00    3   12.00   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    4   0
74  2   23.05   0   1   1   66  4.00    2   5.00    0.00    0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    8   0
51  2   24.44   0   0   0   46  4.00    2   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    5   0
50  2   22.77   0   1   0   38  4.00    3   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    12  0
51  2   25.78   1   1   0   49  4.00    2   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    2   0
68  2   26.08   0   1   0   58  4.00    2   5.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
53  2   19.78   0   0   1   43  4.00    2   12.00   0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
74  2   17.19   0   1   1   65  4.00    2   2.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    9   0
37  2   18.37   0   1   0   32  4.00    2   12.00   0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    5   0
45  2   20.00   0   1   0   31  4.00    2   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    14  0
67  2   19.81   0   1   0   65  4.00    2   5.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    2   0
77  2   25.39   1   1   1   57  4.00    2   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    20  0
63  2   25.81   0   1   1   62  4.00    2   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    1   0
62  2   19.11   0   0   1   50  4.00    1   10.00   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    12  0
62  2   39.45   0   1   1   10  4.00    2   8.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    52  0
63  2   24.24   1   1   1   48  4.00    1   6.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    15  0
69  2   19.33   0   0   1   54  4.00    2   9.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    15  0
68  2   25.71   1   1   1   58  3.00    2   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  0
55  2   20.81   0   0   1   40  4.00    2   12.00   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    15  0
63  2   20.31   0   1   1   61  4.00    2   12.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    2   0
38  2   16.65   0   1   0   10  4.00    1   18.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    28  0
74  2   22.53   1   1   1   66  4.00    1   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    8   0
66  2   24.61   0   0   0   56  2.00    1   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    10  0
39  2   20.70   0   1   1   35  3.00    3   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    4   0
29  2   14.27   0   1   0   24  4.00    3   10.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    5   0
77  2   23.44   0   1   0   57  4.00    3   12.00   0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    20  0
60  2   22.49   0   0   1   50  3.00    2   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
72  2   26.56   0   1   1   62  3.00    2   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    10  0
66  2   19.47   0   0   0   56  4.00    2   8.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  0
42  2   45.00   0   0   1   40  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    1.00    2   0
54  2   20.24   1   1   1   26  3.00    2   0.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    28  0
45  2   29.72   0   1   0   28  3.00    2   5.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    17  0
95  2   21.09   1   0   1   65  3.00    3   6.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    30  0
49  2   17.07   0   0   1   39  3.00    3   5.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  0
63  2   38.93   1   1   0   50  4.00    3   14.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    13  0
57  2   26.04   0   1   1   37  3.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    1.00    20  0
28  2   23.78   0   0   1   25  4.00    2   8.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    3   0
80  2   18.67   1   1   0   78  4.00    3   12.00   0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    2   0
69  2   22.27   0   1   1   59  3.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    10  0
83  2   19.48   0   1   0   82  3.00    3   2.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    1   0
59  2   24.97   0   0   1   53  4.00    3   14.00   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    6   0
48  2   20.58   0   0   1   41  4.00    3   8.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    7   0
57  2   15.23   0   1   1   47  3.00    2   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  1
53  2   20.31   1   1   1   38  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    15  1
29  2   20.00   0   0   0   26  4.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    3   1
61  2   22.04   0   0   1   41  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    20  1
79  2   22.04   0   1   1   69  3.00    3   9.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    10  1
14  2   19.02   0   0   0   12  3.00    3   0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    2   1
26  2   17.86   0   1   0   16  3.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  1
50  2   21.08   0   1   0   40  3.00    3   5.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    10  1
81  2   27.34   0   0   1   56  3.00    3   4.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    25  1
28  2   17.78   0   1   0   23  4.00    3   5.00    1.00    0.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    5   1
63  2   27.34   1   1   1   48  4.00    3   5.00    0.00    1.00    0.00    0.00    0.00    1.00    0.00    0.00    1.00    15  1
81  2   26.67   1   1   1   56  4.00    3   6.00    0.00    1.00    0.00    0.00    0.00    0.00    0.00    0.00    1.00    25  1
34  2   17.78   0   1   0   24  4.00    3   7.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  1
53  2   23.78   1   1   0   43  4.00    2   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    10  1
83  2   21.64   0   1   1   63  4.00    3   0.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    0.00    20  1
59  2   21.48   0   1   1   44  4.00    3   6.00    0.00    1.00    0.00    0.00    0.00    0.00    1.00    0.00    1.00    15  1
END
...