model

 {
 
 for (i in 1:1){                                                              # one subject
 for(j in 1:T) {                                                              # T:number of sampling time (N)
 c[j]~dnorm(mu[j], 100)                                                       # c:drug concentration
 mu[j]<-(D*(1-exp(-cl/v*tin))/(tin*cl*(1-exp(-cl/v*tau))))*exp(-cl/v*ts[j])   # 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 
 }

 }
