
model {
  for (i in 1:N) {
    killprob[i] <- c*((TBL[i]/d)*exp(1-TBL[i]/d))^g
    Kill[i] ~ dbin(killprob[i],10)
  }
  ## priors
  c ~ dunif(0,1)
  d ~ dnorm(0,0.001)I(0,)
  g ~ dnorm(0,0.001)I(0,)
}
