model
{
for (k in 1:n) {

y[k] ~ dnorm(mu[k], tau[k])

mu[k] <- A[i[k]] + B[i[k]] * D[j[k]]

  tau[k]~dunif(0,1000)

    }

  for (i in 1:m){

     B[i]~ dunif(0,1000)

     A[i]~dunif(-1000,1000)
}
for (j in 1:m){
D[j]~ dunif(0,1000)
}
}