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

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

           theta1[i]<-0.237*(1+0.094*age[i])    # population estimate of Pediatric Digxoin clearance
  
           }

           }