model
{
    for(i in 1:n) {
                            for(j in 1:m) {
                                                           N[i, j]~dnorm(u[i], sigma[i])
                                                      }
                     }
             u[i]<-5.03+beta*log(theta[i])
            sigma[i]<-1/(0.052+beta*log(theta[i])-beta1*log(theta[i])
            theta[i]<-(300/s[i])
       }

       beta~dnorm(1.67, 0.00001)
       beta1~dnorm(1.87, 0.00001)
}

list(N=structure(.Data=c(5.603, 5.544, 5.528, 5.630, 5.594, 5.540, 5.581, 5.548, 5.426, 5.567, 5.554, 5.627, 5.630, 5.596, 5.626, 5.028, 5.074, 5.016, 4.894, 4.993, 5.071, 5.024, 5.035, 4.954, 5.039, 5.098, 5.057, 5.092, 5.082, 5.005, 4.784, 4.842, 4.776, 4.813, 4.813, 4.860, 4.798, 4.776, 4.758, 4.770, 4.755, 4.837, 4.736, 4.842, 4.796,  4.477, 4.400, 4.426, 4.462, 4.592, 4.411, 4.447, 4.402, 4.665, 4.475, 4.458, 4.551, 4.551, 4.525, 4.641), .Dim=c(4, 15)), s=c(200, 300, 350, 400),
     c=4, m=15)
list(beta=2, beta1=2)