model

           {     
           for (i in 1:1) {            # one subject
           for(j in 1:T) {             # T:number of sampling time
           c[j]~dnorm(mu[j], 100)   # c:drug concentration
           mu[j]<-1.85*0.85*D/(v_F[X]*(1.85-k))*((1/(1-exp(-k*tau)))*exp(-k*ts[j])-(1/(1-exp(-1.85*tau)))*exp(-1.85*ts[j]))   # oral IR tablet concentration equation at steady state

           }
           cl_F[X]~dnorm(theta1, 100)
           v_F[X]~dnorm(theta2, 0.1)
           k<-cl_F[X]/v_F[X]

           theta1<-(0.037*IBW-0.006*age)*pow(1.284,smoke)*pow(0.751,CHF)      # population estimate of Aminophylline anhydrous clearance
           theta2<-0.42*IBW                                                   # population estimate of Aminophylline anhydrous apparent volume of distribution 

           IBW<-Gender*(50+(2.3*(ht/2.4-60)))+(1-Gender)*(45+(2.3*(ht/2.4-60)))   

           }
           }
