model

           {
           for (i in 1:1) {
           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.1,0.1)
           Cpmax~dgamma(0.1,0.1)
           kc~dgamma(0.1,0.1)
           cl_F~dgamma(0.1,0.1)
           v_F~dnorm(7.5,100)
           }