model

  {
 
  for (i in 1:1) {                                                         # one subject
  c~dnorm(mu,100)                                                          # c:drug concentration; one conc.
  mu<-(D*(1-exp(-cl/v*tin))/(tin*cl*(1-exp(-cl/v*tau))))*exp(-cl/v*ts)     # iv infusion concentration equation at steady state

  cl~dnorm(theta1, 100)
  v~dnorm(theta2, 1)

  theta1<-(0.01+0.0024*(293-2.03*age)*(0.5175-0.01685*Scr)/Scr*(0.86+0.14*Gender))*v      # population estimate of Aminoglycoside clearance
  theta2<-0.25*(0.4*bw+0.6*((0.73*Ht-59.42)*Gender+(0.65*Ht-50.74)*(1-Gender)))           # population estimate of Aminoclycoside apparent volume of distribution 
  }

  }