model{
        for(i in 1:N){
          sigma ~ dnorm(0, tau)
                tmp1 <- 1-exp(-1*(ks+lamda_1s*sigma_s)*1)
                tmp2 <- ks+sigma_s*lamda_1s
                tmp3 <- kx*1
                C0 <- (-1)*(ks*us-sigma_s*lamda_0s)*(1-(1-exp(-1*tmp2))/tmp2)/(1-exp(-1*tmp2))+0.5*sigma_s*sigma_s/((1-exp(-1*tmp2))*(1-exp(-1*tmp2)))*(1-2*(1-exp(-1*tmp2))/tmp2+(1-exp(-2*tmp2))/(2*tmp2))+ (sigma_x*sigma_x*h*h*(1-lamda_t)*(1-lamda_t))/(2*kx*kx)*(1/3-((1-2*(((1-exp(-1*tmp3))/(kx*kx))-(1*exp(-1*tmp3)/kx)))/(1-exp(-1*kx)))+((1-2*((1-exp(-1*tmp3))/kx)+((1-exp(-2*tmp3))/(2*kx)))/((1-exp(-1*kx))*(1-exp(-1*kx)))))
                b1 <- 1
                b2 <- tmp1/(1 - exp(-1*tmp2))
                b3 <- (h*(1-lamda_j)/kx) * (1 - ((1-exp(-1*tmp3))/(1-exp(-1*kx))))
                a <- ((ks*us-sigma_s*lamda_0s)*(1-tmp1/tmp2))/ (1-exp(-1*tmp2)) - (((1/2*sigma_s*sigma_s)/((1-exp(-1*tmp2))*(1-exp(-1*tmp2)))*(1 - 2*tmp1/tmp2 + tmp1/(2*tmp2)))) + (sigma_x*sigma_x*h*h*(1-lamda_t)*(1-lamda_t))/(2*kx*kx)*(1/3*1*1 - ((1*1-2*(((1-exp(-1*tmp3))/(kx*kx))-(1*exp(-1*tmp3)/kx)))/(1-exp(-1*kx))) +((1-2*((1-exp(-1*tmp3))/kx)+((1-exp(-2*tmp3))/(2*kx)))/((1-exp(-1*kx))*(1-exp(-1*kx)))))+C0*1
                y <- a/1+b1/1*rt+b2/1*st+b3/1*xt+sigma
        }
        ks1 ~ dnorm(0.2, 1)
        kx1~ dnorm(0.2, 1)
        us ~ dnorm(0.2, 4)
        lamda_0s ~ dnorm(0, 1)
     lamda_1s ~ dnorm(0, 1)
     lamda_j ~ dnorm(0, 1)
        w1 ~ dnorm(-10, 100)
        w2 ~ dgamma(0.1, 0.1)
        w3 ~ dgamma(0.1, 0.1)
        w4 ~ dgamma(0.1, 0.1)
        h <- 12*exp(w1)/(1+exp(w1))
        ks <- ks1*step(ks1)
        kx <- kx1*step(kx1)
        tau <- 1/w3
        sigma_s <- sqrt(1/w2)
        sigma_x <- sqrt(1/w4)
        sigma_e <- sqrt(tau)
}
list(
rt  = c(1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 1.98, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.25, 2.52, 2.52, 2.52, 2.52, 2.52, 2.52, 2.52, 2.79, 2.79, 3.06, 3.06, 3.33, 3.60, 3.87, 3.87, 3.87, 4.14, 4.14, 4.14, 4.14, 4.14, 4.14, 4.14, 4.14, 4.14, 4.14, 3.74, 2.52, 2.25),
st  = c(0.2304, -0.0107, 0.0014, -0.0583, -0.1866, -0.1891, -0.1385, -0.1313, 0.0580, -0.0653, -0.0411, -0.2056, -0.4754, -0.2451, -0.0550, 0.0624, 0.1942, 0.2158, 0.1467, 0.3681, 0.5161, 0.5513, 0.4180, 0.3560, 0.4876, 0.5460, 0.6929, 0.9877, 0.8995, 0.9905, 1.0134, 1.0507, 0.7504, 0.7610, 0.6639, 0.5703, 0.5872, 0.2145, -0.0217, -0.0684, -0.1506, -0.5164, -0.6565, -0.8284, -0.8292, -0.6157, -0.4141, -0.6072, -0.5636, -0.5510, -0.5374, -0.5218, -0.3282, -0.1958, -0.4018, -0.5050, -0.4913, -0.4462, -0.4149, -0.4523, -0.3810, -0.6142, -0.5004, -0.6855, -0.4840, -0.4877, -0.7030, -0.8081, -0.6037, -0.2616, -0.4948, -0.6449, -0.7989, -0.9342, -0.9595, -0.8556, -0.6973, -0.6030, -0.5810, -0.8152, -0.9108, -0.2559, -0.8939),
xt  = c(0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.27, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.27, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.27, 0.00, 0.27, 0.00, 0.27, 0.27, 0.27, 0.00, 0.00, 0.27, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, -0.41, -1.22, -0.27, -2.25),
y = c(2.21, 1.97, 1.98, 1.92, 1.79, 1.79, 1.84, 1.85, 2.04, 1.91, 1.94, 1.77, 1.50, 1.73, 1.93, 2.04, 2.17, 2.20, 2.13, 2.35, 2.50, 2.53, 2.40, 2.34, 2.47, 2.53, 2.67, 2.97, 2.88, 2.97, 2.99, 3.03, 3.00, 3.01, 2.91, 2.82, 2.84, 2.46, 2.23, 2.18, 2.10, 1.73, 1.59, 1.42, 1.42, 1.63, 1.84, 1.64, 1.69, 1.70, 1.71, 1.73, 1.92, 2.05, 2.12, 2.02, 2.03, 2.07, 2.11, 2.07, 2.14, 2.18, 2.29, 2.37, 2.58, 2.84, 2.90, 3.06, 3.27, 3.61, 3.65, 3.50, 3.34, 3.21, 3.18, 3.28, 3.44, 3.54, 3.56, 3.32, 2.82, 2.26, 1.36),
N = 83)