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]  * Utilisation[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])) 
    
    }
}

    # Priors
    Utilisation[1] ~dbeta(4.965667, 38.23564)   # Adults
    Utilisation[2] ~dbeta(4.183256, 42.99697)   # Children
    Utilisation[3] ~dbeta(4.743089, 42.25244)   # Total

    for (i in 1:3) {
        Incidence[i] ~ dgamma(1.2, 0.002)
    }
    

  Ep <- 0.54
  D <- 259
}