model

           {
           for (i in 1:1) {
           INR~dnorm(mu,1.0E+6)   ## INR~dnorm(mu,1.0E+6)
           mu<-pow(a,(1/0.383))
           a<-((1/((-m*(cl_F/v_F))/(pow(kc,2))*(1-(kc*Dtau/(1-exp(-kc*Dtau))))
               -m/kc*log((D/v_F)/(Cpmax*(1-exp(-cl_F/v_F*Dtau)))))+3.36)/4.368)
            }
           m~dgamma(0.422,0.1)      ## m~dgamma(0.1,0.1)    
           Cpmax~dgamma(5.49,0.1)  ## Cpmax~dgamma(0.1,0.1)
           kc~dgamma(0.98,0.1)     ## kc~dgamma(0.1,0.1)   
           cl_F~dgamma(2.19,0.1)   ## cl_F~dgamma(0.1,0.1) 
           v_F~dnorm(7.5,100)     ## v_F~dnorm(7.5,100)   
           }                      