model

    {
    
    for (i in 1:1){           # one subject
    c~dnorm(mu,1.0E+6)        # c:drug concentration
    mu<-1.2*D/(v_F*(1.2-k))*((1/(1-exp(-k*tau)))*exp(-k*ts)-(1/(1-exp(-1.2*tau)))*exp(-1.2*ts))     # oral concentration equation at steady state
    k<-cl_F/v_F

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

    theta1<-0.0734*pow((D/TBW/tau*24), 0.406)*pow(TBW, 0.694)*pow(1.45,PHT) *pow(1.17,PB)*pow(1.21,VPA)*pow(0.849,E)   # population estimate of Carbamazepine clearance
    theta2<-1.91*TBW                                                                                                   # population estimate of Carbamazepine apparent volume of distribution 
       
    }

    }