model

           {
     
           for (i in 1:1) {      # one subject
           for(j in 1:T) {       # T:number of sampling time
           c[j]~dnorm(mu[j],100)
           mu[j]<-ka*D/(vd_F*(ka-k))*((1/(1-exp(-k*tau)))*exp(-k*ts[j])-(1/(1-exp(-ka*tau)))*exp(-ka*ts[j]))     # oral concentration equation at steady state
           }
           # k<-cl/v
           # v<-0.24*TBW
           k<-cl_F/vd_F          # added by YJ 

           # cl~dnorm(theta1,10000)
           cl_F~dnorm(theta1,1.)
           vd_F~dnorm(theta2,100)         # added by YJ 
           ka  ~dnorm(theta3,100)         # added by YJ 

           theta1<-0.012*pow(TBW,0.715)*pow((D/TBW*24/tau),0.306)*(1+0.359*CBZ)     # population estimate of children valproate clearance
           theta2<-0.24*TBW               # added by YJ
           theta3<-4                      # ka = 4 for valproate oral solution; added by YJ
              
           }

           }