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]*1000/(cl_F[i]*tau[i])   # oral Digoxin concentration equation at steady state

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

           theta1[i]<-(8.03/24)*(1-0.0058*age[i])*bw[i]*pow(scr[i],-0.6)*pow(0.88,(1-Gender[i]))    # population estimate of Digxoin clearance
  
           }

           }