#logit model for dose-response

model{
for(i in 1:k){
x[i] ~ dbin(theta[i], n[i])
logit(theta[i]) <- beta[1]+ beta[2]*dose[i]
}

beta[1] ~ dnorm(0, 0.01)
beta[2] ~ dnorm(0, 0.1)I(0,)

#MTD <- -(log(2)+beta0)/beta1	 
	}