model

           {
     
           for (i in 1:N) {            # N:number of patients
           c[i]~dnorm(mu[i], 1.0E+6)   # c:drug concentration
           mu[i]<-(D[i]/36.9458)/(cl_F[i]*tau[i])   # oral Lthium carbonate concentration equation at steady state

           cl_F[i]~dnorm(theta1[i],  7.032131)

           theta1[i]<-(31.6-0.634*(age[i]-50)*f[i]+(-7.79+0.225*bw[i])/Scr[i])/24  # population estimate of Lithium carbonate clearance
  
           }

           }