model

           {
     
           for (i in 1:1) {          # one subject
           c~dnorm(mu, 1.0E+6)       # c:drug concentration
           mu<-D/(cl_F*bw*tau)*1000  # oral Pediatric Digoxin concentration equation at steady state

           cl_F~dnorm(theta1,100)

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

           }