model { c <- 10000 for (i in 1:N) { ones[i] <- 1 ones[i] ~ dbern(prob[i]) prob[i] <- abs(a)*pow(y[i],a*p-1)/(pow(b[i],a*p)*exp(loggam(p)+loggam(q)-loggam(p+q))*pow(1+pow(y[i]/b[i],a),p+q)*c) b[i]<- exp(mu + alpha[row[i]]+beta[col[i]])*exp(loggam(p)+loggam(q)-loggam(p+1/a)-loggam(q-1/a)) } for (j in 2:5){ beta[j]~ dnorm(0, 0.01) alpha[j] ~ dnorm(0, 0.01) } alpha[1]<- 0-sum(alpha[2:5]) beta[1]<- 0-sum(beta[2:5]) mu ~ dnorm(0, 0.01) a ~ dnorm(0, 0.01)I(-10,-1) q ~ dgamma(0.01,0.001)I(0,90) p ~ dgamma(0.01,0.001)I(1,3.5) } list(a=-2.03, q=1.5,p=1,mu=8,alpha=c(NA,0,0,0,0), beta=c(NA,0,0,0,0)) list(N=15) y[] row[] col[] 100000 1 1 700000 1 2 2000000 1 3 4100000 1 4 7100000 1 5 11300000 1 6 17200000 1 7 22800000 1 8 27700000 1 9 34400000 1 10 40600000 1 11 54900000 1 12 67900000 1 13 77700000 1 14 85600000 1 15 94500000 1 16 102500000 1 17 106500000 1 18 115900000 1 19 119300000 1 20 122600000 1 21 124800000 1 22 130200000 1 23 100000 2 1 700000 2 2 1600000 2 3 2500000 2 4 3800000 2 5 5500000 2 6 8600000 2 7 11200000 2 8 13900000 2 9 19800000 2 10 25900000 2 11 33300000 2 12 41200000 2 13 49700000 2 14 58800000 2 15 68100000 2 16 73500000 2 17 80100000 2 18 88500000 2 19 93800000 2 20 95200000 2 21 97700000 2 22 100000 3 1 700000 3 2 1400000 3 3 2200000 3 4 3300000 3 5 5000000 3 6 7600000 3 7 10600000 3 8 14900000 3 9 20900000 3 10 27900000 3 11 33700000 3 12 40800000 3 13 51000000 3 14 62000000 3 15 67600000 3 16 75700000 3 17 81900000 3 18 86700000 3 19 88000000 3 20 95300000 3 21 1000000 4 1 700000 4 2 1500000 4 3 2600000 4 4 3900000 4 5 5900000 4 6 8000000 4 7 11600000 4 8 16000000 4 9 24000000 4 10 29800000 4 11 35200000 4 12 46200000 4 13 55800000 4 14 62200000 4 15 72700000 4 16 80100000 4 17 86300000 4 18 91700000 4 19 95800000 4 20 0 5 1 600000 5 2 1500000 5 3 2500000 5 4 3500000 5 5 4600000 5 6 7100000 5 7 10800000 5 8 16100000 5 9 20600000 5 10 28000000 5 11 34000000 5 12 46100000 5 13 51400000 5 14 61200000 5 15 69200000 5 16 75900000 5 17 80300000 5 18 88500000 5 19 100000 6 1 600000 6 2 1600000 6 3 2500000 6 4 3400000 6 5 5100000 6 6 7400000 6 7 11600000 6 8 16600000 6 9 24700000 6 10 32100000 6 11 40800000 6 12 48800000 6 13 58600000 6 14 67300000 6 15 79100000 6 16 89100000 6 17 94700000 6 18 100000 7 1 600000 7 2 1400000 7 3 2100000 7 4 3100000 7 5 5000000 7 6 7500000 7 7 11600000 7 8 17400000 7 9 23000000 7 10 30200000 7 11 36400000 7 12 46400000 7 13 52900000 7 14 64400000 7 15 71200000 7 16 76600000 7 17 100000 8 1 600000 8 2 1300000 8 3 2200000 8 4 3300000 8 5 5400000 8 6 8600000 8 7 13400000 8 8 19200000 8 9 26700000 8 10 34700000 8 11 48900000 8 12 58300000 8 13 73200000 8 14 81700000 8 15 93300000 8 16 100000 9 1 600000 9 2 1400000 9 3 2200000 9 4 3300000 9 5 5100000 9 6 9200000 9 7 15100000 9 8 21600000 9 9 27500000 9 10 39900000 9 11 49200000 9 12 63100000 9 13 69200000 9 14 79200000 9 15 0 10 1 500000 10 2 1500000 10 3 2300000 10 4 3400000 10 5 4900000 10 6 9000000 10 7 15700000 10 8 21200000 10 9 31800000 10 10 39700000 10 11 47900000 10 12 56600000 10 13 65300000 10 14 100000 11 1 700000 11 2 1500000 11 3 2400000 11 4 3700000 11 5 6000000 11 6 11300000 11 7 15700000 11 8 23100000 11 9 31400000 11 10 40500000 11 11 47900000 11 12 63100000 11 13 100000 12 1 700000 12 2 1800000 12 3 2700000 12 4 5000000 12 5 7900000 12 6 11700000 12 7 19900000 12 8 28800000 12 9 38300000 12 10 45200000 12 11 53800000 12 12 100000 13 1 800000 13 2 1800000 13 3 3000000 13 4 4500000 13 5 7600000 13 6 13300000 13 7 20400000 13 8 29100000 13 9 36900000 13 10 46800000 13 11 0 14 1 500000 14 2 1100000 14 3 2000000 14 4 2900000 14 5 6100000 14 6 10800000 14 7 18000000 14 8 24800000 14 9 31800000 14 10 100000 15 1 600000 15 2 1400000 15 3 2100000 15 4 4100000 15 5 8500000 15 6 15900000 15 7 22400000 15 8 30100000 15 9 0 16 1 600000 16 2 1300000 16 3 2200000 16 4 4000000 16 5 10200000 16 6 17100000 16 7 24800000 16 8 100000 17 1 500000 17 2 1500000 17 3 3000000 17 4 5000000 17 5 9400000 17 6 17000000 17 7 0 18 1 700000 18 2 1700000 18 3 2700000 18 4 4300000 18 5 9900000 18 6 100000 19 1 600000 19 2 1400000 19 3 2300000 19 4 4300000 19 5 100000 20 1 800000 20 2 1700000 20 3 2800000 20 4 100000 21 1 600000 21 2 1300000 21 3 100000 22 1 700000 22 2 100000 23 1 END