model{

for (i in 1:3) {
    for (j in 1:229) {
    Pos_Total[i,j] ~ dbern(p_Pos[i,j])
    p_PosFinal[i,j] <- (Incidence[i] * Ep * N[i,j] * D/365)/(Et[i,j] * 100000) 
    p_Pos[i,j] <- max(0.0001, min(0.9999, p_PosFinal[i,j]))
    }
}
    for (i in 1:3) {
        Incidence[i] ~ dgamma(1.2, 0.002)
    }
    
  Ep <- 0.54
  D <- 259
}