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]<-(D*(1-exp(-cl/v*tin))/(tin*cl*(1-exp(-cl/v*tau))))*exp(-cl/v*ts[j]) # oral concentration equation at steady state
    }
        cl~dnorm(theta1,10)
        v~dnorm(theta2, 0.1)

        theta1<-1.08*(140-age)*bw/(72*Scr)*pow(0.85,(1-Gender))*60/1000  # population estimate of Vancomycin clearance
        theta2<-0.98*bw                                                  # population estimate of Vancomycin apparent volume of distribution

    }
    }
