model

           {
     
           for (i in 1:1){         # one subject
           c~dnorm(mu,100)         # c:drug concentration
           mu<-(D*24/tau)*Km/(Vmax-D*24/tau)    # oral phenytoin concentration equation at steady state 

           Vmax~dnorm(theta1, 0.0001)    #  Vmax~dnorm(theta1, 0.01) <- original
           Km~dnorm(theta2, 100)         #  Km~dnorm(theta2, 10) <- original

           theta1<-7.22*bw   # population estimate of Phenytoin Vmax
           theta2<-4.44      # population estimate of Phenytoin Km 
              
           }

           }