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 match bounds on MLE fits
  c ~ dunif(0,1)
  d ~ dunif(0,50) ## dnorm(0,0.001)I(0,)
  g ~ dunif(-1,40) ## dnorm(0,0.001)I(0,)
}