model

           {
     
           for (i in 1:1){                 # one subject
           for(j in 1:T) {                 # T:number of sampling time
           c[j]~dnorm(mu[j],1.0E+6)        # c:drug concentration
           mu[j]<-((D/(1.5*cl_F))*(1-exp(-k*1.5))*(exp(-k*(ts[j]-1.5)))*((exp(-k*n[j]*tau))-1))/((exp(-k*tau))-1)    # oral concentration equation at non-steady state
           }
           k<-cl_F/v_F

           cl_F~dnorm(theta1, 500)
           v_F~dnorm(theta2, 1)

           theta1<-(13.8-3.81*OCC)*pow((BW/80),0.301)*pow((Hb/13),0.897)*pow((WBC/16),(-0.105))         # population estimate of Imatinib mesylate clearance
           theta2<-(252-7.82*OCC)*pow((BW/80),0.405)*pow((Hb/13),0.676)*pow((WBC/16),0.07)              # population estimate of Imatinib mesylate apparent volume of distribution 
              
           }

           }
           
           
          
          