model { for (i in 1:n) { for (j in 1:t) { mu[i,j] <- beta0[i] + beta1[i]*(x[j]-x.mean); y[i,j] ~ dnorm(mu[i,j],tau); } beta0[i] ~ dnorm(mu.beta0, tau.beta0); beta1[i] ~ dnorm(mu.beta1, tau.beta1); } tau ~ dgamma(0.01, 0.01) mu.beta0 ~ dnorm(0, 0.01) tau.beta0 ~ dgamma(0.01, 0.01) mu.beta1 ~ dnorm(0, 0.01) tau.beta1 ~ dgamma(0.01, 0.01) x.mean <- mean(x[]); } list(n=6, t=44, y = structure(.Data = c(1.1933,1.2173,1.2544,1.2847,1.3198,1.3351,1.3601,1.4116,1.4517,1.4781,1.5126,1.5306,1.5427,1.5639,1.5798,1.586,1.6107,1.6435,1.6713,1.7154,1.7558,1.7931,1.8195,1.8477,1.8827,1.9398,1.976,2.0128,2.0441,2.0698,2.0979,2.1284,2.1632,2.2119,2.2651,2.325,2.3587,2.4056,2.4522,2.5051,2.5607,2.5729,2.5866,2.6127, 87.973,90.021,90.12,91.18,89.832,90.668,89.85,91.276,89.964,91.566,91.405,91.691,89.341,90.31,89.324,89.436,88.815,89.92,90.48,92.367,91.637,94.256,95.021,96.547,95.45,97.657,97.985,99.383,98.213,99.85,99.877,101.169,99.922,102.186,102.657,104.522,103.445,106.05,105.963,107.644,106.425,108.741,108.897,110.26, 9.392,13.637,14.392,13.301,10.411,13.057,13.621,13.354,11.056,14.559,14.067,12.286,9.658,13.659,12.351,12.307,10.508,16.287,16.697,15.523,13.54,19.19,18.827,17.215,13.848,20.319,20.403,19.179,16.79,24.4,24.215,22.042,18.182,22.67,21.204,19.746,16.656,24.552,23.775,22.449,18.145,26.344,25.278,22.511, 42.96,47.516,43.888,40.298,40.016,39.073,40.484,38.703,41.919,45.398,46.635,39.97,41.554,48.384,46.041,46.411,47.558,59.342,57.487,57.165,62.2,74.678,69.074,67.03,70.453,85.102,85.142,71.203,73.38,88.139,92.899,81.838,75.122,90.515,89.907,77.163,84.508,99.279,94.556,87.039,87.696,102.934,102.095,84.497, 9.203,10.174,11.078,12.334,10.228,10.394,11.246,13.071,10.954,11.41,11.892,13.204,10.439,11.148,11.556,13.37,11.114,12.043,13.165,15.534,13.002,14.602,15.684,19.088,15.436,16.117,16.782,20.491,16.383,18.18,19.782,23.605,18.483,19.659,21.042,24.203,20.242,22.143,22.87,27.486,22.399,22.412,22.223,25.648, 30.363,25.812,26.026,38.081,21.909,26.615,27.039,40.17,23.641,30.477,30.065,43.805,24.256,30.976,30.579,45.471,25.891,33.05,33.385,50.311,29.177,36.855,35.418,53.126,30.325,37.977,37.028,54.369,31.594,37.617,36.297,52.457,30.718,41.09,40.706,60.491,34.723,42.871,42.072,64.101,36.516,46.782,47.006,69.966), .Dim = c(6,44) ), x = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22, 23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44) ) list(beta0 = c(1,1,1,1,1,1), beta1 = c(1,1,1,1,1,1), tau = 1, mu.beta0 = 2, tau.beta0 = 1, mu.beta1 = 2, tau.beta1 = 1)