最近在自学WINBUGS系统，写了个模型（如下显示），但在load data 的过程中总是显示“ data loaded ( variables not in the model: GROUPS, TRIALS, TREATMENTS)（如图），请各位高手帮帮忙，谢谢啦！


model  
{
for (i in 1:GROUPS ) { logit(p)<-mu[s]+ delta  * (1-equals(t,b))
    r~dbin(p,n) }               #binomial likelihood
# random effect distributions of trial-specific LORs
# means of distributions
   delta~dnorm (md, tau)
   md <- d[t - d[b]
# vague priors for 24 trial baselines
for (j in 1:TRIALS ) { mu[j]~dnorm(0, 0.0001) }

# vague priors for basic parameters
d[1] <-0
for (k in 2:TREATMENTS ) {d[k]~dnorm(0, 0.0001) }

# vague prior for RE standard deviation
sd~dunif(0,2)
tau< -1/pow(sd,2)
}

# DATA
list(GROUPS=50, TRIALS=24, TREATMENTS=4 ),
s[ ]     t[ ]     r[ ]     n[ ]     b[ ]
1       1        79     702   1
1       2        77     694   1
2       1        18     671   1
2       2        21     535   1
3       1        8       161   1
3       2       19      146   1
4       1       75      731   1
4       3       363    714   1
5       1       2        106   1
5       3       9        205   1
6       1       58      549   1
6       3       237    1561 1
7       1       0        33     1
7       3       9        48     1
8       1       3        100   1
8       3       31      98     1
9       1       1        31     1
9       3       26      95     1
10     1       6        39     1
10     3       17      77     1
11     1       64      642   1
11     3       107    761   1
12     1       5        62     1
12     3       8        90     1
13     1       20      234   1
13     3       34      237   1
14     1       95      1107 1
14     3      134     1031 1
15     1      15       187   1
15     3      35       504   1
16     1      78       584   1
16     3      73       675   1
17     1      69       1177 1
17     3      54       888   1
18     1      9         140   1
18     3      23       140   1
19     1      0         20     1
19     4      9         20     1
20     2      20       49     2
20     3      16       43     2
21     2      11       78     2
21     3      12       85     2
21     4       29      170   2
22     2       7        66     2
22     4       32      127   2
23     3       12      76     3
23     4       20      74     3
24     3       9        55     3
24     4       3        26     3
END
list(p=c( 0.5, NA, NA, NA), t=c( 1, 2, 3, 4 ), b=c( 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3 ) ])